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SECTION  1 
INTRODUCTION 


1.1  BACKGROUND 

For  more  than  two  decades  it  has  been  well  known  that  propaga¬ 
tion  through  the  natural  ionosphere  can  cause  degradation  of  satellite 
signals  (Yeh  and  Swenson,  1959;  Skinner,  et  al.,  1971;  Pope  and  Fritz, 
1971;  Taur,  1976;  Whitney  and  Basu,  1977;  Johnson  and  Lee,  1978;  Mullen, 
et  al.,  1978;  Paulson  and  Hopkins,  1978;  Fremouw,  et  al.,  1978).  It  has 
also  been  widely  reported  for  a  similar  period  of  time  that  the  ionosphere 
is  subject  to  even  greater  disturbances  following  high  altitude  nuclear 
explosions  (Matsushita,  1959;  Keyes  and  Tinsley,  1962;  Rothwell,  et  al., 
1963;  Maeda,  et  al.,  1964;  Arendt  and  Soicher,  1964;  Zinn,  et  al.,  1966; 
Boquist  and  Snyder,  1967;  Hoerlin,  1976;  Glasstone  and  Dolan,  1977;  King 
and  Fleming,  1980).  In  both  natural  and  artificially  disturbed  iono¬ 
spheres,  electron  density  irregularities  can  cause  random  variations  in 
the  amplitude  and  phase  of  a  propagating  wave.  Such  signal  variations  are 
called  scintillation  or  fading  and  have  been  observed  on  satellite  links 
through  the  natural  ionosphere  at  frequencies  up  to  the  GHz  range  (Taur, 
1976;  Fremouw,  et  al.  1978). 

If  all  frequency  components  of  the  received  scintillating  signal 
vary  essentially  identically  with  time,  the  propagation  channel  is  re¬ 
ferred  to  as  nonselect ive  or  as  a  flat  fading  channel.  When  the  scintil¬ 
lations  exhibit  statistical  decorrelation  at  different  frequencies  within 
the  signal  bandwidth,  the  channel  is  referred  to  as  frequency  selective. 

Frequency  selective  scintillations  are  therefore  encountered 
when  the  bandwidth  utilized  by  a  communication  link  exceeds  the  coherence 
bandwidth  of  the  ionospheric  channel.  This  situation  is  more  likely  to 
occur  when  the  system  bandwidth  is  large,  as  in  spread  spectrum  systems. 
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Such  systems  typically  emoloy  either  direct  sequence  spectrum  spreading 
over  a  large  instantaneous  bandwidth,  or  carrier  frequency  hopping  of  a 
small  instantaneous  signal  spectrum  over  a  large  bandwidth.  The  effects 
of  frequency  selective  scintillation  are  quite  different  for  direct 
sequence  and  frequency  hopped  systems  (Bogusch,  et  al.,  1981).  Only  the 
case  of  an  instantaneous  wide  bandwidth  signal  is  considered  in  this  work 
although  the  appropriate  results  presented  here  may  easily  be  applied  to 
the  latter  case. 

Frequency  selective  scintillation  is  much  more  likely  to  occur 
if  the  ionosphere  is  more  highly  disturbed,  as  for  example  by  nuclear 
explosions  or  by  chemical  releases  (Knepp  and  Bogusch,  1979).  Increased 
electron  concentrations  and  irregular  structure  of  the  ionization  can  lead 
to  intense  signal  scintillations  and  potentially  significant  frequency 
selective  effects  at  frequencies  as  high  as  the  7-8  GHz  SHF  band  (Knepp, 
1977).  Consequently,  the  effects  of  frequency  selective  propagation  dis¬ 
turbances  are  of  importance  to  spread  specLrum  satellite  systems  that  may 
have  to  operate  in  highly  disturbed  environments. 

In  this  report,  several  techniques  are  presented  to  determine 
the  effects  of  ionization  irregularities  on  a  propagating  wide  bandwidth 
signal  subject  to  frequency  selective  scintillation. 

1.2  MRS  PROPAGATIOM  SIMULATION 

The  multiple  phase-screen  (MPS)  propagation  simulation  (Knepp, 
1977;  Wittwer,  1978)  is  the  most  general  technique  to  obtain  a  solution  to 
the  parabolic  wave  equation.  In  this  solution  technique,  the  ionized 
mediim  is  divided  into  a  finite  number  of  layers.  The  field  Fluctuation 
through  each  layer  is  calculated  by  replacing  the  layer  by  a  phase¬ 
changing  screen  located  at  its  center,  whose  statistical  properties  are 
determined  from  the  statistics  of  the  electron  density  fluctuations  within 
Lhe  layer.  At  each  phase-screen,  the  statistical  phase  is  added  to  the 
electric  field  phase  and  the  wave  is  then  propagated  to  the  next  phase- 
"Screen  or  to  the  observer.  For  wide  bandwidth  waveforms,  the  numerical 
solution  is  obtained  at  a  number  of  discrete  frequencies  centered  about  a 
carrier  frequency  and  then  a  time-domain  solution  is  obtained  by  the  use 
of  Fourier  transform  techniques. 
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The  MPS  propagation  technique  is  quite  general,  and  may  be 
easily  applied  to  problems  involving  numerous,  separated,  layers  of  ioni¬ 
zation  characterized  by  spatially  varying  electron  density  power  spectra. 
Section  2  of  this  report  describes  the  MPS  numerical  propagation  simula¬ 
tion  and  its  application  to  wide  bandwidth  signal  generation. 

The  MPS  propagation  simulation  can  handle  all  levels  of  iono¬ 
spheric  disturbances  from  the  least  severe,  where  only  minor  phase  fluc¬ 
tuations  occur,  to  the  most  severe  cases  of  frequency  selective  Rayleigh 
fading.  Since  the  MPS  simulation  provides  a  direct  solution  to  the 
parabolic  wave  equation,  the  results  are  exact  given  a  certain  description 
of  the' propagation  environment.  Thus  the  applicability  of  the  MPS  propa¬ 
gation  simulation  depends  upon  the  accuracy  of  the  description  of  the 
scattering  medium.  In  most  of  the  MPS  results  presented  in  this  report, 
the  propagation  environment  is  described  by  a  one-dimensional  K"^ 
power-law  phase  power  spectrum.  This  spectrum  corresponds  to  a 
three-dimensional  power  spectrum  for  electron  density  fluctuations  and  is 
representative  of  many  in-situ  measurements  (Dyson,  et  al,  1974;  Phelps 
and  Sagalyn,  1976)  as  well  as  numerical  simulation  of  striation  spectra 
(Scannapieco,  et  al.,  1976). 

Section  2  of  this  report  describes  the  MPS  propagation  simula¬ 
tion  in  detail  as  it  is  applied  to  calculate  realizations  of  received 
time-domain  signals. 

The  word  "realization"  is  used  in  this  report  to  refer  to  a 
specific  sample  of  a  statistical  process.  For  example,  realizations  of 
the  phase-changing  screen  are  generated  in  a  random  manner  based  upon  a 
sequence  of  pseudo-random  numbers  generated  numerically.  However,  once 
these  "random"  numbers  are  generated,  the  phase-screen  is  completely 
defined  and  known.  Similarly,  the  wave  field  propagated  through  this 
phase-screen  is  subject  to  exact  calculation.  However,  a  different 
sequence  of  pseudo-random  nunbers  will  yield  a  different  phase-screen 
realization  and  hence  a  different  realization  of  the  received  electric 
field . 
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An  important  function  of  the  MPS  propagation  simulation  is  to 
serve  as  an  intermediate  step  in  the  analysis  of  the  effects  on  receivers 
of  a  disturbed,  ionized  propagation  environment.  Thus  realizations  of  the 
received  signal  amplitude  and  phase  are  important  results  of  the  simula¬ 
tion  and  in  turn  serve  as  direct  input  to  detailed  receiver  simulations 
(Bogusch,  ot  al.,  1981;  Knepp  and  Bogusch,  1979). 

In  Section  2  the  MPS  propagation  simulation  is  applied  to  pre¬ 
dict  the  effect  of  a  large,  striated  barium  cloud  on  a  wide  bandwidth 
signal  passing  through  the  disturbed  ionization.  These  predictions  are 
based  on  an  analysis  of  several  earlier  bariim  releases  (Prettie  and 
Marshall,  1978)  from  which  a  model  of  barium  cloud  irregularity  structure 
was  obtained.  The  barium  cloud  is  then  represented  as  a  single  phase- 
screen  in  the  MPS  simulation.  A  large  Gaussian  phase  lens  is  used  to 
model  the  mean  or  deterministic  ionization  effects  from  the  barium  cloud 
and  small  level  stochastic  phase  perturbations  are  added  to  represent  the 
ionization  fluctuations.  This  model  is  quite  general  and  has  been  applied 
to  describe  several  other  observed  barium  clouds  (Knepp  and  Bogusch, 

1979).  Results  for  the  received  wide  bandwidth  signal  after  propagation 
through  the  barium  cloud  are  presented  and  analyzed.  Scattering  from 
various  components  of  the  barium  cloud  is  considered  independently  and 
then  in  combination  to  give  a  useful  and  intuitive  understanding  of  the 
elements  of  the  scattering  process  that  are  important  for  wide  bandwidth 
signals. 


The  Defense  Nuclear  Agency  barium  release  experiment  discussed 
above  took  place  in  December  1980  and  the  results  for  the  wide  bandwidth 
propagation  experiment  were  provided  by  Dr.  James  Marshall  (Marshall, 
1982).  The  agreement  between  the  1979  predictions  and  the  1980  experiment 
is  quite  remarkable  and  is  discussed  in  Section  2. 

Major  emphasis  in  Section  2  is  also  given  to  a  comparison  of  the 
MPS  results  to  theoretical  results.  As  stated  earlier,  the  MPS  propaga¬ 
tion  simulation  is  the  only  effective  method  to  provide  signal  realiza¬ 
tions  applicable  to  a  wide  range  of  propagation  conditions.  For  compari¬ 
son  purposes,  these  signal  realizations  may  be  analyzed  to  obtain  average 
quantities  for  which  analytic  solutions  or  approximations  can  be  found. 


18 


MPS  results  for  the  two-position,  two-frequency  mutual  coherence 
function  r(Ax,Au))  are  compared  to  a  theoretical,  strong  scatter  approxima¬ 
tion  developed  in  Section  4.  Although  the  theoretical  calculation  is  more 
restrictive  than  the  MPS  propagation  technique,  good  agreement  between 
numerical  and  theoretical  results  is  shown.  As  an  additional  comparison, 
the  time  domain  MPS  results  for  the  mean  time  delay  and  time  delay  jitter 
of  a  wide  bandwidth  PN  spread  spectrum  signal  are  compared  to  theoretical 
moment  method  calculations  (Yeh  and  Liu,  1977). 

1.3  STATISTICAL  SIGNAL  GENERATION 

In  Section  3  a  second  analytical/numerical  technique  is  describ¬ 
ed  to  generate  realizations  of  the  received  signal  after  propagation  of  a 
wide  bandwidth  waveform  through  a  layer  of  strongly  turbulent  media. 

These  statistical  signal  realizations  are  generated  to  have  Rayleigh  amp¬ 
litude  statistics  (Fante,  1975)  and  to  have  spatial  and  frequency  correla¬ 
tion  properties  which  obey  the  parabolic  wave  equation  in  the  strong  scat¬ 
ter  limit.  Although  these  statistical  signal  realizations  apply  only  to 
the  case  of  strong  scattering,  their  generation  requires  only  a  fraction 
of  the  computer  resources  required  for  signal  generation  by  wide  bandwidth 
multiple  phase-screen  calculations.  Furthermore,  the  limitation  to  the 
strong  scatter  regime  is  not  of  great  significance  since  frequency  selec¬ 
tive  scintillation  only  occurs  during  strong  scattering  conditions. 

This  statistical  technique  follows  the  basic  formalism  described 
in  Wittwer  (1980)  and  is  based  upon  the  solution  for  the  two-frequency 
mutual  coherence  function,  r,  for  spherical  wave  propagation  with  trans¬ 
mitter  and  receiver  located  on  opposite  sides  of  a  thick  finite  layer  of 
ionized  electron  density  irregularities.  An  analytic  solution  is  obtained 
for  r  in  the  strong  scatter  regime  by  use  of  the  quadratic  approximation 
for  the  phase  structure-function.  The  thin  phase-screen  approximation  to 
the  thick  layer  is  then  utilized  and  great  simplification  to  the  analytic 
expression  for  r  is  obtained.  The  relationship  between  the  impulse 
response  function  of  the  propagation  channel  and  the  two-frequency  mutual 
coherence  function  and  its  Fourier  transform  is  then  used  to  directly 
obtain  statistical  realizations  of  wide  bandwidth  signals. 
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The  accuracy  of  the  thin  phase-screen  approximation  is  discussed 
and  several  comparisons  of  statistical  signals  with  signals  obtained 
directly  from  the  MPS  propagation  simulation  are  shown  in  Section  3. 

The  analytic  solution  for  r  which  is  the  foundation  of  this 
statistical  technique  to  generate  realizations  is  presented  in  Section  4. 

1.4  SOLUTION  FOR  TWO-FREQUENCY  MUTUAL  COHERENCE  FUNCTION 

An  analytic  solution  is  obtained  in  Section  4  for  the  two-posi¬ 
tion,  two-frequency  mutual  coherence  function  for  spherical  wave  propaga¬ 
tion.  This  solution  provides  the  basis  of  the  statistical  signal  genera¬ 
tion  technique  presented  in  Section  3  of  this  report.  It  is  assumed  that 
strong  scattering  conditions  prevail  and  tliat  the  quadratic  approximation 
to  the  phase  structure-function  is  valid.  Here  the  ionized  scattering 
region  occupies  a  thick,  finite  layer  with  transmitter  and  receiver 
located  in  free  space  on  opposite  sides  of  the  layer.  General  analytic 
solutions  are  derived  for  two  cases.  In  the  first  case  the  random  slab  is 
represented  by  a  one-dimensional  power  spectrum  of  electron-density  fluc¬ 
tuations  corresponding  to  propagation  through  elongated  irregularities  as 
would  occur  for  an  equatorial  satellite  link  to  a  ground  station.  In  the 
second  case  the  random  slab  is  represented  by  isotropic  ionization  irregu¬ 
larities  which  corresponds  to  the  physical  geometry  of  propagation  along 
the  direction  of  the  earth's  magnetic  field. 

The  case  of  isotropic  irregularities  represents  a  worst  case 
condition  while  the  case  of  infinitely  elongated  irregularities  leads  to 
less  severe  propagation  conditions.  Both  cases  taken  together  represent 
the  extremes  of  the  range  of  results  to  be  expected  for  propagation 
through  ionospheric  fluctuations,  or  ionization  ir»*egularities  caused  by 
bariun  cloud  instabilities  or  by  nuclear  detonations. 

For  both  cases  the  complex  general  analytic  results  are  simpli¬ 
fied  by  use  of  the  tnin  phase-screen  approximation  (zero  slab  thickness) 
to  obtain  useful  analytic  expressions  for  r  as  well  as  the  resulting  power 
impulse  response  function.  It  is  shown  that  the  impulse  response  to  an 
incident  power  delta  function  reduces  to  an  exponential  form  in  the  strong 
scattering  limit  and  to  a  Gaussian  form  in  the  weak  scattering  limit.  The 
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Gaussian  form  corresponds  to  pulse  wander  and  dispersion  while  the  expo¬ 
nential  form  corresponds  to  diffractive  pulse  spreading  produced  by  multi- 
path  effects. 

Since  the  thin  phase-screen  results  developed  in  Section  4  are 
used  in  the  statistical  signal  generation  technique  described  in  Section 
3,  the  accuracy  of  the  thin-screen  approximation  is  a  matter  of  some 
importance.  This  approximation  is  considered  in  some  detail  in  Section  4 
where  comparisons  of  the  thick  layer  cases  to  the  thin  phase-screen  case 
are  presented  for  the  S^  scintillation  index,  the  mean  time  delay,  <t>, 
and  the  time  delay  jitter,  <j^.  It  is  found  that  the  thin  phase-screen 
approximation  gives  exact  results  for  <t>  and  accurate  results  for  S^  and 


The  simplification  resulting  in  the  thin  phase-screen  approxima- 
I  ion  for  tlie  two-frequency  mutual  coherence  function  for  spherical  waves 
^Iso  leads  to  simplified  expressions  for  two  very  important  parameters 
that  describe  the  propagation  medium.  These  parameters  are  the  signal 
decorrel ation  distance  (and  time)  and  the  coherence  bandwidth.  Both 
parameters  are  defined  and  discussed  in  Section  4  and  it  is  shown  that  the 
signal  decorrelation  time  and  the  coherence  bandwidth  both  obey  the  prin¬ 
ciple  of  reciprocity. 
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SECTION  2 

NPS  NUMERICAL  PROPAGATION  SIMULATION 


2.1  PHASE-SCREEN  TECHNIQUES 

Random  variations  in  the  amplitude  and  phase  of  a  propagating 
wave,  called  scintillations,  are  caused  by  propagation  through  a  region  in 
which  the  electron  density  is  irregularly  structured.  In  ionospheric 
plasmas,  electron  density  structure  is  produced  by  plasma  instabilities 
which  cause  the  ionization  to  break  up  into  long  filaments,  or  striations, 
aligned  with  the  earth's  magnetic  field  lines  (Linson  and  Workman,  1970). 
The  instabil ity  mechanisms  accounting  for  ionospheric  irregularities  have 
been  the  topic  of  debate  for  many  years  and  are  summarized  in  a  recent 
review  (Basu  and  Kelley,  1978)  .  Roughly  speaking,  one  may  describe  stria¬ 
tions  as  long  sheets  or  rods  of  relatively  high  electron  density  imbedded 
in  a  background  of  lower  electron  density. 

Figure  2-1  illustrates  the  geometry  of  the  problem.  Propagation 
of  satellite  signals  through  a  large  striated  region  presents  the  problem 
of  radio  wave  propagation  through  a  thick  medium  composed  of  random  fluc¬ 
tuations  in  the  index-of-refraction.  Consider  for  a  moment  a  plane, 
unmodulated  carrier  wave  traversing  the  striated  region.  The  wave  first 
suffers  random  phase  perturbations  due  to  variations  in  the  phase  velocity 
within  the  medium.  These  phase  variations  in  the  propagating  wavefront 
introduce  small  random  changes  in  the  direction  of  propagation  of  the 
wave.  Thus  portions  of  the  once  plane  wavefront  now  propagate  in  different 
directions  relative  to  other  portions.  As  the  wave  propagates  farther, 
diffraction  or  angular  scattering  causes  constructive  and  destructive 
interference  which  introduces  fluctuations  in  amplitude  as  well  as  phase. 
These  time  varying  amplitude  and  phase  fluctuations  represent  an  undesired 
complex  modulation  of  the  carrier. 

The  actual  satellite  signal  encompasses  a  spectrim  of  frequen¬ 
cies  because  of  the  transmitted  modulation.  This  transmitted  signal  can 
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be  Fourier  decomposed,  and  the  propagation  of  each  spectral  component  can 
be  analyzed.  Because  the  phase  velocity  in  an  ionized  medium  dPoends  on 
wave  frequency,  each  spectral  component  experiences  somewhat  different 
phase  perturbations  and  hence  exhibits  different  amplitude  and  phase  scin¬ 
tillations  after  propagating  to  the  receiver.  Thus  the  received  signal 
spectrun  is  distorted  by  both  the  temporal  (time  selective)  and  dispersive 
(frequency  selective)  properties  of  the  ionospheric  channel. 

Radio-wave  propagation  through  an  extended  region  of  irregular¬ 
ities  has  been  studied  for  many  years  by  numerous  researchers  (Bramley, 
1954;  Ratcliffe,  1956;  Mercier,  1962;  Tatarskii,  1967  and  1971;  Ishimaru, 
1978  and  references  therein).  Although  partial  solutions  are  available  in 
both  weak  (Tatarskii,  1967)  and  strong  scattering  cases  (Fante,  1975  and 
1981),  no  general  analytical  solution  has  been  obtained  even  for  a  single 
frequency,  much  less  for  a  wide-band  signal.  However,  the  problem  may  be 
solved  by  the  use  of  an  analytical/numerical  technique  known  as  the 
multiple  phase-screen  (MRS)  method  (Knepp,  1977;  Wittwer,  1978), 

The  use  of  phase-screen  techniques  to  calculate  the  propagation 
of  electromagnetic  waves  was  formulated  as  early  as  1956  by  Ratcliffe.  A 
number  of  researchers  have  used  phase-screen  techniques  in  the  solution  of 
problems  involving  propagation  through  random  media  (Mercier,  1961; 
Salpetei",  1967;  Jokippi,  1970;  Buckley,  1975;  Rino,  1980).  Recently, 
pliase-screen  techniques  have  been  used  in  the  simulation  of  adaptive 
optical  systems  (Brown,  1975)  and  in  the  propagation  of  high  energy  laser 
beams  subject  to  thermal  blooming  (Fleck,  et  al.,  1976).  The  use  of 
phase-screen  techniques  to  calculate  the  propagation  of  wide  bandwidth 
signals  is  new. 

The  analytical  foundation  of  the  MRS  technique  was  neatly  set 
forth  by  Ratcliffe  in  1956.  The  striated  ionosphere  is  represented  by  a 
series  of  thin  diffracting  screens  at  various  points  along  the  propagation 
path.  Each  screen  is  perpendicular  to  tne  path  and  consists  of  a  random 
phase  equal  to  the  integrated  phase  change  due  to  electron  content  within 
the  local  region  represented  by  the  screen.  An  incident  electromagnetic 
wave  emerges  from  the  screen  with  this  random  phase  shift  superimposed. 
Thus,  at  any  position  along  the  path,  the  wave  can  be  represented  by  a 
complex  field  whose  amplitude  and  phase  are  functions  of  position  in  a 
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plane  perpendicular  io  the  nominal  direction  of  propagation.  The  rippled 
wave  field  can  be  thought  of  as  a  collection  of  uniform  plane  waves 
traveling  in  slightly  different  directions.  This  collection  of  plane 
waves,  called  the  angular  spectrin  of  the  wave  field,  is  given  by  the 
spatial  Fourier  transform  of  the  field  distribution  over  the  diffracting 
screen.  The  angular  spectrin  can  be  propagated  to  the  next  diffracting 
screen  by  introducing  the  free-space  phase-shift.  Transforming  back  to 
spatial  coordinates  yields  the  incident  electric  field  distribution,  which 
is  then  perturbed  by  the  random  phase-shift  in  this  screen.  The  modified 
angular  spectrin  is  obtained  and  propagated  to  the  next  screen,  and  the 
process  is  repeated  until  the  field  distribution  is  obtained  at  the 
receiver  plane. 

In  this  section  the  MPS  propagation  simulation  and  its  limita¬ 
tions  are  described  with  some  results  included  from  an  earlier  report 
(Knepp,  1977).  Time-domain  results  after  propagation  through  a  strong 
dispersive  Gaussian  lens  are  presented  and  shown  to  be  identical  to  analy¬ 
tic  results.  Also  shown  are  predictions  of  the  expected  dispersive  and 
scattering  effects  due  to  propagation  of  a  pseudo-random  phase-shift  keyed 
modulated  waveform  through  a  finite  sized  barium  cloud.  Measured  results 
from  a  later  experiment  have  been  provided  by  Dr.  James  Marshall  and  are 
shown  to  be  very  similar  to  the  pre-experimental  predictions. 

2.2  FORMULATION 

In  the  parabolic  approximation,  let  the  electric  field  be 

written 


E(x,z,u))  =  U(x,z,a))  exp(-ikz)  (2-1) 

so  that  the  parabolic  wave  equation  (Tatarskii,  1971)  becomes 

-  2ik  IH  +  2k2An(x,z,a))U  =  0  (2-2) 

3x2  3Z 

where  k  is  the  wavenumber  and  ^n  is  the  variation  in  the  index-of-refrac 
tion.  For  ionization  irregularities 
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(2-3) 


where  ANg  is  the  variation  in  the  electron  density,  n^  is  the  critical 
electron  density  and  rg  is  the  classical  electron  radius  (2.82x10"^^). 

A  derivation  of  Equation  2-2  is  presented  in  Section  3  for  the  more  gener¬ 
al  case  where  the  propagation  mediun  has  a  background  ionization  level. 

In  Equation  2-2  it  is  assumed  that  the  mean  or  background  ionization  level 
is  zero.  In  writing  Equation  2-2  it  is  also  assumed  that  the  z-direction 
is  the  direction  of  propagation.  The  striations  or  irregularities  are 
assumed  infinitely  elongated  in  the  y-direction  as  would  effectively  be 
the  case  for  field  aligned  ionization  and  a  satellite  communication  link 
between  a  synchronous  satellite  and  an  equatorial  ground  station.  These 
assumptions  limit  the  application  of  the  MPS  simulation  to  a  two-dimen¬ 
sional  geometry  with  no  variation  in  the  y-direction.  The  simulation 
could  easily  be  extended  to  handle  a  three-dimensional  geometry  but  at 
considerably  greater  cost  in  computer  resources. 

Now  assume  that  the  ionized  region  has  been  divided  into  a 
number  of  thin  layers,  with  each  layer  perpendicular  to  the  direction  of 
propagation.  Consider  a  layer  of  thickness  az  centered  at  zero  z.  For 
small  AZ,  the  equation  for  propagation  through  this  thin  layer  is  obtained 
from  Equation  2-2  with  the  first  term  neglected.  The  remaining  differen¬ 
tial  equation  is  easily  solved. 


U(x,_,u)  =  U(x,- _,(d)  exp  {-Ik  f  An(x,z' ,w)dz' ' 

2  2  ' 

-Az/2 


(2-4) 


The  exponential  quantity  is  simply  the  geometric  optics  phase  change 
imparted  aftT  propagation  through  the  layer.  Equation  2-4  states  that 
transversal  of  a  thin  layer  by  a  propagating  wave  is  accomplished  by 
simply  adding  the  phase  change  associated  with  the  layer  to  the  phase  of 
the  electric  field  at  the  entrance  to  the  layer. 

Now  collapse  the  irregularities  in  every  thin  layer  to  thin 
phase-screens  whose  effect  is  represented  mathematically  by  Equation  2-4. 


.2,8 


Now  the  ionized  medium  is  represented  by  a  series  of  thin  phase-screens 
separated  by  free-space.  Equation  2-2  is  valid  for  free-space  propagation 
in  the  region  between  the  phase-screens  if  the  last  term  of  the  equation 
is  neglected.  The  resulting  equation  is  easily  solved  using  Fourier 
transform  techniques.  The  solution  for  free-space  propagation  from  z^  to 
Z2  can  be  written 

U(K,Z2,u)=—  ?U{x,Z2,a))e''"^’‘dx 

=  0{K,Zi,a»)  exp[iK2(z2-Zi)/2k]  (2-5) 

U(x,Z2.a))  =  JU{K,Z2,0))e’'^V 

•00 

00 

*  J"  U(K,Zi,a))  exp[iK2(22-2i)/2k+iKx]dK  (2-6) 

Equation  2-1  may  be  inserted  into  Equation  2-6  to  write  the  equation 
governing  free-space  propagation  in  terms  of  the  electric  field 

E(x,Z2,a))  =  jfE(K,Zi,u))  exp [iK2(z2-Zi)/2k+iKx ]dK  (2-7) 


Equations  2-4  and  2-7  may  be  used  to  propagate  an  electromagnetic  wave 
from  the  transmitter  to  the  receiver  as  follows.  Assume  that  the  incident 
wave  is  a  plane  wave  as  it  approaches  the  first  phase-screen.  Immediately 
after  passing  through  this  phase-screen  only  the  phase  of  the  electric 
field  is  affected  according  to  Equation  2-4.  The  suitably  modified  elec¬ 
tric  field  is  then  Fourier  transformed  to  obtain  E(K,z,a))  at  the  exit  of 
the  first  phase-screen.  Free-space  propagation  to  the  next  phase-screen 
(located  at  Z2)  is  then  accomplished  according  to  Equation  2-7.  Equation 
2-7  is  implemented  numerically  in  two  steps;  first  a  multiplication  in 
K-space  by  exp(iK2(z2-Zj)/2k)  and  then  by  an  efficient  fast  Fourier  trans¬ 
form.  This  brings  the  wave  to  the  next  phase-screen  where  the  above 
process  is  repeated  until  the  receiver  plane  is  reached.  Thus  a  nimerical- 
ly  efficient  solution  to  the  parabolic  wave  equation  is  obtained. 
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In  effect  this  solution  is  equivalent  to  replacing  the  scatter¬ 
ing  mediim  by  a  number  of  diffracting  phase-screens  separated  by  free- 
space.  For  example,  if  three  phase-screens  are  used  to  represent  the 
geometry  of  Figure  2-1,  the  resulting  MPS  geometry  is  shown  in  Figure  2-2. 

This  formulation  in  two  dimensions  assumes  that  no  variations 
occur  in  the  third  or  y  dimension.  This  is  often  a  good  approximation  for 
ionospheric  propagation  because  the  striated  ionization  is  often  elongated 
parallel  to  the  earth's  magnetic  field  (Paulson  and  Hopkins,  1978;  Fremouw 
et  al.,  1978).  The  formulation  can  easily  be  extended  to  three  dimen¬ 
sions,  at  the  expense  of  a  significant  increase  in  numerical  computation. 

2.2.1  Phase  Power  Spectral  Density 

Specific  realizations  of  the  random  phase  ^{x)  that  define  a 
phase-screen  are  obtained  by  sampling  a  distribution  of  phase-shifts 
having  statistical  properties  determined  the  statistics  of  the  electron 
density  irregularities.  These  statistics  are  specified  by  the  spatial 
power  spectral  density  (PSD)  of  the  irregularities  or,  equivalently,  by 
their  spatial  autocorrelation  function. 

The  relationship  between  the  PSD  of  the  phase  and  the  PSD  of  the 
in-situ  electron  density  irregularities  is  obtained  as  follows.  From 
Equation  2-4  the  phase  that  defines  a  phase-screen  is  given  by 

Az/2 

't'(x)  =  -  k  f  An(x,z')dz'  (2-8) 

-az/2 

where  the  frequency  dependence  is  temporarily  omitted  for  convenience. 
Using  relation  2-3  the  phase  autocorrelation  function  can  be  formed 

B^(C)  =  <<>(x)(fr*(x+5)> 

Az/2  az/2 

"  /  /  <4N  (x,z)aN  (x+c,z')>  dzdz'  (2-9) 

"  -AZ/2  -az/2  ®  ® 
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The  quantity  within  the  angle  brackets  is  recognized  as  the  correlation 
function  of  electron  density  fluctuations  5,z-z' ).  The  double  inte¬ 
gral  here  can  be  reduced  to  a  single  integration  by  a  change  of  variables 
as  explained  in  Papoulis,  1965,  p.  325  with  the  result 

B^{?)  =  r^x^AZ  (l-|z'|/Az)B^^(C,z')dz'  (2-10) 


If  Az  is  greater  than  the  correlation  length  of  the  electron  density 
fluctuations,  the  contribution  of  the  second  term  in  the  integration  will 
always  be  negligible  so  that  the  limits  can  be  extended 


B^(?)  =  r2x2AZ 


f  (C,z')dz' 
-00  e 


(2-11) 


Evaluation  of  this  expression  at  5  =  0  gives  for  the  variance  of  the  phase 
fluctuations 


f  (2-12) 

"  e 

Multiplication  of  Equation  2-11  by  exp(-iK5)/2n  and  integration  with 
respect  to  c  gives 


I-  f  B  (Oe’^'^^dc  =  (r2x2Az/2TT)  jT  JTb^  ( C,z' )e"'''^^dz'd5 

cTT  ^  .QO  «.oo  0 


=  s  JK)  (2-13) 

The  term  on  the  left  hand  side  of  the  equation  is  recognized  as  the  phase 
PSD  and  thus  the  one-dimensional  phase  PSD  is  simply  related  to  the  inte¬ 
gral  of  the  two-dimensional  autocorrelation  function  of  in-situ  electron 
density  fluctuations.  Now  the  relationships  between  one-dimensional  PSD's 
and  two-dimensional  PSD's  and  autocorrelation  functions  are 
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(2-14) 


Bi(x)  =  jrSi(K)e^'^’‘dK 
Si(K)  /Bi(x)e"’'^'‘dx  (2-15) 

blT  -wQO 

f  S^(K^,K Je  X  y^dKxdKy  (2-16) 

•  »  -iK  x-iK  y 

S2(K^.Ky)  =  J_  J  JB2(x,y)e  ^dxdy  (2-17) 

41^^  .oa  -CO 

With  these  transform  relations  in  mind  it  is  seen  that  Equation  2-13  has 
the  form  of  Equation  2-17  and 

S  (K)  =  2Trr|A2AZ  {K^,K^=0)  (2-18) 

"  e 

Equations  2-11  and  2-18  are  general  relationships  between  PSD's  and  auto¬ 
correlation  functions  of  phase  and  electron  density  fluctuations. 

2.2.2  Power  Law  Phase  PSD 

A  number  of  different  phase  PSD's  are  available  for  use  in  the 
MPS  propagation  simulation.  The  most  useful  PSD  is  the  power  law  form 
with  phase  variance,  outer  scale,  and  inner  scale  specified  as  input.  The 
power  law  form  is  representative  of  many  in-situ  observations  (Dyson,  et 
al.,  1974;  Phelps  and  Sagalyn,  1976)  as  well  as  nunerical  simulations  of 
striation  spectra  (Scannapieco,  et  al.,  1976). 

For  this  case  the  two-dimensional  electron  density  PSD  has  the 
form  (Shkdrofsky,  1968) 
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where  is  the  Bessel  function  of  the  third  kind.  The  PSD  has  a 


(K2  + 


behavior  for  values  of  scale  size  between  the  inner  scale  and  the 
outer  scale  Lq.  The  corresponding  two-dimensional  autocorrelation 
function  is  given  as 

- (m-2)/2  - 

/l:Vx2  +  z2  +  ff?  j  K  (L:  Vx2  +  z2  +  je?  \ 

R  \  ° _ W  (ni-2)/2\o _ 

’  (£./L  ){'"-2)/2k  (£7L  ) 

T  0  (m.2)  2  1  0 


(2-20) 


The  corresponding  one-dimensional  phase  PSD  is  obtained  with  the  aid  of 
Equations  2-12  and  2-18  as 


where 


/  - V  m/  2 

K  (£,/L  )  (£.^2  +  1/L2 

(m.i)/2  •  0  \  1  X  of 


(2-21) 


,f  =  (2.«,Lo)>”  rhhi 

♦  "«  '  '  K(™-2)/2  («/L«) 


(2-22) 


For  the  often  used  case  of  a  K"^  one-dimensional  phase  PSD,  m  has  a  value 
of  three  and  in  the  limit  that  Lo>£i 


2(V)^  L^jAZ  o2 


(2-23) 


2.2.3  Gaussian  Phase  PSD 


A  second  useful  phase  PSD  is  the  Gaussian  form 


(2-24) 
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which  has  the  one-dimensional  autocorrelation  function 


B^(x)  =  a2  exp[.x2/L2]  (2-25) 

For  the  Gaussian  PSD  the  relationship  between  the  phase  variance  and  the 
electron  density  variance  is 

<p  '’g 

2.2.4  Phase-Screen  Generation 

In  this  subsection,  the  numerical  technique  to  generate  a  phase- 
screen  realization  for  the  numerical  MPS  propagation  code  is  discussed. 

The  goal  here  is  to  generate  a  stationary,  random  function  (|)(nAx)  which 
represents  the  phase  evaluated  along  the  MPS  grid.  The  spacing  along  the 
MPS  grid  is  ax  and  n  is  an  integer  representing  the  n'th  point  on  the 
grid.  Typically  2^'*  (16384)  discrete  points  are  used  to  represent  the 
functions  in  the  x-direction  transverse  to  the  direction  of  propagation. 

As  will  be  shown,  <ti{nAx)  is  generated  from  initial  knowledge  of  its  PSD. 

In  continuous  notation,  the  phase  may  be  written  as  the  Fourier  transform 

00  , 

<>(x)  =  /i(K)e’'^^dK  (2-27) 

-00 

In  the  discrete*  case,  as  used  in  the  MPS  code 
N-1 

(t(nAx)  =  I  i(mAK)e^^’'"’"^^AK  n=0,...,N-l  (2-28) 

m=0 


Now  if  the  Fourier  transform  of  the  phase  was  available,  a  phase-screen 
realization  could  be  easily  generated  by  using  Equation  2-28.  For  the 
moment  let  us  choose  as  the  Fourier  transform  of  the  phase  the  quantity 


*  In  this  section  all  discrete  sums  are  taken  over  a  range  of  the  index 
from  0  to  N-1. 
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(2-29) 


i(niAK)  =  rn,[S(mAK)L/2Tr] 

where  S(!iiaK)  represents  discrete  values  of  the  known  desired  phase  PSD. 

In  the  following  it  is  proven  that  this  choice  is  correct.  Here  aK  =  2tt/L, 
L  is  the  length  of  the  MPS  grid,  ax  =  L/N  and  therefore  axaK  =  2it/N.  r^, 
is  a  complex  number  given  as  the  sum  of  two  independent  Gaussian  random 
variables  with  zero  mean  and  variance;-,  of  unity 


■•n,  ■=  /1/2  (Sw  *  <ga)  (2-30) 

Successive  values  of  g^n,  and  gj^  may  be  obtained  numerically  by  samp¬ 
ling  from  a  pseudo-random  sequence  of  timbers  with  a  Gaussian  distri¬ 
bution.  The  factor  of  /Uz  is  included  so  that 


<r_r„>  = 
m  n 


mn 


(2-31) 


It  is  apparent  that  with  the  above  choice  f'or  r„,,  (|)(nAx)  is  the  sum  of  a 
sequence  of  Gaussian  variates  and  thus  it',  real  and  imaginary  parts  both 
have  a  Gaussian  or  normal  probability  dist - ibution.  Since  the  phase  of  an 
individual  phase-screen  is  a  real  quantity,  one  may  choose  either  the  real 
or  imaginary  part  of  (j.  computed  in  this  manner.  Or  since  the  real  and 
imaginary  parts  are  independent,  each  discrete  Fourier  transform  may  be 
used  to  generate  two  phase-screens  with  a  savings  in  computer  time. 


Now  since  the  r^,  are  independent,  it  is  apparent  that  the 
Fourier  coefficients  defined  in  Equation  2-29  are  also  independent.  Thus 
the  phase  power  spectriim  has  the  general  form 

<i{mAK)i*(nAK)>  «  «n,nS{mAK)  (2-32) 

which  is  a  requirement  if  the  phase  is  to  be  stationary  (Tatarskii, 

1971).  That  is,  if  Equation  2-32  is  satisfied,  then  the  phase  autocorre¬ 
lation  function  <(|i(mAx) (()*(nAx)>  depends  only  on  the  separation  distance 
(m-n)Ax  as  it  must  for  a  stationary  randor  variable. 

In  order  to  prove  the  validity  of  the  choice  of  Equation  2-29  to 
give  a  realization  of  the  phase-screen  it  is  convenient  to  calculate  the 


t. 
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discrete  autocorrelation  function.  Under  the  ergodic  hypothesis,  ensemble 
averages  or  expectations  of  random  fields  may  be  replaced  by  spatial 
averages.  Thus  the  phase  autocorrelation  function  may  be  written 


B.(5)  =  <(|.(x+5)<^*(x)> 

9 

1  L 

=  ^  f  (|i(x+c)(|>*(x)dx  (2-33) 

In  the  discrete  case  of  interest  here 

B  (kAx)  =-l  <t)(nAx+kAx) (ti*{nAx) AX  (2-34) 

^  L  n 

Now  the  Fourier  transform  of  (|)  and  ({>*  given  by  combining  Equations  2-28 
and  2-29  may  be  used  in  Equation  2-34  to  obtain 

B.(kAx)=-I  I  I(L/2iT)[S(mAK)S*(m'AK)]^/2 
L  n  m  m' 

X  r  r  ,  (2-35) 

m  m 


Now 

r  i2Trmn/N  -i2itm'n/N 
I  e  e 

n 

So  that  Equation  2-35  becomes 

B  (kAx)  -  I  S(mAK)e’^’^''/^  |r^|2  tK  (2-37) 

A  comparison  of  the  above  equation  to  the  continuous  relationship  between 
PSD  and  autocorrelation  function  as  given  by  Equation  2-14  shows  that  the 
PSD  of  the  numerically  generated  phase  is  S(mAK)  |rn,|  It  is  apparent 
that  different  values  of  the  index  m  correspond  to  different  spatial - 
frequency  components  of  the  PSD.  Since  |rn,l^  is  the  sun  of  the  squares 
of  two  Gaussian  variates,  each  of  the  Fourier  components  of  the  PSD  of  an 


m' 


(2-36) 
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individual  phase-screen  is  a  chi-squared  variate  with  two  degrees  of  free¬ 
dom  and  a  mean  value  of  S(m/J<).  Thus  for  any  given  phase-screen  realiza¬ 
tion,  the  phase  PSD  will  not,  in  general,  be  identical  to  the  desired 
PSD.  However,  the  average  PSD  of  many  such  phase  realizations  may  be 
obtained  by  taking  the  expected  value  of  Equation  2-37.  Since  <|rn,(^>  =  1 

AVG[B  (kAx)]  =  I  S{mAK)e’^’'"’'''^'^AK  (2-38) 

m 

Therefore  the  average  PSD  obtained  from  many  phase-screen  realizations  is 
the  desired  PSD,  S(mAK),  and  Equation  2-29  is  correct. 

2.2.5  Criteria  For  M’S  Application 

The  application  of  the  MPS  propagation  code  requires  that  the 
electric  field  and  the  phase  be  specified  at  a  discrete  number  of  grid 
points.  The  number  and  spacing  of  these  points  must  in  general  satisfy 
the  following  criteria:  1)  the  phase  distribution  of  a  phase-screen  must 
adequately  represent  the  actual  phase,  2)  the  wave  should  propagate  with¬ 
out  angular  aliasing,  and  3)  edge  effects  or  angular  scattering  off  the 
end  of  the  grid  should  be  minimal. 

2. 2. 5.1  Phase  Representation 


Adequate  phase  representation  is  assured  if  the  phase-screen 
length  L  is  at  least  5  to  10  times  as  large  as  the  phase  decorrelation 
distance  or  equivalently  the  outer  scale. 

L  >  5  L 

0 

The  spacing  between  grid  points,  ax,  should  be  several  times  smaller  than 
the  inner  scale.  Thus 

AX  <  Z./3 
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';h6  change  in  phase  from  one  grid  point  to  the  next  should  be  less  than  it 
to  satisfy  the  Nyquist  sampling  theorem.  Mathematically 

<{i(^2^  ”  ^11  »  X2  “  Xj  —  AX 

In  terms  of  the  distance  between  sample  points  ax,  this  constraint  can  be 
v-ritten 


This  inequality  can  be  written  in  an  rms  sense  as  follows.  Under  the 
ergodic  hypothesis  discussed  previ''usly  the  autocorrelation  function  can 
be  written 


B  (C)  <(t!(x+c)4>*{x)>  =-  J  (|>(x+0<t>*{x)dx 

L  0 


(2-39) 


Taxing  the  second  derivative  with  respect  to  c  and  integrating  by  parts 


/d'»(x+0 

\  dx  dx  / 


(2-40) 


Ecuacion  2-40  states  the  relationship  between  the  autocorrelation  function 
the  derivative  of  the  phase  and  the  second  derivative  of  the  phase 
autocorrelation  function.  In  a  mean  sqrare  sense  then 


d^B  (c) 
<!> 


d^2 


C=0 


Thus  in  an  rms  sense  the  above  li:„it  on  the  grid  spacing  ax  may  be 
expressed  as 


AX  <  TI  (  "  - ^ - 


\ 


riC‘ 


1/2 


C=0 
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For  a  Gaussian  phase  PSD  the  second  derivative  term  may  be  calculated  from 

Equation  2-25  as  so  that  the  necessary  condition  is 

0 

AX  <  ttL  //2  0, 

0  (p 

2. 2. 5. 2  Wave  Propagation 

In  order  to  adequately  represent  propagation  in  free-space  by 
the  use  of  the  Fourier  transform  relationship  of  Equation  2-7,  it  is 
necessary  that  the  various  functions  involved  be  accurately  sampled.  To 
satisfy  the  Nyquist  sampling  criterion  the  difference  in  the  function 

k2az 

"TF 

must  be  less  than  tt  when  evaluated  from  one  value  of  K  to  the  next.  Since 
the  maximun  value  of  K  is  wN/L  where  N  is  the  nunber  of  grid  points  and  L 
is  the  grid  length,  the  necessary  condition  is 

-  -  - ! - - -  <  TT 

2kL2  2kL2 
or 

tt^NaZ  . 

-  V  TI 

kL2 

or 

AZ  <  2Lax/x 

This  condition  may  be  relaxed  somewhat  in  practice  if  the  phase  PSD 
is  very  small  at  large  values  of  the  wavenumber,  K.  In  that  case  the  high 
wavenumber  values  are  cut  off  by  the  phase  power  spectrum  and  small 
inaccuracies  in  the  phase  are  relatively  unimportant. 

2 . 2 . 5 . 3  Edge  Effects 

Because  of  the  discrete  nature  of  the  MPS  grid  representation, 
energy  leaving  one  side  of  the  grid  appears  on  the  other  side.  As  the 
propagation  distance  z  increases  this  effect  becomes  more  important.  This 
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effect  must  be  controlled  for  a  propagation  solution  to  be  valid, 
the  scattering  angle  is  given  by 


Since 


e  = 

k  dx 


(2-41) 


where  ^  is  the  phase  at  a  phase-screen,  the  energy  scattered  at  an  angle  e 
travels  a  distance  ze  perpendicular  to  the  direction  of  propagation  after 
propagating  a  distance  z.  To  adequately  insure  against  edge  effects  it  is 
necessary  that  the  MPS  grid  size  L  be  greater  than  ze  for  each  propagation 
step.  Thus 

L  >i 

k  dx 

This  expression  may  be  represented  in  a  mean  square  sense  as  in  Section 
2. 2. 5.1  as  , 

.  7  I  d^B^(g)  I  \i/2 


L  >  i 


For  a  Gaussian  phase  PSD,  the  second  derivative  is  so  that  the 

necessary  condition  is 


/2  ttL. 


to  minimize  edge  effects. 

2.3  APPLICATION  TO  A  GAUSSIAN  LENS 

In  this  section  the  numerical  M!’S  propagation  code  is  applied  to 
the  problem  of  scattering  by  a  single,  deterministic  Gaussian  phase  lens. 
For  this  example  of  scattering  from  a  non-random  structure,  theoretical 
results  are  obtained  using  the  Fresnel -Kirchhoff  integral  equation  and  are 
shown  to  be  identical  to  the  numerical  MPS  results. 

For  this  example  a  single  Gaussian  phase  lens  given  by 


<|,(x)  =  -  <t,Q  exp(-x2/r2) 


(2-42) 
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is  used  at  the  first  phase-screen  location  (z=0)  with  (t)Q  chosen  as  10 
radians  and  set  equal  to  the  wavelength  x.  Since  (t)(x)  is  negative 
this  Gaussian  lens  corresponds  to  a  decrease  in  the  ionization  level  below 
the  ambient  and  acts  as  a  convergent  or  focusing  lens.  A  positive  phase 
would  correspond  to  an  increase  in  ionization  and  K)uld  act  as  a  divergent 
lens. 

The  intensity  I  =  jEl^  was  observed  as  a  function  of  distance 
from  the  lens,  z/x  as  indicated  in  Figure  2-3.  In  this  figure  the 
observation  screens  are  aligned,  so  that  features  of  the  diffraction 
pattern  can  be  easily  followed  for  changing  propagation  distance  z/x. 

The  focal  length  for  a  Gaussian  lens  is  given  by  (Salpeter, 

1967) 

F  =  kr2/2<t,^  (2-43) 

where  k  is  the  wavenumber  Zn/x.  Taking  r  =  x  and  *  =10  radians,  the 
focal  length  F/x  is  0.31  which  corresponds  quite  well  with  the  value  of 
z/x  =  0.5  ,  shown  in  the  figure  where  the  intensity  in  the  diffraction 
pattern  at  x  =  0  builds  up  to  a  maximum.  At  values  of  z/x  greater  than 
the  focal  length,  the  diffraction  pattern  exhibits  increasingly  more  com¬ 
plex  patterns  associated  with  interfering  rays  coming  from  the  edges  of 
the  lens,  rather  than  from  the  center. 

An  analytic  relationship  between  the  electric  field  at  z  =  Zj 
and  z  =  Z2  is  given  by  the  Fresnel -Kirchhoff  integral  (Ratcliffe,  1956) 

E(x,Z2)  =  (-i2,,(Z2-z,)/k]-‘'2 

X  exp(-ik(x-e)2/2(z2-Zi)  I  E(5,Zi)dc  (2-44) 

where  E(x,Zi)  is  the  electric  field  as  a  function  cf  x  in  the  z  =  z^ 
plane.  For  an  initial  electric  field  in  the  z  =  0  plane  given  by 

E(x,0)  =  exp[-i<j>^  exp(-x2/r2)  ] 
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Equation  2-44  may  be  written  as  an  infinite  series  by  expanding 
exp[-i(t)oexp(-x2/rQ2)  ]  in  a  Taylor  series  as 

00  n 

e*  =  J  *  (2-45) 

n=0 

The  resulting  integral  over  ^  may  then  be  analytically  performed 
with  the  use  of  standard  integral  tables  (Gradshteyn  and  Ryzhik,  1965). 

The  resulting  electric  field  in  the  z  =  22  '"sy  be  expressed  as 

E(x,z)  =  exp[ikz+ikxV2z] 

/ikz 


This  series  is  easily  summed  numerically  and  the  results  for  several 
values  of  z/A  are  shown  in  Figure  2-3  to  be  identical  to  the  MRS  code 
calcul ation  for  z/x  _<  2. 

For  z/x  =  10  the  multiple  phase-screen  calculations  deviate  from 
the  theoretical  results  because  energy  which  has  left  one  side  of  the  grid 
is  coming  back  into  the  other  side  (this  is  the  well  known  wrap-around  or 
aliasing  effect  in  discrete  Fourier  transforms).  This  point  illustrates 
that  care  must  be  taken  in  the  application  of  the  MRS  techniques  to  avoid 
erroneous  numerical  results. 

2.4  APPLICATION  TO  A  RANDOM  PHASE-SCREEN 

This  section  contains  a  second  application  of  the  multiple 
phase-screen  propagation  code  to  a  problem  for  which  analytic  results  are 
available  -  propagation  through  a  thin  phase-screen  characterized  by  a 
Gaussian  power  spectral  density.  For  this  example,  only  one  phase-screen 
is  used  to  characterize  the  random  medium.  The  screen  is  located  at  z  =  0 
and  is  generated  with  a  phase  power  spectral  density  (PSD)  of 
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Figure  2-3(a}.  Diffraction  pattern  of  a  Gaussian  lens— nunerical 
results  from  simulation  on  left-hand  side  versus 
theoretical  results  on  right-hand  side. 
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(2-47) 
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where  the  scale  size  Lq  is  chosen  as  X  and  the  phase  standard  deviation 
0^  takes  on  the  values  0.1,  1.0,  and  10.0  radians.  The  MPS  grid  here  is 
comprised  of  2048  cells  which  represent  a  total  length  of  30x. 

In  the  MPS  simulation  a  plane  wave  then  propagates  through  the 
single  phase-screen  and  is  observed  at  a  number  of  different  receiver 
plane  locations  corresponding  to  different  values  of  the  propagation  dis¬ 
tance  from  the  phase-screen.  At  each  receiver  plane  location  defined  by  a 
value  of  z^,  the  received  electric  field  E(x,Zp)  is  observed. 

To  obtain  statistical  results  for  the  received  signal,  the  simu¬ 
lation  is  exercised  ten  times  with  ten  different  phase-screen  realizations 
(each  screen  is  generated  using  a  different  sequence  of  random  nimbers  as 
discussed  previously).  Figure  2-4  shows  a  plot  of  one  of  the  ten  phase- 
screen  realizations  as  a  function  of  x/Lq,  the  normalized  distance  along 
the  MPS  grid  perpendicular  to  the  direction  of  propagation.  No+e  the 


Figure  2-4.  Realization  of  phase  for  a  Gaussian  PSD. 
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rather  smooth  plot  with  very  little  rapid  phase  variation.  This  plot  is 
characteristic  of  a  Gaussian  PSD.  The  figure  shows  the  phase  over  the 
entire  MPS  grid;  note  that  the  phase  is  generated  to  be  continuous  at  the 
edges  of  the  grid.  This  behavior  is  necessary  to  avoid  errors  associated 
with  taking  the  Fourier  transforms  of  discontinuous  functions.  Phase- 
screen  realizations  generated  by  the  MPS  code  are  always  continuous  as 
demonstrated  here.  Similarly,  the  received  electric  field  also  possesses 
this  same  type  of  continuity  across  the  boundaries  of  the  MPS  grid. 

The  code  has  input  options  to  provide  for  some  statistical 
analysis  cf  the  generated  phase-screens.  Figures  2-5  and  2-6  show  a 
comparison  of  the  desired  phase-screen  PSD  and  autocorrelation  functions 
versus  the  computed  PSD  and  autocorrelation  functions  averaged  over  ten 
realizations.  The  phase  standard  deviation,  a^,  was  1  radian  in  the 
plots.  As  can  be  seen  in  the  plots,  the  agreement  between  the  desired  PSD 
and  autocorrelation  functions  and  their  numerically  generated  equivalents 
is  quite  good. 

Figure  2-7  shows  another  example  of  a  phase-screen  realization 
which  is  generated  using  the  same  string  of  pseudo-random  numbers  as  in 
the  previous  case.  Here  the  phase  PSD  is  a  K"^  power  law  form  with  outer 
scale  equal  to  the  wavelength.  Note  that  the  large  scale  structure  in 
both  screens  is  quite  similar  since  both  are  based  on  the  same  random 
number  sequence.  The  appearance  of  the  small  scale  structure  in  Figure 
2-7  is  caused  by  the  relatively  slow  fall-off  of  the  power  law  PSD  with 
increasing  wavenumber  (decreasing  scale  size).  Figure  2-8  shows  a  com¬ 
parison  of  the  MPS  generated  phase  PSD  averaged  over  ten  realizations  with 
the  desired  power  law  PSD.  Again,  in  keeping  with  the  results  shown  in 
the  previous  example,  the  agreement  between  nunerical  and  desired  results 
is  quite  good. 

An  additional  check  on  the  results  of  the  MPS  code  for  a  Gaus¬ 
sian  PSD  is  made  possible  by  comparing  analytic  results  for  the  scintil¬ 
lation  index  with  theoretical  results.  The  S^  scintillation  index  is 
defined  as  the  normalized  variance  of  the  received  power  lEp,  where  E  is 
the  received  complex  voltage 

S.'-  <(|eP-<|E|2>)2)  /<|E|2>2  (2-48) 
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CORRELATION  FUNCTION 


Figure  2-5.  Comparison  of  numerically  generated  phase  PSD  with 
desired  Gaussian  PSD. 


Figure  2-6.  Comparison  of  correlation  function  of  numerically  generated 
phase  with  desired  analytic  correlation  function. 


Figure  2-9,  scintillation  index  for  Gaussian  spectrum.  The  solid  lines 
are  analytic  approximations  valid  for  the  weak-scatter  case. 


For  weak  scattering  is  given  in  the  thin  phase-screen  approximation 
(Salpeter,  1957)  as 


00 

=4/S  (K)  sin2[K2z/2k]dK  (2-49) 

-00  ^ 


This  integral  is  easily  evaluated  for  the  case  of  a  Gaussian  PSD  and  is 
plotted  as  a  function  of  z/x  in  Figure  2-9  for  the  three  values  of  of 
0.1,  1.0,  and  10.0  radians.  The  nunerical  results  are  shown  as  circles  in 
this  figure  and  connected  with  dashed  lines.  Note  the  excellent  agreement 
for  small  values  of  and  small  values  of  the  propagation  distance  z/x 
where  the  weak  scattering  theory  holds.  For  the  larger  values  of  z/x,  and 
>  1.0  radian  the  weak  scattering  result  of  Salpeter  is  no  longer  valid; 
in  this  case  the  scintillation  index  saturates  at  a  value  of  unity  as  is 
predicted  on  theoretical  grounds  (Fante,  1975)  and  as  is  predicted  by  this 
and  other  numerical  MPS  calculations. 
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2.5 


WIDE  BANDWIDTH  NPS  CALCULATIONS 


In  the  preceding  formulation  and  examples,  the  MPS  propagation 
simulation  was  applied  to  a  single  wave  frequency.  For  dispersive  or  fre¬ 
quency  selective  ionospheric  communication  channels,  it  is  necessary  to 
first  determine  the  Fourier  components  or  spectrun  of  the  transmitted 
signal.  Then  each  of  these  Fourier  components  is  propagated  through  an 
identical  sequence  of  phase-screens.  If  (t)£(x)  is  the  phase  distribution 
for  the  center  or  carrier  frequency,  toj.,  then  the  corresponding  phase 
shift  at  another  frequency  in  the  signal  bandwidth  is 

*  (x)  =  *  (x)  (2-50) 

“  0) 

Once  each  signal  Fourier  component  has  been  transmitted,  then  a  Fourier 
transform  is  applied  to  obtain  the  resulting  time-domain  signal  at  the 
receiver. 

Now  if  the  transmitted  waveform  consists  of  a  complex  envelope 
modulated  on  a  carrier  it  is  shown  in  Section  3  that  the  received  waveform 
can  be  expressed  in  terms  of  a  complex  modulation  envelope  as 

v(x,2^,t)  =  Re  je(x,Zj.,t-tj)e  °  °  (2-51) 

where  the  complex  envelope  is  given  by  the  Fourier  transform 

CO 

e(x,z^,T)  =  (1/2ti)  J" M(v)  U(x,z^,v+a)p)e^'^^dv  (2-52) 

«ao 

M(v)  is  the  spectrum  of  the  transmitted  signal  envelope  and  U(x,Zp,v)  is 
the  solution  to  the  parabolic  wave  equation  for  the  case  of  irregularities 
infinitely  elongated  in  the  y-direction.  For  simplification  in  this 
section  it  is  assumed  that  mean  ionization  effects  are  negligible  in  com¬ 
parison  to  stochastic  effects.  In  that  case  U(x,z^,v)  is  a  realization 
of  the  solution  to  the  parabolic  wave  equation  given  by  the  MPS  simulation 
for  a  particular  value  of  frequency  u.  0q  is  the  total  mean  phase  shift 
(at  the  carrier  frequency)  and  t^j  is  the  mean  group  time  delay.  These 
quantities  are  discussed  in  detail  in  Section  3  of  this  report. 
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In  the  following  sections,  the  MPS  propagation  simulation  is 
applied  to  obtain  time-domain  results  for  wide  bandwidth  signals  after 
propagation  through  ionization  irregularities.  Two  interesting  cases  are 
considered.  In  the  first  case  time-domain  results  for  scattering  by  a 
strong  Gaussian  lens  are  presented  and  compared  to  theoretical  results  in 
order  to  verify  the  MPS  time-domain  calculations.  In  the  process  some 
interesting  properties  of  the  scattering  process  are  noted  and  discussed. 
In  the  second  case  the  MPS  code  is  used  to  predict  the  scattering  proper¬ 
ties  of  an  ionized  barium  cloud  through  which  a  GPS-like  pseudo-noise 
spread-spectrim  navigation  signal  is  transmitted  on  a  carrier  at  100  MHz. 

2.5.1  Time-Domain  Propagation  Through  a  Strong  Gaussian  Lens 

As  an  example  of  the  MPS  propagation  simulation  when  frequency 
selective  effects  are  important,  the  time-domain  diffraction  pattern  of  a 
Gaussian  lens  is  numerically  determined  and  compared  to  theory.  In  this 
example  the  single  Gaussian  phase  lens  is  given  by 

f 

♦(x,f)  .  -£^exp  (.(x-x^)Vr2)  (2-53) 


where  f^.  is  the  carrier  or  center  frequency  and  is  100  MHz  here,  (j)Q  is 
70  radians,  the  lens  is  located  at  the  center  of  the  50  km  MPS  grid  and 
has  a  half-width  rg  of  400  m.  Since  the  phase  here  is  positive  this 
Gaussian  lens  corresponds  to  an  increase  in  ionization  above  the  ambient 
and  thus  is  a  defocusing  or  divergent  lens. 


The  transmitted  waveform  consists  of  a  single,  band-limited  tri¬ 
angular  pulse  modulated  on  a  carrier  at  100  MHz.  The  spectrum  of  the 
transmitted  pulse  is  given  by 


where  T^.  is  roughly  twice  the  total  duration  of  this  pulse. 
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If  the  plilse  were  not  band-limited  to  ±  l/T^.,  the  time  domain 
waveform  would  be  giv'en  as 


NsirKojT  /2\2 

- :4-  1  e  (2-55) 

»T  /2\,/ 

•  00  ^ 

\, 


(2-56) 


where  oj  =  2Trf  and  the  subscript  <»  refers  to  the  infinite  bandwidth  case. 
In  the  case  under  consideration  where  the  spectrum  is  limited  to  a  total 
bandwidth  oi  2/T(.,  corresponding  tc  the  first  nulls  of  the  spectrum 
M(f),  the  resulting  transmitted  time-domain  waveform  lacks  the  sharp  tri 
angle  edgos  at  t  =  0  and  ±  and  is  a  smooth  function  with  continuous 
first  derivatives  at  these  points.  Also  a  loss  of  approximately  0.89  dB 
occurs  at  T  =  0  for  the  band-limited  case  as  compared  to  the  infinite 
bandwidth  case. 


The  case  of  a  single  transmitted  triangular  waveform  described 
above  has  an  exact  analog  to  the  case  of  a  rang.ig  system  that  utilizes  a 
pseudo-noise  spread  spectrum  ranging  signal.  In  this  type  of  ranging 
system  a  high  bandwidth  phase-shift  keyed  signal  is  transmitted  as  a 
random  sequence  of  pulses  called  a  pseudo-noise  (PN)  code.  A  PN  spread 
spectrum  receiver  is  then  used  to  correlate  a  local,  receiver  generated 
code  with  the  actual  received  code.  This  correlation  operation  is  per¬ 
formed  in  the  receiver  code  correlator  hardware.  The  code  correlator 
output  can  be  shown  (Knepp  and  Bogusch,  1979;  Bogusch,  et  al.,  1981)  to  be 
identical  to  Equation  2-51  where  Equations  2-52  and  2-54  define  the  appro¬ 
priate  quantit.es  required  in  Equation  2-51. 


Therefore  the  time-domain  results  to  be  presented  in  this 
section  can  be  interpreted  in  either  of  two  equally  correct  manners.  On 
the  one  hand  the  results  for  the  received  signal  may  be  interpreted  as  the 
received  envelope  after  propagation  of  a  single,  band-limited  modulated 
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triangular  pulse  through  a  disturbed  channel.  On  the  other  hand,  these 
same  results  may  be  interpreted  as  the  output  of  the  code  correlator  of  a 
spread  spectrum  receiver.  In  the  latter  interpretation  T(.  is  the  dura¬ 
tion  of  a  chip  of  the  pseudo-noise  code  and  l/T^.  is  thus  the  chip  rate. 

In  this  example  the  chip  rate  is  taken  as  10.23x10®  bps  which  is  the  rate 
used  in  the  Global  Positioning  System  (GPS).  Thus  the  value  of  is 
approximately  97.75  nsec  and  the  two-sided  bandwidth  2/T{-  which  effec¬ 
tively  limits  the  spectral  region  of  interest  is  20.46  MHz  for  this 
example. 


7or  the  case  of  no  scintillation,  U(x,z^, v+u^^)  is  unity  and 
the  received  signal  is  a  smooth  band-limited  triangle.  That  is,  for  no 
scintillation 

U(x,z^,v+a)Q)  =  1  (2-57) 

and  from  Equation  2-52 

e(x,2  ,t)  =  L  /M(v)e*“'dv  .  (1-|t|/\)bl  (2-58) 

tTT  ..00 

where  the  subscript  BL  indicates  that  the  triangular  waveform  is  band- 
limited. 


For  a  non-selective  fading  case,  where  the  MPS  solution  to  the 
parabolic  wave  equation,  U(x,z^,v+^3)  is  independent  of  frequency  v 
but  the  amplitude  and  phase  of  U  are  functions  of  x,  U  can  be  written  as 

U(x,z^,v+a)Q)  =  U(x,z^)  (2-59) 

and  the  received  complex  envelope  is  then 

e(x,z^,T)  =  U(x,z^)(l-|T|/T^)gL  (2-60) 

In  this  case  the  complex  envelope  retains  the  triangular  form  of  the 
transmitted  waveform  in  delay,  but  varies  in  amplitude  and  phase  as  a 
function  of  x,  the  distance  along  the  MPS  grid. 
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Note  that  the  received  signal  here  is  a  function  of  x,  the  dis¬ 
tance  along  the  MPS  grid,  z^,  the  propagation  distance,  and  x,  the  time 
delay.  For  a  moving  receiver,  or  for  motion  of  the  striations,  the  func¬ 
tional  dependence  on  x  is  converted  into  time  dependence  by  inclusion  of 
an  effective  velocity  which  converts  motion  of  the  line-of-sight  path  into 
motion  along  the  MPS  grid.  At  a  single  position  x  or  a  single  instant  in 
time  x/Vgff,  it  is  possible  to  receive  a  signal  that  is  non-zero  over  a 
range  of  delay  values,  x.  This  situation  is  the  case  when  angular  scatter¬ 
ing  acts  to  cause  multipath  effects  where  the  signal  propagates  over  a 
number  of  slightly  different  paths  and  thus  arrives  at  the  receiver  at 
different  delays  corresponding  to  the  different  paths.  The  difference  be¬ 
tween  time  and  time-delay  in  this  context  is  discussed  further  in  Section 
3. 2. 3.1. 


In  the  selective  fading  case  U(x,Zj.,\H-a^)  is  a  random  func¬ 
tion  of  both  X  and  v  and  the  triangular  waveform  suffers  distortion  since 
each  spectral  component  of  the  signal  experiences  different  degrees  of 
fading.  The  top  portion  of  Figure  2-10  shows  an  example  of  the  amplitude 
of  the  received  triangular  pulse  for  the  case  of  non-selective  fading. 

The  bottom  portions  of  the  figure  show  examples  of  the  received  signal 
during  selective  fading  conditions.  The  resulting  distorted  waveforms  may 
be  interpreted  as  being  the  superposition  of  a  number  of  triangular 
signals  which  have  propagated  over  different  geometric  paths  to  reach  the 
receiver  at  a  number  of  different  time  delays.  Thus  the  distortion  shown 
is  always  extended  or  elongated  in  the  direction  of  increasing  or  positive 
time  delay. 

Most  of  the  results  to  be  presented  for  this  deterministic  exam¬ 
ple  are  obtained  using  the  Fresnel -Kirchhoff  Equation  2-46  valid  for  a 
single  Gaussian  lens  discussed  previously.  These  results  were  compared  to 
MPS  calculations  and  found  to  agree.  With  this  in  mind  the  two  different 
but  equivalent  calculation  methods  will  no  longer  be  distinguished. 

Now  in  the  example  under  consideration,  a  single  phase-screen 
defined  by  16384  grid  points  over  a  distance  of  50  km  is  used  to  represent 
the  large  Gaussian  lens  described  by  Equation  2-53.  A  plane  wave 

E(x,z=0)  =  1  (2-61) 


53 


AMPLITUDE 


TIME  DELAY  (CHIPS) 


Figure  2-10.  Example  of  envelope  of  wideband  signal  after 
propagation  through  nonselective  and  frequency 
selective  propagation  channels. 


is  incident  on  the  lens.  64  distinct  frequency  components  are  used  to 
represent  the  signal  spectrum  over  the  null-to-null  bandwidth  of  20.46  MHz 
centered  at  the  carrier  frequency  of  100  MHz.  Each  of  these  frequency 
components  of  the  signal  is  then  propagated  through  the  environment  which 
consists  of  the  single,  strong  Gaussian  lens.  The  solution  U{x,Zj.,v+(i:^) 
to  the  parabolic  wave  equation  is  obtained  via  MPS  simulation  at  64  values 
of  the  frequency  v  and  several  values  of  the  propagation  distance 
To  obtain  time-domain  results  for  the  received  signal  envelope  Equation 
2-52  is  applied  directly. 

First  consider  the  propagation  of  the  single  discrete  frequency 
component  at  100  MHz.  Figures  2-11  to  2-15  show  the  received  signal 
amplitude  and  phase  for  the  100  MHz  frequency  component  of  the  spectrin  at 
distances  of  3,  5,  10,  50,  and  100  km  from  the  lens  location.  It  can  be 
seen  that  the  effect  of  this  defocus ing  lens  is  to  cause  energy  to  scatter 
away  from  the  center  of  the  lens  toward  the  edges.  Thus  at  distances  of  3 
and  5  km  the  energy  can  be  seen  to  focus  near  the  edges  of  the  lens  with  a 
reduced  signal  level  immediately  below  the  lens  center  at  a  distance  of 
zero  on  the  figure.  The  geometric  focal  distance  of  a  convergent  lens  of 
this  size  is  2.4  km  (Equation  2-43)  which  corresponds  reasonably  well  with 
the  edge  focusing  observed  at  5  km.  At  a  propagation  distance  of  100  km 
it  can  be  seen  from  Figure  2-15  that  the  signal  level  has  faded  16  dB  over 
a  large  2  km  region  just  below  the  center  of  the  lens.  At  100  km  the  usual 
ringing  which  occurs  in  any  diffraction  pattern  is  evident.  Also  note  from 
Figure  Z-’S  that  Ine  signal  level  returns  to  0  dB  for  distances  greater 
than  10  km  from  the  lens  center.  This  corresponds  to  the  incident  plane 
wave  which  has  propagated  past  the  lens  and  remains  undisturbed  by  tne 
lens. 


Figures  2-16  to  2-20  show  the  magnitude  of  the  envelope  of 
the  received  time-domain  signal  |e(x,Zr,T)l  for  propagation  distances  of 
3,  10,  50,  and  100  km  from  the  Gaussian  lens.  In  these  three-dimensional 
figures  the  time-domain  signal  received  at  each  succeeding  value  of  x  is 
plotted  behind  the  preceding  received  time-domain  signal.  Thus  the  MPS 
x-axis  is  now  directed  into  the  figure. 

In  Figure  2-16  a  4  km  interval  centered  immediately  below  the 
lens  center  is  shown.  Near  the  two  ends  of  the  x  interval  shown  the 
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Figure  2-13.  Signal  amplitude  and  phase  at  100  MHz,  10  km  from 
a  strong  Gaussian  lens. 
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received  signal  is  identical  to  the  transmitted  signal  and  has  the  smooth 
shape  of  a  band-limited  triangular  waveform.  Near  the  central  portion  of 
the  figure  the  received  signal  is  delayed  relative  to  the  unperturbed 
signal.  This  effect  is  caused  by  the  effectively  greater  ionization 
sampled  by  rays  which  travel  through  the  center  of  the  lens  and  hence 
experience  increased  time  delay.  Figure  2-16  also  exhibits  the  effects  of 
focusing  which  are  shown  as  increases  in  signal  level  near  the  edges  of 
the  lens.  This  behavior  was  discussed  previously  in  conjunction  with 
Figure  2-11. 

As  the  propagation  distance  increases  the  effects  of  diffraction 
become  more  evident.  Figure  2-19  shows  the  received  signal  at  a  propaga¬ 
tion  distance  of  100  km.  The  two  main  effects  of  the  Gaussian  lens  are 
evident  in  this  figure.  First,  in  the  central  region  of  the  figure  cor¬ 
responding  to  the  region  immediately  below  the  lens  center,  there  is  a 
minimim  in  received  power  caused  by  the  outward  scattering  of  waves  away 
from  the  lens  center.  Secondly,  two  secondary  waves  delayed  up  to  12 
chips  (1  chip  =  97.75  nsec)  with  respect  to  the  unperturbed  wave  are 
apparent.  These  two  secondary  waves  which  appear  on  the  figure  as  large 
intersecting  semicircles  correspona  to  two  outward  propagating  spherical 
waves  which  originate  from  points  on  the  two  edges  of  the  original 
Gaussian  lens,  i.e.,  at  x  =  ±  400  m. 

For  the  MRS  calculation  of  the  scattering  properties  of  this 
Gaussian  lens,  16384  grid  points  are  used  to  represent  a  spatial  region 
50  km  in  extent  in  the  x-direction.  In  Figure  2-19  only  101  of  these  grid 
points  are  shown  over  a  region  of  20  km.  Thus  much  of  the  fine  detail 
available  in  the  calculation  is  not  shown  in  Figure  2-19. 

Figure  2-20  shows  a  close-up  of  the  MRS  results  for  a  propaga¬ 
tion  distance  of  100  km.  The  3  km  region  which  extends  to  the  center  of 
the  MRS  grid  immediately  below  the  center  of  the  Gaussian  lens  is  shown. 
The  fine  details  of  the  diffraction  pattern  at  zero  delay  are  shown  as 
well  as  the  reduced  signal  level  directly  below  the  center  of  the  lens. 
Also  note  that  the  time  delay  of  the  wave  that  propagates  directly  through 
the  center  of  the  lens  is  correctly  given  by  (j)Q/2Trf  which  for  a  phase  of 
70  radians  and  a  frequency  of  100  MHz  is  1.11x10"^  sec  or  1.14  chips  as 
shown  in  the  figure. 
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Figure  2-16.  Envelope  of  time-domain  waveform  3  km  from  a 
strong  Gaussian  lens. 
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Figure  2-17.  Envelope  of  time-domain  waveform  10  km  from  a 
strong  Gaussian  lens. 
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Figure  2>18.  Envelope  of  time-domain  waveform  50  km  from  a 
strong  Gaussian  lens. 
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Figure  2-19.  Envelope  of  time-domain  waveform  100  km  from  a 
strong  Gaussian  lens. 
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This  comprehensive  example  involves  direct  comparisons  of  the 
time-domain  results  from  the  MPS  propagation  simulatioi'.  with  theoretical 
results  for  the  case  of  a  strong  dispersive  Gaussian  lens.  Although  the 
results  are  presented  here  as  a  straightforward  example,  many  direct  com¬ 
parisons  were  performed  during  the  work.  Agreement  between  numerical  MPS 
and  theoretical  results  was  obtained  in  all  cases. 

2.5.2  Time-Domain  Propagation  Through  a  Striated  Barlun  Cloud 

In  this  section  the  MPS  propagation  simulation  is  applied  to 
predict  the  effect  of  a  large,  striated  bariun  cloud  on  a  wide  bandwidth 
spread  spectrum  signal.  Here  the  propagation  environment  caused  by  a 
striated,  hour  old,  bariun  cloud  is  modeled  as  a  single  phase-screen.  The 
phase-screen  model  is  based  on  analysis  oi’  the  data  taken  during  the 
STRESS  bariun  release  series  (Prettie  and  Marshall,  1978)  48  minu*‘es  after 
the  release  of  barium  cloud  Esther.  The  propagation  results  are  intended 
to  provide  predictions  of  the  effects  to  be  experienced  during  the  PLACES* 
barium  release  experiment  scheduled  for  winter  1980.  Back  propagation 
data  provided  by  Dr.  Prettie  were  utilized  to  provide  a  phase-screen 
model*  of  Esther  at  48  minutes  that  is  comprised  of  two  parts:  1)  a  large 
deterministic  Gaussian  cloud  which  represents  the  unstriated  portion  of 
the  overall  barium  cloud  and,  2)  the  striations  represented  by  a  K"^ 
power-law  phase  PSD  with  phase  stanuard  deviation  of  33  radians  at 
the  100  MHz  carrier  frequency,  with  an  outer  scale  of  390  m  and  an  inner 
scale  of  10  m.  The  deterministic  Gaussian  lens  is  given  by 

280fr 

(p(x,f)  =  — ^  exp{-[(x-25)/5]^}  radians  (2-62) 

fj.  is  the  carrier  frequency  of  100  MHz  and  the  distance  x  is  given  in 
kilometers.  Figure  2-21  shows  the  phase-screen  realization  of  the  bariun 
cloud  used  in  the  MPS  simulation.  The  entire  50  km  length  of  the  MPS  grid 
is  shown  in  this  figure.  The  component  of  the  phase  due  to  the 


*  The  details  of  this  model  are  fully  discussed  in  "Predictions  of  GPS 
X-set  Performance  During  the  PLACES  Experiment,"  by  Knepp  and  Bogusch, 
1979,  PLACES  is  the  acronyn  for  Position  Location  and  Communication 
Effects  Simulation. 
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deterministic  Gaussian  lens  is  plotted  separately  in  the  figure.  Note 
that  the  phase-screen  shown  represents  only  a  single  realization  of  the 
bariun  cloud  and  thus  leads  to  only  one  realization  of  the  received 
signal.  To  obtain  statistically  significant  results,  a  nunber  of  differ¬ 
ent  phase-screen  realizations  would  have  to  be  generated  using  different 
random  number  sequences  and  the  results  of  the  propagation  calculations 
would  have  to  be  averaged. 

In  the  PLACES  experiment  the  transmitted  signal  is  a  spread- 
spectrum  pseudo-noise  code  at  a  chip  rate  of  10.23  Pbps.  The  chip  dura¬ 
tion  T(.  is  again  97.75  nsec  and  the  receiver  is  band-limited  to  the 
null-to-null  bandw'dth  of  the  transmitted  code,  20.46  MHz.  Thus  the 
three-dimensional  plots  of  the  received  signal  to  be  shown  in  this  example 
are  interpreted  as  the  ouv.put  of  the  receiver  code  correlator  although  no 
error  is  committed  using  the  interpretation  of  transmitted  and  received 
triangular  pulses  discussed  previously.  As  before  16384  MPS  grid  points 
are  used  to  represent  an  actual  grid  50  km  in  length.  Again  64  distinct 
frequency  components  are  used  to  represent  the  20.46  MHz  bandwidth. 
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For  this  example  results  are  obtained  with  the  deterministic 
Gaussian  lens  alone  as  well  as  with  the  entire  bariun  cloud  model.  Figure 
2-22  shows  the  time-domain  signal  scattered  by  the  deterministic  Gaussian 
lens  alone  at  a  propagation  distance  of  100  km.  As  the  central  point 
directly  below  the  center  of  the  barium  cloud  is  approached  from  either 
side,  increasing  delay  is  experienced  because  of  the  increasing  ionization 
(in  the  form  of  increased  phase)  along  the  propagation  path.  Directly 
below  the  cloud  center  the  delay  should  be 

_i_  = - 280 - ^  ^^5  g  (2-63) 

2irfc  2it  •  100  X  10^ 

or  4.56  chips  as  shown  on  the  plot.  Also  note  the  increased  pulse  spread¬ 
ing  in  the  time-domain  which  is  apparent  for  that  part  of  the  received 
signal  subject  to  the  largest  delay.  This  dispersive  spreading  is  caused 
by  the  increased  ionization  near  the  center  of  the  Gaussian  lens.  Finally 
this  large  barium  cloud  is  also  acting  as  a  divergent  lens  to  direct  rays 
away  from  the  lens  center  towards  the  edges  of  the  lens.  This  slight  edge 
focusing  is  apparent  in  the  increased  signal  level  near  the  barium  cloud 
boundaries  about  5  km  from  the  center  point. 

Now  consider  the  scattering  effects  of  the  complete  PLACES  bar¬ 
ium  cloud  on  the  propagation  of  a  10.23  hfcps  PN  spread  spectrum  signal. 
Figure  2-23  shows  the  output  of  the  receiver  code  correlator  for  the  total 
PLACES  barium  cloud  model  (deterministic  cloud  plus  striations).  The 
entire  50  km  MPS  grid  is  shown  in  the  figure  by  displaying  the  code 
correlator  output  at  101  equally  spaced  locations  along  the  x-direction  of 
the  MPS  grid.  Note  the  increased  time  delay  at  the  central  portion  of  the 
cloud  caused  by  mean  electron  density  effects  as  modeled  by  the  large 
Gaussian  lens.  Since  all  the  striations  are  confined  to  the  central 
region  occupied  by  the  barium  cloud  and  each  striation  acts  roughly  as  a 
single  Gaussian  lens,  it  is  easy  to  explain  the  appearance  of  the  undis¬ 
turbed  region  of  delay  space  directly  beneath  the  center  of  the  cloud  in 
terms  of  the  overlapping  of  a  finite  number  of  scattering  patterns  from 
strong  Gaussian  lenses.  The  effect  of  a  single  strong  Gaussian  lens  is 
discussed  in  the  preceding  s'jboection. 
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Figure  2-22.  Envelope  of  time-domain  signal  after  propagation 
through  deterministic  barium  cloud. 
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Figure  2-23.  Envelope  of  received  time-domain  signal  after 
propagation  through  an  ionized*  striated 
barium  cloud. 
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Only  101  of  the  16384  MRS  grid  points  are  displayed  in  Figure 
2-23  so  that  there  appears  to  be  little  continuity  between  successive 
individual  code  correlation  functions.  This  strictly  visual  shortcoming 
is  remedied  in  the  next  figure. 

Figure  2-24  shows  the  amplitude  of  the  code  correlator  output 
over  the  range  from  26.50  km  to  27.72  km  corresponding  to  the  scale  shown 
on  the  plot  of  the  bariim  cloud  phase.  Figure  2-21.  The  cloud  center  is 
located  at  25  km  on  the  MRS  grid  and  so  the  1.22  km  region  shown  in  Figure 
2-24  is  located  just  past  the  barium  cloud  center.  Note  how  the  overall 
signal  time  delay  spread  is  increasing  as  the  distance  increases.  This  is 
due  to  the  overlapping  of  the  scattering  patterns  from  the  individual 
striations  which  make  up  the  random  part  of  the  barium  cloud. 

Results  from  the  actual  PLACES  experiment  were  provided  by  Dr. 
James  Marshall  (Marshall,  1982)  and  are  shown  in  Figure  2-25.  The  actual 
experiment  and  these  results  are  described  in  the  above  reference  and  will 
not  be  discussed  here.  However,  it  is  noted  that  all  the  essential  details 
of  the  experimental  results  were  accurately  predicted  in  advance  by  use  of 
numerical  MRS  propagation  simulation  techniques.  A  comparison  of  Figures 
2-23  and  2-25  shows  that  the  major  scattering  features  including  the 
effects  of  the  mean  barium  cloud  and  of  the  combined  scattering  pattern  of 
the  irregularities  is  accurately  predicted.  The  total  predicted  time 
delay  of  around  30  chips  is  also  in  agreement  in  both  figures.  These 
results  indicate  that  the  MRS  simulation  can  be  used  to  accurately  deter¬ 
mine  the  effects  of  variations  in  the  quantities  that  describe  the  barium 
cloud  including  mean  ionization  and  electron  density  PSD. 

2.6  HOMOGENEOUS  FLUCTUATIONS  -  r  COMPARED  TO  THEORY 

In  this  section  numerical  results  from  the  MRS  propagation  code 
are  compared  to  theoretical  approximations  for  the  two-frequency  mutual 
coherence  function  and  for  the  mean  time  delay  and  time  delay  jitter.  Here 
the  theoretical  results  for  the  two-position,  two-frequency  mutual  coher¬ 
ence  function  are  developed  in  Section  4  for  the  case  of  a  homogeneous 
thick  scattering  layer.  As  explained  in  Section  4,  the  theoretical  results 
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Figure  2-24.  Close-up  of  central  portion  of  Figure  2-23. 
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-25.  Received  PN  code  correlator  output,  PLACES  experiment 
first  beacon  rocket. 
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are  derived  on  the  basis  of  a  strong  scattering  assumption  using  a  quad¬ 
ratic  expansion  for  the  phase  structure-function  (Sreenivasiah,  et  a1, 
1976;  Sreenivasiah  and  Ishimaru,  1979). 

First  a  detailed  comparison  of  the  analytic  two-position,  two- 
frequency  mutual  coherence  function  with  several  interesting  MRS  cases  is 
presented.  Then  the  results  from  a  large  number  of  MRS  cases  are  compared 
to  theoretical  results  on  the  basis  of  the  e"^  point  of  r(Ax=0,  Af ) . 

Consider  the  case  of  plane  wave  propagation  of  a  wide  bandwidth 
signal  through  a  strong  scattering  layer.  The  carrier  frequency  is  2.25 
GHz,  the  phase  standard  deviation  is  350  radians  with  a  K"^  power-law  RSD 
between  the  outer  scale  of  1  km  and  the  inner  scale  of  10  m.  The  layer 
thickness  is  1000  km  and,  after  passage  through  the  layer,  the  disturbed 
plane  wave  propagates  an  additional  1000  km  in  free-space.  Ten  equally 
spaced  single  phase-screens  are  used  in  the  MRS  simulation  to  represent 
the  thick  layer.  Each  grid  consists  of  16384  cells  or  points  which  repre¬ 
sent  a  spatial  extent  of  10  km. 

A  total  of  32  discrete  frequencies  are  used  to  represent  the 
signal  spectrum  over  a  two-sided  bandwidth  of  4.096  MHz.  This  bandwidth 
corresponds  to  the  null-to-null  bandwidth  of  a  pseudo-noise  (RN)  spread 
spectrum  code  with  a  chip  rate  of  2.048x10^  sec“^  and  with  a  chip  duration 
of  488.3  nsec.  Equivalently  this  null-to-null  bandwidth  may  be  used  to 
represent  a  single  modulated  triangular  pulse  band-limited  to  4.096  MHz 
with  the  form 

m{T)  =  1  -  (t|/Tc,  |t|<Tc  (2-64) 

where  T^.  is  488.3  nsec. 

Figures  2-26  through  2-31  compare  MRS  results  for  the  two- 
frequency  mutual  coherence  function  r(Ax,Af)  with  the  analytic  results 
given  by  Equation  4-72.  Figure  2-26  shows  the  magnitude  of  r  as  a  func¬ 
tion  of  Af  for  parametric  values  of  ax  which  correspond  to  multiples  of 
the  distance  between  cells  on  the  numerical  MRS  grid.  The  agreement 
between  the  MRS  and  the  analytic  results  is  seen  to  be  quite  good  for  this 
case  with  some  bias  toward  a  higher  calculated  correlation  than  predicted 
analytical ly. 
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Figure  2-28.  Comparison  of  phase  of  r,  single  MRS  realization 
versus  theory,  ax  =  0.61  m. 


Figure  2-29.  Comparison  of  phase  of  r,  single  MRS  realization 
versus  theory,  ax  =  1.22  m. 
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Figure  2-30.  Comparison  of  phase  of  ft  single  MRS  realization 
versus  theory,  ax  =  1.83  m. 


Figure  2-31.  Comparison  of  phase  of  r,  single  statistical  realization 
versus  theory,  ax  =  2.44  m. 
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One  of  the  assumptions  involved  in  the  strong  scattering  theory 
discussed  in  Section  4  is  that  the  bandwidth  of  interest  is  small  relative 
to  the  center  frequency  so  that  the  refractive  index  may  be  taken  to  be 
independent  of  frequency  within  the  frequency  band  of  interest.  As  a 
result  of  this  assumption  the  analytic  two-frequency  mutual  coherence 
function  is  symmetric  about  the  center  frequency.  As  can  be  seen,  the  MPS 
result  for  r  is  not  synmetric  about  zero  Af  but  is  close  enough  to  justify 
the  analytic  approximation.  Since  the  index-of-refraction  is  a  function 
of  frequency,  synmetry  about  zero  Af  does  not  actually  occur;  but,  as  is 
shown  in  the  figure,  this  symmetry  can  be  a  valid  approximation. 

Figures  2-27  through  2-31  compare  the  MPS  results  for  the  phase 
of  r(Ax,Af)  with  the  analytic  results.  The  numerical  MPS  results  are  shown 
at  the  same  ax  values  as  shown  in  the  amplitude  results  just  compared. 
Again  the  comparison  is  quite  good  especially  in  the  most  important  region 
near  the  carrier  frequency  where  Af  is  small.  This  region  is  important  to 
time-domain  calculations  since  in  any  modulation  format  most  power  is 
transmitted  at  or  very  close  to  the  carrier  frequency  with  a  gradual 
decrease  in  transmitted  power  spectral  density  for  frequencies  farther 
displaced  from  the  carrier  frequency. 

Figure  2-32  shows  an  example  of  the  time-domain  results  for  this 
homogeneous  MPS  case.  The  amplitude  of  the  received  triangular  waveform, 
|e(x,Zr,T)|  is  plotted.  101  curves  are  shown  in  the  figure  as  a 
three-dimensional  plot  with  curves  for  successive  x-values  plotted  behind 
those  for  prior  x-values.  Each  individual  curve  represents  the  received 
time-domain  signal  observed  at  an  individual  MPS  cell,  separated  from  the 
next  by  0.61035  m.  This  plot  shows  only  the  first  101  cells  of  the  total 
of  16384  generated  but  is  quite  useful  to  display  the  time-domain  results 
because  of  the  homogeneous  geometry. 

Figure  2-33  shows  the  average  pulse  shape  defined  as  the  average 
over  X  of  the  amplitude  of  the  received  complex  envelope 

<|e(x,Z|.,T)|>j^ 

where  the  subscript  on  the  angle  brackets  indicates  that  the  averaging  is 
performed  over  x.  The  curve  shown  is  obtained  by  averaging  only  the  first 
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TIME  DELAY  (CHIPS) 


Figure  2-32.  Envelope  of  received  time-domain  waveform— first 
101  x-values  from  M»S  simulation. 


Figure  2-33.  Average  received  time-domain  envelope--MPS  signal 
realization. 
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2048  points  of  the  total  16384  point  MPS  realization.  This  case  is  dis¬ 
cussed  in  more  detail  in  Section  3  where  a  generated  statistical  signal 
realization  is  compared  to  the  MPS  results  shown  here. 

Figures  2-34  to  2-37  show  additional  comparisons  of  MPS  calcu¬ 
lations  of  the  magnitude  of  r  with  the  analytic,  strong  scattering  results 
of  Section  4.  The  cases  shown  are  calculated  for  a  satellite  communica¬ 
tion  system  with  an  80  MHz  bandwidth  spread  spectrim  signal  at  a  carrier 
frequency  of  7.5  GHz.  All  four  cases  have  a  one-dimensional  phase 
power  spectrum. 

In  Figures  2-34  to  2-36,  the  MPS  grid  length  is  30  km  with 
16384  grid  points,  the  striation  layer  thickness,  L,  is  14000  km  with  a 
free-space  propagation  distance  z^-L  of  1000  km  (see  Figure  2-2).  In 
all  three  figures  the  inner  scale  is  10  m.  In  Figure  2-34  is  300 
radians  and  the  outer  scale  is  3  km.  Figure  2-35  shows  results  for  a 
similar  case  to  Figure  2-34  but  with  a  smaller  phase  standard  deviation, 
a^,  of  150  radians.  In  the  comparison  shown  in  Figure  2-36,  is  100 
radians  and  the  outer  scale  is  1  km.  The  results  shown  in  Figure  2-37  are 
for  a  layer  thickness  of  1000  km  centered  1500  km  from  the  receiver 
plane.  Here  the  outer  scale  is  10  km,  the  inner  scale  is  10  m  and  the  MPS 
grid  length  is  100  km  with  16384  grid  points. 

The  agreement  between  the  MPS  results  and  the  analytic  results 
is  quite  good  for  all  these  cases  with  the  largest  difference,  on  the 
order  of  20  to  30  percent  near  the  maximmi  frequency  deviation  tf  of  40 
MHz  shown  for  Figure  2-37.  However,  this  case  shows  evidence  of  signal 
focusing  with  a  measured  S^  scintillation  index  of  1.15  characteristic  of 
non-saturated  scintillation;  as  such  the  strong  scattering  assumption 
involved  in  using  a  two  term  expansion  for  the  phase  structure-function  is 
inval id . 


In  order  to  make  possible  a  comparison  of  a  large  nunber  of  MPS 
calculations  with  analytic  results,  one  feature  of  r(Ax=0,f)  is  selected 
for  comparison  purposes.  Figure  2-38  shows  the  comparison  of  the  selected 
measure  f^.  The  quantity  f^  is  defined  such  that,  if  the  magnitude  of 
r(Af,Ax=0)  is  larger  than  1/e  for  ^  Af  £  fmax*  value  of 
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Figure  2-34.  Comparison  of  amplitude  of  r>  single 
realization  versus  theory. 
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Figure  2-35.  Comparison  of  amplitude  of  r,  single  MRS 
realization  versus  theory. 
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Figure  2-38.  Comparison  of  two-sided  bandwidth,  f^,  many  MPS  calculations 
versus  theory. 

is  fniax'^min-  ^  measiirf^  of  the  two- 

sided  bandwidth.  Table  1  is  a  list  of  the  two-sided  bandwidth  results 
presented  in  the  figure  and  shows  the  pertinent  MPS  parameters  that  are 
used.  In  all  cases  a  K"^  phase  PSD  is  used,  16384  points  comprise  the  MPS 
grid  and  the  inner  scale  is  fixed  at  10  m.  Lq  is  the  outer  scale,  L  is 
the  scattering  layer  thickness  and  z^-L  is  the  free-space  propagation 
distance  between  the  receiver  plane  and  the  nearest  edge  of  the  scattering 
medium.  Figure  2-2  depicts  the  MPS  geometry.  Because  of  the  large  nunber 
of  varying  quantities  that  are  used  here,  statistical  measures  that 
describe  the  overall  deviation  are  rather  meaningless.  Suffice  it  to  say 
that  the  results  are  in  quite  nood  agreement.  This  agreement  is  somewhat 
unexpected  because  of  the  apparently  severe  limitations  of  tiie  analytic 
calculation.  However,  for  the  most  part,  the  MPS  cases  described  here  do 
indeed  represent  strong  scattering  situations  and  seem  to  agree  quite  well 
with  the  theory. 
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Table  2-1 


Analytic  versus  numerical  MPS  comparison  of  1/e  points 
of  r{Ax=0,f). 


f 

0^ 

L 

L 

z  “L 

f 

0 

o 

r 

r 

MPS 

THEORY 

0.1  GHz 

33.0  rad 

0.39  km 

10  km 

100  km 

1.2  MHz 

1.2  MHz 

0.  15 

22.0 

0.39 

10 

100 

6.2 

5.5 

0.20 

16.5 

0.39 

10 

100 

16.4 

14.3 

1.38 

300 

5 

1000 

1000 

9.6 

8.2 

1.38 

500 

5 

1000 

1000 

4.4 

3.9 

2.25 

350 

1 

1000 

1000 

2.3 

1.8 

7.5 

100 

1 

14000 

1000 

46.5 

42.4 

7.5 

150 

3 

14000 

1000 

84.0 

82.5 

7.5 

300 

3 

14000 

1000 

31.5 

28.4 

7.5 

300 

3* 

14000 

1000 

39.0 

37.8 

7.5 

300 

10 

1000 

1000 

82.0 

70.2 

*  For  this  case  the  Inner  scale  Is  100  m.  For  all  other  cases  here  the  Inner  scale  Is  10  m. 


2.7  COMPARISON  OF  <t>  AND  TO  THEORY 

An  important  measure  of  the  effects  of  scintillation  on  propaga¬ 
tion  of  wide  bandwidth  signals  is  the  mean  time  delay  and  the  time  delay 
standard  deviation.  These  quantities  are  direct  measures  of  the  effects 
of  time  delay  jitter  in  the  mean  sense  and  also  serve  as  useful  character¬ 
izations  of  the  propagation  environment. 
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where  M(v)  is  the  transmitted  signal  spectrum  and  U(x,z^,v+u^)  is  a 
realization  of  the  solution  to  the  parabolic  wave  equation  at  a  frequency 
v+uq.  It  is  apparent  that  U(x,Zj.,(oq)  is  the  received  signal  amplitude  and 
phase  at  the  carrier  frequency.  In  the  comparisons  to  be  shown  in  this 
section,  the  transmitted  signal  is  always  a  band-limited  triangular  wave 
with  a  spectrin  given  by  Equation  2-54. 

The  ensemble  average  defined  in  Equations  2-65  to  2-67  is 
obtained  in  the  case  of  a  ninerical  MPS  calculation  by  averaging  the 
results  over  the  MPS  calculation  grid,  which  here  means  averaging  over  all 
x-values.  As  can  be  seen,  <t>  and  are  both  functions  of  the  receiver 
location  z^,  but  this  dependence  is  omitted  for  notational  convenience. 

The  numerical  MPS  results  are  to  be  compared  to  results  derived 
following  the  moment  method  formulation  of  Yeh  and  Liu  (1977),  with  the 
single  exception  that  the  ionization  irregularities  are  assumed  elongated 
infinitely  in  the  y-direction  so  as  to  maintain  the  fVS  geometry.  The 
moment  method  results  can  be  derived  as 
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(2-69) 


Xa2(2z^-L)|A2| 


a2  =  t2(0)  + 

T  ,.2 


_4  ,2  a2 

+  -i - A_  (12z2  -  16LZ^  +  6L2) 

24.2,2a2 

+  ^  (3z2  .  3Lz  +  l2)  (2-70) 

o  2  A  2  '  ' 

2ii^  ^o“o 

where  L  is  the  thickness  of  the  ionized  region  and  the  free-space 
propagation  distance  to  the  receiver  is  Zj.-L.  t2(0)  is  the  spread  of 
the  undisturbed  time-domain  waveform  which  here  is 


/  T^|rn(T)|^dT 

t2(0)=l!! -  (2-71) 

/|m(T)|^dT 


which  for  the  case  of  a  triangular  pulse 

m(T)  =  1  -  |t|/Tj.  ,  |t|<T^  (2-72) 

becomes 


This  result  is  a  slight  overestimate  of  the  actual  undisturbed  waveform 
spread  because  it  assumes  (in  effect)  an  infinite  bandwidth.  However,  the 
approximation  should  be  reasonably  accurate  since  most  of  the  power  is 
contained  within  the  first  nul 1-to-f irst  null  bandwidth  used  here. 
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The  quantities  Aq  and  A2  are  defined  by  Equations  4-44  and  4-45 
in  Section  4.1.3.  A^  is  the  coefficient  of  the  next  term  in  the  expansion 
of  the  phase  structure-function  and  is  defined  similarly  to  Aq  and  A2.  For 
the  case  of  a  K"^  Bessel  function  phase  PSD  with  outer  scale  Lq  and 
inner  scale  Ji,, 


Aj/A,  =  «n(t,/LJ/2L2 

(2-74) 

A./Ao  .  1/(8a?L2) 

(2-75) 

Note  that  the  result  for  <t>  derived  using  the  moment  method  of  Yeh  and 
Liu  is  identical  to  the  thick  layer  result  derived  in  Section  4.4.1  and 
presented  in  Equations  4-122  and  4-125.  The  second  and  third  term  of  the 
result  for  is,  in  the  thin  layer  approximation  {L=0),  identical  to  the 
results  of  Equation  4-135  for  a  thin  phase-screen.  The  A^  term  is,  of 
course,  not  included  in  the  theoretical  results  of  Section  4  where  only  a 
two  term  expansion  of  the  phase  structure-function  is  used. 

Figures  2-39  and  2-40  provide  a  comparison  of  the  MPS  results 
with  the  theoretical  results  obtained  following  the  moment  method  calcula¬ 
tion  of  Yeh  and  Liu.  The  results  shown  here  encompass  all  MPS  calculations 
to  date  and  involve  five  values  of  outer  scale,  three  values  of  inner 
scale,  a  range  of  carrier  frequencies  from  100  MHz  to  7.5  GHz  and  several 
different  geometries.  In  all  cases  a  K"^  phase  power  spectrim  is  assumed. 
Table  2  summarizes  the  important  calculation  parameters  used  and  gives  the 
values  plotted. 

Fairly  good  agreement  is  noted  in  all  cases  with  a  tendency  for 
the  analytic  results  to  be  somewhat  larger  than  the  MPS  calculations.  The 
100  MHz  carrier  frequency  calculation,  listed  first  in  Table  2,  is  omitted 
from  the  figure  because  it  would  have  required  a  much  larger  scale  for  the 
gain  of  one  additional  point.  The  lack  of  agreement  at  the  higher  values 
of  <T>  and  may  be  caused  by  the  impact  of  the  narrow  bandwidth  approxi¬ 
mation  used  in  the  theoretical  calculation  which  becomes  less  appropriate 
with  increasing  scintillation  severity  or  increasing  bandwidth. 
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Figure  2-39.  Comparison  of  mean  time  delay,  <t>,  many 
MRS  realizations  versus  theory. 


Figure  2-40.  Comparison  of  time  delay  jitter,  o^.,  many 
MRS  realizations  versus  theory. 


Table  2-2.  Comparison  of  analytic  and  nunerical  HPS  results  for 
mean  time  delay  and  time  delay  jitter. 


f 

T 

0, 

L 

L 

z  -L 

<T>/T 

0  /T 

0 

c 

<t> 

O 

1 

r 

c 

T 

c 

MRS 

THEORY 

MPS 

THEORY 

0.1  GHz 

97.5  nsec 

33.0 

0.39  kffl 

10  ni 

10  km 

100  km 

8.0 

9.5 

9.6 

14.2 

0.15 

97.6 

22.0 

0.39 

10 

10 

100 

2.0 

1.9 

2.4 

3.1 

0.20 

97.6 

16.5 

0.39 

10 

10 

100 

0.6 

0.7 

1.1 

1.1 

1.38 

97.6 

100 

3 

10 

14000 

1000 

0.85 

1.1 

1.1 

1.V 

1.38 

97.6 

100 

5 

10 

1000 

1000 

0.06 

0.08 

0.34 

0.37 

1.38 

97.6 

300 

5 

10 

1000 

1000 

0.65 

0.7 

0.8  7 

1.1 

1.38 

97.6 

500 

5 

10 

1000 

1000 

1.5 

1.8 

2.0 

2.6 

2.25 

488.3 

350 

1 

10 

1000 

1000 

1.1 

1.3 

1.4 

1.8 

7.5 

25.0 

30 

1 

10 

14000 

1000 

0.04 

0.09 

0.32 

0.36 

7.5 

25.0 

100 

1 

10 

14000 

1000 

1.0 

1.0 

1.1 

1.3 

7.5 

25.0 

100 

3 

10 

14000 

1000 

0.03 

0.14 

0.37 

0.39 

7.5 

25.0 

150 

3 

10 

14000 

1000 

0. 29 

0.31 

0.48 

0. 54 

7.5 

25.0 

300 

3 

10 

14000 

1000 

1.1 

1.3 

1.1 

1.6 

7.5 

25.0 

300 

3 

30 

14000 

1000 

0.97 

1.0 

1.3 

1.3 

7.5 

25.0 

300 

3 

100 

14000 

1000 

0.84 

0.76 

0.87 

0.97 
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SECTION  3 

STATISTICAL  SIGNAL  GENERATION 


3.1  STATISTICAL  SIGNAL  REALIZATIONS 

In  this  section  a  second  technique  is  described  to  generate 
realizations  of  wide  bandwidth  signals  after  propagation  through  the  iono¬ 
sphere.  This  technique  is  referred  to  as  statistical  signal  generation 
and  relies  on  the  properties  of  the  solution  to  the  parabolic  wave  equat¬ 
ion  under  strong  scattering  conditions.  This  reliance  renders  statistical 
signal  generation  less  general  than  the  MRS  propagation  algorithm  discus¬ 
sed  in  Section  2. 

The  basic  formalism  to  generate  statistical  signal  realizations 
was  developed  by  Wittwer  (1979  and  1980)  for  the  case  of  isotropic  irregu¬ 
larities.  This  section  utilizes  much  of  the  initial  formalism  developed 
at  that  time. 

This  statistical  signal  generation  technique  is  based  on  the 
solution  for  the  two-position,  two-frequency  mutual  coherence  function 
r(Ax,Au))  for  spherical  waves  described  in  Section  4.  In  that  section  the 
solution  for  r  in  the  strong  scattering  limit  is  discussed  in  detail  and 
its  Fourier  transform  is  presented.  This  Fourier  transform,  S{K,t),  is 
known  as  the  generalized  power  spectrum  (Wittwer,  1979)  and  its  knowledge 
is  a  requirement  for  statistical  signal  generation.  The  formulation  from 
Section  4  for  S(K,t)  in  the  thin  phase-screen  limit  is  used  in  this  sec¬ 
tion  as  the  basis  of  statistical  signal  generation.  The  thin  phase-screen 
approximation  is  necessary  to  provide  tractable  expressions  of  S(K,t)  for 
numerical  evaluation.  As  shown  in  detail  in  the  next  section,  the  thin 
phase-screen  approximation  is  highly  accurate  in  all  cases  as  an  approxi¬ 
mation  to  the  geometry  of  a  satellite  communication  link. 

Hence,  the  major  limitation  on  the  accuracy  of  statistical  sig¬ 
nal  generation  is  the  strong  scattering  assumption  which  provides  the 
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basis  for  the  simplified  solution  for  the  two-frequency  mutual  coherence 
function.  In  the  case  of  weak  scattering,  the  MPS  methods  of  Section  2 
may  easily  be  applied.  For  strong  scattering  cases,  often  of  interest  in 
the  study  of  the  effects  of  nuclear  detonations  on  communications  links, 
the  statistical  signal  generation  technique  is  more  efficient. 

Although  statistical  signal  realizations  are  useful  only  for  the 
case  of  strong  scattering,  their  generation  requires  only  a  fraction  of 
the  computer  resources  required  for  signal  generation  by  wide  bandwidth 
MPS  calculations. 

Statistical  signal  generation  is  applicable  to  the  case  of  elon¬ 
gated  irregularities  as  handled  by  the  MPS  simulation  or  to  the  more 
general  case  of  isotropic  irregularities.  Here  the  formulation  is  de¬ 
veloped  for  both  cases.  Two  detailed  examples  are  presented  for  the  case 
of  elongated  irregularities  and  are  compared  to  exact  MPS  calculations 
that  are  described  in  Section  2. 

3.2  FORMULATION 

The  scalar  Helmholtz  wave  equation  may  be  written  as 

(v2+k2e^(r))  E{r,a,)  =  0  (3-1) 


where  kg  is  the  free-space  wavenumber,  ty  is  the  relative  permittivity  and 
E(T,a))  is  the  electric  field  at  the  angular  frequency  w.  The  exp(i(ji)t) 
time  dependence  has  been  suppressed.  If  the  permittivity,  e:^(T’),  is  a  ran¬ 
dom  quantity  with  an  average  value  plus  a  small,  zero-mean  random  part, 
then 

tjr)  =  <£(r)>  /l*  ££<rL\  (3-2) 

"  \  <4n>l 


Now  if  the  dielectric  permittivity  fluctuations  are  caused  by  electron 
density  fluctuations,  then 


<e(r)> 


<Ng(F)> 


n 

c 


1  - 


(3-3) 
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(3-4) 


Ae(r)  =  - 


AN„(r) 


ANg(r) 

<''e> 


where  <Ng>  is  the  mean  electron  density,  n^  is  the  critical  electron 
density  and  wp  is  the  plasma  frequency.  These  three  quantities  are 
related  by 


0)2  =  4Trc2r  <N  > 
p  e  e 


(3-5) 


n 


c 


(3-6) 


where  rg  is  the  classical  electron  radius  (2.82xl0"^5m)  and  c  is  the 
speed  of  light  in  free-space.  Combining  Equations  3-1  to  3-4,  the  scalar 
Helmholtz  equation  for  a  mediim  subject  to  electron  density  fluctuations 
combined  with  mean  ionization  is  obtained 


V^E  +  k2 


<N  > 
e 


1  -  lOp/o)^ 


E  =  0 


(3-7) 


Now  k  is  the  mean  wavenumber  for  propagation  in  an  ionized  medium 


k  =  “  <e>172  =  w  (1  _  (3_g) 

c  c  p 

where  w  is  the  radian  frequency  of  the  transmitted  waveform.  The  mean 
wavenumber  k  is  a  function  of  propagation  distance  z,  but  this  dependence 
is  omitted  for  convenience  of  notation  in  the  following. 

Figure  3-1  shows  the  geometry  of  the  problem  where  a  spherical 
wave  originates  from  a  transmitter  located  at  -z^  and  propagates  in 
free-space  in  the  positive  z  direction,  until  it  is  incident  on  a  layer  of 
ionization  irregularities  of  thickness  L.  After  emerging  from  the  layer 
of  irregularities  at  z  =  L,  the  wave  then  propagates  in  free-space  to  a 
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receiver  located  in  the  plane  z  =  z^.  This  geometry  is  chosen  to  model 
transionospheric  satellite  communications  links  through  ionization  pro¬ 
duced  by  natural  causes,  bariuti  releases,  or  high  altitude  nuclear 
detonations. 


is  the  transverse  Laplacian. 
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The  complex  amplitude  U  can  vary  markedly  over  distances  no 
smaller  than  the  scale  size  of  the  electron  density  fluctuations  that 
is,  the  worst  case  variation  of  U  in  the  direction  of  propagation  can  be 
characterized  by  U-  exp(-z2/£?,„  )  where  jl >  l.  Thus  the  second  deri- 

rt  niclX  fiy  OiqX  S 

vative  a'^U/az^  is  of  the  order  of  .  On  the  other  hand,  the  third 

ni3x 

term  of  Equation  3-10  is  of  the  order  of  U/X£„,„.  Therefore,  for  X«Z^  , 

iTidx  rnsix 

or  X«Jl^  the  second  term  of  Equation  3-10  may  be  ignored  with  respect  to 

the  third  term  and  the  parabolic  wave  equation  is  obtained  as: 


9  ..1  <N  >  m2 

vfu  -  2ik  -  k2 — ! -  U  =  0  (3-11) 

az  1  -  0)2/ 0)2 

The  solution  to  tne  parabolic  wave  equation  at  the  receiver 
plane  is  UCp,z^,o)}.  As  discussed  in  Section  2,  if  the  ionization  irregu¬ 
larities  are  infinitely  elongated  in  the  y-direction,  there  is  no  y  varia¬ 
tion.  Then  the  multiple  phase-screen  propagation  technique  can  be  used  to 
obtain  Monte  Carlo  "realizations"  of  U(x,z^,o)),  the  solution  to  the 
parabolic  wave  equation. 

3.2.1  Wide  Bandwidth  Signals 

Now,  if  instead  of  a  monochromatic  wave,  the  transmitted  signal 
is  composed  of  a  spectrim  of  Fourier  components,  S(u)),  then  the  time- 
domain  signal  at  the  receiver  plane  is 


v(p,z^,t)  =L  jfs(u))  U(T.z„,u))  e'®^“^e'“^d, 

2ll  _oo 


(3-12) 


where 


g(  w) 


k(z' )  dz’ 


“  /  (1  -  u)2(z')/a)2]'^^  dz' 

C  .fj  P 


(3-13) 


91 


The  propagation  of  wide  bandwidth,  modulated  waveforms  is  of 
interest  in  this  work.  Thus  let  the  real  transmitted  waveform  be 
expressed  as 


s(t)  =  Re  I m(t)e  j 

=  ^m(t)e  °  +  ^m*(t)e  °  (3-14) 

where  Wq  is  the  carrier  angular  frequency  and  m(t)  is  the  modulation 
waveform.  For  a  pseudo-noise  (PN)  code,  m(t)  consists  of  a  long  sequence 
of  pseudo-random  polarity  changes  (binary  PSK-modulated  signal).  The 
frequency-domain  representation  of  the  transmitted  waveform  is  given  by 

S(a))  =  /  s(t)e"''“^dt 


i( u)  -(o)t 


=  I-  f  m(t)e  °  dt  +  1  r  m*(t)e 
^  ^  2 


-  i(  U)  +(;is)t 


dt 


=  i 

2  0  2 


(3-15) 


whei  e  M(u))  is  the  Fourier  transform  o.  the  modulation  waveform: 


M(a))  =  f  m(t)e"''“^  dt 
-00 


(3-16) 


Combining  Equations  3-12  and  3-15,  the  received  waveform  v  can  be 
expressed  as 

v(p>z  ,t)  =  L-  f  )  U('p,z  ,(jj)  e^®^“^du) 

^  Att  0  ^ 


f  M*(-u)-u,q)  U(p,z^,a))  e^®^“^da) 


47T 


(3-17) 


(3-23) 


i 


In  the  second  term  of  Equation  3-17,  expand  e(o))  in  a  Taylor  series  about 

-wq- 


(to+oj  9"(-U-) 

6_(a))  =  +  (ojtcoQle' (-Uq)  +  - - +  ... 

(-to-o)  )2  e"{u  ) 

=  -  e{uQ)  -  (-uj-toQ)e'(u)iQ)  -  — — ^ — — - +  ...  (3-24) 

For  notational  simplicity,  let  Qim^)  =  0'(ujj)  =  6^  and  so  forth.  Re 

taining  the  first  four  terms  in  the  Taylor  series,  the  received  waveform 
can  be  written  as: 

00 

v(‘p,z^»t)  '  —  f  U(?,z^,aj)  exp {i 0jj+i¥( u-to^)+iwt }d(0 

«»oo 


1 

/  M*(-0)-(0q)  U(p,Z^,0))  exp{-ieQ-i¥(-a>-u)Q)-iu)t}da> 

(3-25) 

where 

^(v)  =  v0'  +  -  v^e"  +  -  v^e'"  (3-26) 

0  2^6^ 

In  the  first  integral  let  v  =  u-Up  and  in  the  second  integral  use  the  sub¬ 
stitution  V  =  -oj-oJq.  Then  Equation  3-25  for  the  received  waveform  becomes 

i  ( uot+  Sq )  r  _ 

v(p,z^.t)  = — - -  J  M(v)  U(-p,z^,vtwQ)  exp(i0(v)+iv(t-tjj)}dv 

-00 


g-i(uot+eo) 

4it 


/ 


M*(v)  U(p,z^,-v-u)q)  exp{-ie(v)-iv(t-t^)}dv 


(3-27) 


C\A 


+ 


where  Qq  can  be  seen  to  be  the  total  mean  phase-shift  (caused  by  disperr, 
Sion  alone)  at  the  carrier  frequency  and  t^j  is  the  mean  group  time  delay 
(i.e.,  the  group  delay  at  the  carrier  frequency): 


(3-33) 


I  : 

: 

oj  2itcr  fL 

0-  a  -  -2.  R  +• - r_L  (radians) 

°  C 


2ncrgN^ 


(seconds) 


where 


R  »  Zt  + 


z 


r 


<N^(z')>  dz' 
e 


(3-34) 


R  is  the  total  propagation  distance  and  is  known  as  the  total 
electron  content. 


The  well  known  Doppler  shift  due  to  range  and  ionization  varia¬ 
tions  is 

f  si _ 0 

27T  dt 


c 


ere  N, 


(3-35) 


The  first  terms  of  Equations  3-33  and  3-34  are  simply  the  free-space 
phase-shift  and  time  delay  which  are  proportional  to  the  total  propagation 
distance  R.  Variations  in  this  distance  with  respect  to  time  produce  the 
well  known  Doppler  shift  and  effects  of  time  delay.  The  second  term  of 
Equations  3-33  to  3-35  represents  contributions  to  propagation  phase,  time 
delay  and  Doppler  shift  caused  by^mean  ionization.  Note  that  increasing 
total  electron  content  (positive  Nj)  causes  positive  Doppler  shift  and 
increasing  time  delay.  This  is  opposite  to  the  effect  of  range  variation, 
where  increasing  time  delay  is  associated  with  negative  Doppler  shift. 
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It  is  evident  from  the  parabolic  wave  equation,  3-10,  and  from 
the  definition  of  k,  Equation  3-8,  that 


=  U*(p,z^,-u)) 


With  this  fact  in  mind,  it  is  seen  that  the  two  terms  of  Equation  3-27  are 
complex  conjugates  of  one  another  so  that  the  received  signal  is  real,  as 
it  should  be.  This  waveform  can  be  written  in  terms  of  a  complex 
modulation  envelope, 


v(p,2^,t)  =  1  e(p,z^,t-t^)e 


=  Re  je(p,z^,t-tj)e  °  °  j  (3-36) 

where  the  received  complex  envelope  is 

09 

S  Hv)  U(^,z^,v+u)Q)e<p(i'iJ( v)+i vtjdv  (3-37) 

c7! 

Note  that  Equation  3-36  has  the  same  general  form  as  Equation  3-14  -  a 
wave  modulated  on  a  carrier  at  angular  frequency  uq  with  an  additional 
phase-shift  and  time  delay  as  given  by  Equations  3-28  and  3-29. 

In  Equation  3-27  the  received  complex  envelope  e(^,Zf.,T)  is 
given  as  the  integral  over  frequency  of  the  Fourier  transform  of  the 
transmitted  complex  envelope  multiplied  by  the  solution  to  the  parabolic 
wave  equation  U.  Since  the  value  of  the  resultant  received  signal  v  is  a 
function  of  the  carrier  frequency  uq,  the  Uq  dependence  must  appear 
explicitly  in  Equation  3-37.  Although  the  uq  dependence  is  not  explic¬ 
itly  shown  in  writing  e(pr,zr,T),  the  received  complex  envc''')pe  does 
indeed  depend  upon  the  carrier  frequency  as  shown  explicitly  in  Equation 
3-37. 
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3.2.3  Impulse  Response  Function 

Now  let  us  assume  that  the  complex  received  electric  field 
e(?,z^,T)  represents  the  response  of  a  linear  system  to  some  input  signal 
m(T).  Then  in  general  one  can  write 

e(p,z^,T)  =7m(f)  h(p.z^.T-t')dt'  (3-38) 

which  simply  states  that  the  output  signal  is  the  convolution  of  the  input 
signal  with  the  impulse  response,  m(T)*h{T).  Now  with  the  following 
substitution  obtained  directly  from  Equation  3-16 

M{v)  =  7  m(t')e'^''^'dt' 

••00 

Equation  3-37  can  be  written  as 

e(p,z^.T)  »  i-  7'"(^') 

cIT  mm 

Now  comparison  of  Equations  3-38  and  3-39  yields  the  important 
result  for  the  impulse  response  function 

h(p,z^,T)  = J'U('p,Zj.,\H-u^)e^^^''^e^^’'dv  (3-40) 

ZlT  ••09 

Equation  3-40  gives  a  direct  result  for  the  impulse  response  of  an  ionized 
mediim  in  terms  of  the  frequency  domain  solution  to  the  parabolic  wave 
equation.  It  is  apparent  from  Equation  3-40  that  the  impulse  response  h 
does  depend  upon  the  carrier  frequency  wq  although  this  dependence  is 
not  explicitly  shown  in  writing  h. 

Notice  that  Equation  3-40  represents  a  Fourier  transform  rela¬ 
tionship  for  the  impulse  response  function.  The  transform  of  the  impulse 
response  function  is  composed  of  the  product  of  two  terms.  One  term, 
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U(p,z^iv+u)Q),  is  the  solution  to  the  parabolic  wave  equation  and  repre¬ 
sents  random  effects  due  to  fluctuating  plasma.  The  second  term 

is  the  transfer  function  of  a  smooth  ionized  medium  and  can  cause  pulse 
dispersion.  Since  the  two  transfer  functions  appear  in  a  product  relation¬ 
ship,  it  is  convenient  to  separate  the  generation  of  signal  structures  in¬ 
to  a  two  step  process.  The  first  step  consists  of  generating  realizations 
of  the  impulse  response  function  with  only  fluctuating  plasma  effects. 

The  second  step  adds  the  effect  of  the  gross  ionization. 

With  the  above  thought  in  mind,  let  us  consider  a  method  to  gen¬ 
erate  a  realization  of  the  impulse  response  function  for  the  case  where 
^(v)  »  0.  In  addition,  consider  for  the  moment  the  generation  of  the  im¬ 
pulse  response  function  for  the  case  where  there  is  no  y-variation.  This 
situation  corresponds  to  the  case  of  elongated  irregularities  considered 
in  Section  4  and  to  the  two-dimensional  MPS  propagation  simulation  dis¬ 
cussed  in  Section  2.  Now  the  impulse  response  function  satisfies  the 
equation 


h(x,z^,T)  ^  /  U(x,z^,v+u)Q)e^'’^dv  (3-41) 

The  following  initial  development  is  limited  to  this  tivo-dlmensional  case; 
however  we  shall  return  to  the  general  three-dimensional  case  later. 

3. 2. 3.1  Distinction  Between  Time  (t)  and  Time  Delay  (t) 


Recall  that  the  x  or  'p  dimension  corresponds  to  a  displacement 
on  the  receiver  plane  as  shown  in  Figure  3-1.  In  an  actual  physical  situa¬ 
tion,  the  x-variation  is  converted  into  variation  with  time  by  effective 
motion  of  the  propagation  path  through  the  irregularities.  Thus, 

t  =  ^  (sec)  (3-42) 

''eff 
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where  the  effective  velocity  is  a  function  of  many  system  and  environ¬ 
mental  parameters  including  wind  velocities,  platform  velocity,  plasma 
velocity  and  propagation  path  geometry.  This  effective  motion  may  be 
caused  by  motion  of  the  irregularities,  by  motion  of  the  line-of-sight 
through  the  irregularities  or  by  a  combination  of  both.  Note  that  the 
time  (t)  variation  of  the  fluctuations  is  assumed  to  be  very  slow  relative 
to  the  time  (t)  variation  of  the  modulated  waveform.  In  other  words,  the 
bandwidth  of  the  fluctuations  is  small  compared  to  the  bandwidth  of  the 
modulation.  This  consideration  is  important  since  the  solution  to  the 
parabolic  wave  equation  is  obtained  simultaneously  at  a  number  of  frequen¬ 
cies  over  the  signal  bandwidth.  For  the  solution  to  be  accurate,  the 
mediun  cannot  vary  over  the  time  duration  of  the  modulation  waveform. 

3.2.4  Statistics  of  a  Realization 

At  this  point  it  is  useful  to  consider  the  technique  discussed 
in  Section  2.2.4  to  generate  phase-screen  realizations  as  it  might  be 
applied  to  generate  realizations  of  the  complex  electric  field. 

First  a  pseudo-random  realization  of  the  Fourier  transform  of 
the  electric  field  is  chosen  as 


E{mAK)  =  rJS(mAK)L/2n)V2 


(3-43) 


A 

where  E  is  a  realization  of  the  Fourier  transform  of  the  electric  field,  S 
is  the  desired  power  spectral  density  (PSD),  L  is  the  length  of  the  grid 
and  r^,  is  a  complex  Gaussian  random  variable  such  that 


<r  r*>  =  5 
m  n  nm 


(3-44) 


The  electric  field  itself  is  then  obtained  by  taking  the  discrete  Fourier 
transform  of  Equation  3-43 


E(ux)  «  Y  E(mAK)e^^’'*^^'^AK 
m=0 


t=0,...N-l  (3-45) 
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The  most  important  properties  of  the  electric  field  are  immedi¬ 
ately  evident  from  the  preceding  brief  discussion.  First,  since  the 
electric  field  is  a  linear  combination  of  a  number  of  Gaussian  variatee, 
both  its  real  and  imaginary  parts  possess  Gaussian  probabilfty  distribu¬ 
tions.  Thus  the  realization  of  received  signal  amplitude  is  the  square 
root  of  the  sun  of  the  squares  of  the  inphase  and  quadrature  .phase  parts, 
and 


{(Real  E)2  +  (Imaginary 

possesses  a  Rayleigh  probability  distribution. 

In  the  following,  a  formulation  is  developed  to  generate  real¬ 
izations  of  the  received  electric  field  which  are  restricted  to  possess 
Rayleigh  amplitude  distributions.  Fortunately  this  limitation  is  not  at 
all  restrictive  for  the  case  of  a  nuclear  disturbance  since  large  spatial 
regions  can  experience  Rayleigh  fading  statistics. 

Finally  it  is  evident  from  Equations  3-43  and  3-44  that  the 
electric  field  power  spectrum  has  the  form 

<E(mAK)  E*(nAK)>  -  6^^  S(m6K)  (3-46) 

which  is  a  requirement  if  the  electric  field  is  to  be  stationary 
(Tatarskii,  1971).  Thus  the  resulting  autocorrelation  function  generated 
by  the  method  described  in  Section  2.2.4  is  dependent  only  on  the  distance 
between  measurement  points. 

3.3  INPULSE  RESPONSE  FUNCTION  GENERATION 

In  this  section  techniques  are  developed  to  generate  realiza¬ 
tions  of  the  impulse  response  function  h(x,Zp,T)  defined  in  Section 
3.2.3.  Once  this  function  is  available,  Equation  3-38  or  its  equivalent 
may  be  applied  to  obtain  a  realization  of  the  complex  received  electric 
field.  It  is  shown  here  that  it  is  possible  to  generate  the  Fourier 
transform  of  the  impulse  response  function  U(x,z^,v)  in  essentially  the 


101 


same  manner  as  realizations  of  the  phase  are  generated.  Note  that  the 
Independent  variable  pairs  v-t  and  x-K  are  used  to  unambiguously  denote 
Fourier  transform  pairs. 

3.3.1  Realizations  of  the  Parabolic  Wave  Equation  Solution 

In  this  section  a  realization  of  the  Impulse  response  function 
Is  obtained  by  first  generating  a  realization  of  Its  Fourier  transform 
U(x,v).  Thus  given  the  Fourier  transform  pair, 

U(x,v)  =  f  h(x,T)e“^'’'^dT  (3-47) 

«•» 

h(x,T)  »  4-  ru(x,v)e^'’^dv  (3-48) 

2it 

first  generate  a  realization  of  U(x,v)  as  Illustrated  In  the  following. 
Then  take  the  transform  to  obtain  a  realization  of  the  Impulse  response 
function  h(x,T).  Note  that  Equation  3-48  Is  Equation  3-41  with  the  Zp 
and  uq  dependence  removed  for  convenience.  Both  of  these  dependencies 
play  an  Important  role  In  the  results  but  are  unnecessary  complications  to 
the  formulation. 

Now  U(x,v)  can  be  written  as  the  Fourier  transform 

U(x,v)  -/  jrh(K.T)e^'^V^''%dT  (3-49) 

•  00  *00 

where  h(K,T)  is  recognized  as  the  Fourier  transform  of  the  Impulse 
response  function  with  respect  to  x.  h(K,T)  Is  not  the  Impulse  response 
function  which  Is  denoted  as  h(x,T).  In  the  discrete  case  Equation  3-49 
may  be  written 

U(u«.pw)  .  .Y  h(m;«.s4,)  sKat 

s**0 

a»0,...N-1;  p»0,...M-1  (3-50) 
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where  N  discrete  points  are  used  to  represent  a  total  distance  of  L  so 
that  the  point  spacing  ax  =  L/N.  By  the  Nyquist  sampling  theorenij  the 
minimum  spatial  frequency  is  aQ  =  1/L.  The  minimum  spatial  radian 
frequency  is  then  2it/L.  Similarly  T  is  the  total  time  duration  of  the  M 
samples  used  to  represent  the  delay  space  functional  dependence.  Thus  the 
delay  increment  is  at  =  T/M.  By  the  Nyquist  sampling  theorem  the  minimum 
frequency  is  Af  =  1/T,  or  the  minimum  radian  frequency  is  2it/T.  The 
total*  two-sided  spatial  bandwidth  is  B^  -  1/ax*  and  the  total* 
two-sided  delay  bandwidth  is  B^j  =  1/at.  For  later  convenience  also 
define  the  total  radian  bandwidth  as  By  »  2ttB(J.  With  the 
substitutions  just  discussed  Equation  3-50  becomes 


U()iAx*pAv)  =  I  I  h(mAK*SAT)  e  g-2itps/M 

m=0  s=0 

t»0*...N-l;  p=0*...M-l  (3-51) 


Now  choose  as  a  realization  of  the  Fourier  transform  the  quantity 


h(mAK*sAT) 


JL  S(mAK*SAT)  r, 


(3-52) 


where  S(K*t)  is  the  generalized  power  spectrum  given  by  either  Equation 
4-97  or  4-99  in  Section  4.  is  a  complex  random  variable  such  that 


■■ms  *  \  (9lm*‘92m><9ls*'92s) 


(3-53) 


where  the  g's  are  independent*  Gaussian  random  variables  with  zero  mean 
and  unity  variance.  With  this  definition  for  the  random  quantity  r^g  it 
is  easy  to  show  that 


ms  nr  mn  sr 


(3-54) 


ms 


(3-55) 


Now  it  must  be  shown  that  the  choice  for  the  Fourier  transform. 
Equation  3-52,  is  correct.  Let  us  calculate  the  autocorrelation  function 
of  U(x,v)  and  show  that  it  is  indeed  the  two-frequency  mutual  coherence 
function.  The  continuous  autocorrelation  function  of  U(x,v)  is  given  as 

B^j(x,v)  =  <U(x'+x,v'+v)U*(x' ,  v*  )> 

,  L  B,/2 

=  ^  f  f  U(x'+x,v'  +  v)U*(x',v')dx'dv'  (3-56) 

•■“v  *0  -b;;/2 

where  ergodicity  is  assumed  and  ensemble  averages  are  replaced  by  averages 
over  distance  and  frequency.  Now  the  preceding  equation  may  be  written 
discretely  as 

By(x,v)  =  B^j(jiAx,pAv) 

,  N-1  M-1 

=  ^  I  I  U(l<AX+m,qA\H-pAv)  U*(l<AX,qAv)  AXAv 
LB^  k=0  q=0 


£=0,...N-1;  psO,...M-l 


(3-57) 


From  Equation  3-51  the  Fourier  transform  of  U  may  be  substituted  into 
Equation  3-57  to  obtain 

,  N-1  M-1  N-1  M-1  N-1  M-1 

B^j(jiax,pav)  =  _  I  I  I  I  I  I 

k=0  q=0  m=0  s=0  n=0  r=0 

X  h(mAK,sAT)  exp{i2n(k+2)m/N  -  i2Ti(q+p)s/M} 

X  h*(nAK,rAT)  exp{-i2Tikn/N  +  i2Tiqr/M}  axav(aKat)2 


2*0,... N-1;  p=0,...M-l 


(3-58) 
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Now  substitute  the  realization  of  the  Fourier  transform  given  by  Equation 
3-52  into  Equation  3-58  to  obtain 


'  B^j()iax,pav) 


1  N-1  M-1  N-1  M-1 

lilt 

{2it)2  k=0  q=0  m=0  s=0 


N-1 


M-1 


X  {S(mAK»SAT)S*(nAK,rAT)}^/^  ’"ms’'*nr 
X  exp{i2it(k+ji)ni/N  -  i2Ti(q+p)s/M} 


X  exp{-i2irkn/N  +  i2TTqr/M}  axav(aKat)2 

Now  since 

JZMn  ,-i2itkn/N  _ 

)  6  6  N6„„ 

k=0 

and 


(3-59) 

(3-60) 


Y  e-^2Ttqs/M  g+i2iTqr/M  , 

q=0  * 


(3-61) 


Equation  3-59  may  easily  be  reduced  to 

B,,(jiax,pav)  =  I  I  S(mAK,SAT) 

^  (2tj)2  m=0  s=0 

X  exp{i2TTjjm/N  -  i2iips/M}  (aKAt)^ 


N-1  M-1 

=  I  I  S(mAK,SAT)  Ir  |2  expli2TT«m/N  -  i2Tips/M}  /KAt 
m=0  s=0  ' 

(3-62) 

It  is  apparent  that  a  Fourier  transform  relationship  exists  between  the 
correlation  function  Bj  and  the  function  S(mAK,SAT) |rn,sl Thus  the 
power  spectrum  of  the  single  generated  realization  U(aax,pav)  is 

S(mAK,SAT)  |r^5|2 
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i»«i»iiaaaTO'ga6Tgi^ 


where  S(K,t)  is  the  desired  generalized  power  spectrum  discussed  in 
Section  4  of  this  report. 

It  is  apparent  that  each  of  the  Fourier  components  of  the  power 
spectrum  is  independent  of  Fourier  components  at  other  values  of  K  and  t 
because  of  the  choice  of  the  random  number  r^^g.  Because  of  this  delta 
correlation  property  of  the  power  spectrum,  the  autocorrelation  function  ^ 
is  only  a  function  of  ax  and  av  as  it  should  be  in  accordance  with 
the  description  of  the  two-position,  two-frequency  mutual  coherence  func¬ 
tion,  r,  given  in  Section  4.  Since  the  random  numbers  rn,s  are  Gaussian 
variates,  |rms|^  is  a  chi-squared  variate  with  two  degrees  of  freedom. 

Since  <|rn)sr>  =  1  the  mean  value  of  the  power  spectrum  is  the  desired 
value.  In  other  words,  each  realization  of  U(jiax,pav)  is  dependent  upon  a 
sequence  of  random  numbers,  r^^,  and  thus  has  a  spectrum  which  is  not 
identically  equal  to  the  desired  spectrum.  However,  the  average  PSD  of  a 
large  number  of  such  realizations  is  given  as 

Avg  (  S(mAK,SAT)|rj^5|2}  =  S(mAK,sAT)  (3-63) 

which  is  exactly  the  result  desired. 

The  average  autocorrelation  function  determined  by  averaging 
over  a  number  of  numerically  generated  realizations  is  obtained  by  taking 
the  average 'Of  Equation  3-62 

N-1  M-1 

Avg  {B,,(  AAx,pAv) }  ^  I  I  S(mAK,sAT)  exp{i2n]!m/N-i2tips/M}AKAT 
^  m=0  s=0 

(3-64) 

This  equation  may  be  rewritten  in  its  continuous  form  as 

Avg{B^j(x,v)}  =7  7  S(K,t)  e"’’'''"  dKdT  =  r(x,v)  (3-65) 

•00  mOO 

which  is  identical  to  Equations  4-84  or  4-87  and  proves  that  the  formula 
given  in  Equation  3-52  for  the  Fourier  transform  of  a  single  realization 
is  indeed  correct  as  claimed. 
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At  this  point  it  has  been  proven  that  Equations  3-51  and  3-52 
may  be  used  to  generate  a  realization  of  the  solution  U(x,v)  to  the  para¬ 
bolic  wave  equation.  This  realization  has  r(x,v)  as  its  autocorrelation 
function  and  S{K,t),  the  generalized  power  spectrum,  as  the  Fourier  trans¬ 
form  of  r.  In  order  to  facilitate  numerical  generation  of  realizations, 
the  formulation  for  S(K,t)  has  been  specialized  with  little  loss  of  gener¬ 
ality  to  the  case  of  scattering  by  a  thin  phase-screen. 

The  equations  for  S(K,t)  are  given  in  Section  4  for  two  cases. 
Equation  4-97  is  valid  for  the  case  of  irregularities  that  are  infinitely 
elongated  perpendicular  to  the  direction  of  propagation.  Equation  4-99  is 
valid  for  isotropic  irregularities.  Appendix  A  is  devoted  to  numerical 
evaluation  of  the  two  equations  and  shows  a  number  of  plots  of  S(K,t)  for 
an  interesting  selection  of  parameters. 

The  specialization  to  the  thin  phase-screen  case  mentioned  above 
is  discussed  in  detail  in  Section  4  where  results  are  presented  for  both 
thick  and  thin  layers  to  show  the  effect  of  layer  thickness  on  the  S4 
scintillation  index,  the  mean  time  delay  and  the  time  delay  jitter.  In 
general,  for  satellite  communication  links  through  ionospheric  disturb¬ 
ances,  the  thin  phase-screen  approximation  is  quite  accurate  and  enjoys 
a  distinct  advantage  in  usefulness  over  the  thick  layer  solution. 

3.3.2  Realizations  of  the  Impulse  Response  Function 

At  this  point  a  method  has  been  described  to  generate  discrete 
realizations  of  U(x,v),  the  solution  to  the  parabolic  wave  equation, 
directly  from  the  computed  two-frequency  mutual  coherence  function.  In 
order  to  generate  realizations  of  the  impulse  response  function  h(x,T)  it 
is  necessary  to  apply  the  Fourier  transform  relationships  given  by 
Equations  3-47  and  3-48.  In  the  discrete  case,  Equation  3-48  becomes 

h()iAx,sAT)  =  —  I  U(  aAx,pAv)  (3-66) 

p=0 

t=0,...N-l;  s=0,...M-l 
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Substitution  of  Equation  3-52  into  Equation  3-66  yields 


h()iAx,SAt)  =  |-  III  h(niAK,rAT)  aKatAv 

2ir  j„=o  r=0  p=0 


but  since 


£=0,...N-1;  s=0,...M-1 


(3-67) 


Y  e«>P<  *-'■)/«=  MS 
p=0 


sr 


(3-68) 


Equation  3-67  reduces  to 


N-1 

h()iAx,SAT)  =—  I  h(niAK,3AT)  aKatAv 

m=0 


t=0,...N-l;  s=0,...M-l 


(3-69) 


Now  use  Equation  3-52  to  specify  a  random  realization  of  the  Fourier 
transform  and  simplify  by  using  the  relationships  for  aK,  at,  and  av  to 
obtain 

N-1 


h(JlAx,sAT)  =  (BJl)^/^  I  [S(mAK,sAT)  r 

m=0 


i2TT«m/N 


ms 


(3-70) 


Jl=0,...N-l;  s=0,...M-l 


This  equation  is  the  final  statement  of  the  method  to  generate 
realizations  of  the  impulse  response  function  after  transionospheric  prop 
agation.  Note  the  interesting  property  observed  here  that  the  values  of 
the  impulse  response  function  are  uncorrelated  (or  delta-correlated)  in 
delay.  This  property  occurs  because  the  random  variable  rp,s  obeys  the 
correlation  relationship  specified  by  Equation  3-54. 


Again  it  is  quite  important  to  note  that  the  impulse  response 
function  generated  in  this  manner  is  a  complex  Gaussian  variate  whose 
amplitude  {h(x,T)h*(x,T)  follows  a  Rayleigh  probability  distribution 
function.  The  importance  of  this  fact  is  illustrated  in  the  following 
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section  where  two  examples  of  statistically  generated  signals  are  given 
and  compared  to  multiple  phase-screen  calculations. 

3.4  RAYLEIGH  COMPONENT  OF  THE  COMPOSITE  SIGNAL 

This  section  describes  the  steps  which  must  be  taken  to  generate 
statistical  realizations  of  the  impulse  response  function.  Several  exam¬ 
ples  are  discussed,  and  results  compared  to  MPS  results  presented  in 
Section  2  of  this  report. 

One  very  important  problem  remains  to  be  resolved  before  statis¬ 
tical  signal  realizations  can  be  generated.  The  problem  is  caused  by  the 
fact,  discussed  in  the  previous  section,  that  the  choice  of  Gaussian  ran¬ 
dom  variables  results  in  Rayleigh  amplitude  statistics  for  the  generated 
signal.  The  reason  for  the  problem  is  obvious  -  the  analytically  calcu¬ 
lated  generalized  power  spectrum  S{K,t)  includes  all  scattering  effects, 
but  the  limitation  imposed  by  requiring  that  the  generated  signal  have 
Rayleigh  statistics  effectively  requires  the  user  to  generate  only  the 
Rayleigh  portion  of  the  signal. 

The  solution  to  this  problem  is  pragmatic.  Arbitrarily  divide 
the  received  signal  into  two  disjoint  parts.  One  part,  called  the  disper¬ 
sive  component,  represents  the  effects  of  large  phase  structure  on  the 
signal.  The  other  signal  component,  called  the  Rayleigh  signal  component, 
is  that  part  of  the  received  signal  which  obeys  Rayleigh  statistics  and 
represents  severe  diffractive  or  angular  scattering  effects.  These  two 
signal  components  have  some  similarities  to  the  focus  and  scattering 
signal  components  upon  which  Fremouw's  (Fremouw,  et  al.,  1978)  two-com¬ 
ponent  signal  statistical  model  is  based. 

Now  for  the  case  of  a  phase  power  spectrim  which  represents  a 
range  of  scale  sizes,  the  scattering  processes  are  well  understood.  Large 
scale  sizes,  such  that  the  focal  length  of  the  irregularities  is  greater 
than  the  propagation  distance,  cause  large  phase  effects  but  cause  little 
amplitude  scintillation.  This  is  another  way  of  saying  that  larger  scale 
sizes  cause  dispersion  and  not  diffraction.  Small  scale  sizes,  whose 
focal  lengths  are  less  than  the  propagation  distance,  contribute  to  the 
scintillation  index  because  these  small  sizes  are  responsible  for  diffrac¬ 
tion  or  angular  scattering  effects. 
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Now  it  is  also  well  known  (Fante,  1975)  that  for  saturated  scin¬ 
tillation,  the  signal  amplitude  statistics  are  Rayleigh  and  the  scin¬ 
tillation  index  is  tnus  unity.  It  is  also  known  that  weak  scattering 
predictions  of  the  scintillation  index  in  either  the  Born  approximation  or 
the  Rytov  approximation  give  accurate  results  for  S4  £  0.5  (Crane,  1977). 

Now  assume  for  the  moment  that  saturated,  strong  scattering 
occurs  for  S4  >  0.5.  That  is,  with  the  heuristic  assumption  that  for 
$4  >  0.5,  the  amplitude  statistics  are  basically  Rayleigh  and  the  scintil¬ 
lation  index  is  effectively  unity,  the  S4  scintillation  index  can  be 
accurately  used  to  separate  the  Rayleigh  signal  component  from  the  disper¬ 
sive  component.  That  is,  since  the  scintillation  index  can  be  expressed 
as  the  integral  over  wavenumber  of  the  phase  PSD  times  a  Fresnel  filter 
factor,  the  integral  may  be  used  to  define  a  scale  size  which  separates 
dispersive  effects  from  diffractive  effects. 

In  this  section  the  discussion  is  limited  to  a  K"^  one-dimen¬ 
sional  phase  PSD.  The  integral  expression  for  S4  is  used  to  define  an 
outer  scale  size  for  the  Rayleigh  signal  component,  L[^,  which  is  in  turn 
used  to  define  a  new  phase  PSD  with  phase  standard  deviation  0^.  The 
Rayleigh  phase  standard  deviation,  and  the  Rayleigh  outer  scale 
size,  Lr,  are  then  used  to  define  the  generalized  power  spectrum  for 
purposes  of  statistical  signal  generation.  The  effect  of  this  ad  hoc 
procedure  is  to  generate  only  that  portion  of  the  signal  that  obeys 
Rayleigh  statistics.  The  rest  of  this  subsection  describes  the  details 
involved  in  obtaining  the  effective  PSD  for  the  Rayleigh  component  of  the 
received  signal , 

In  the  thin  phase-screen  approximation  (Salpeter,  1967)  the  S4 
scintillation  index  is  given  as  an  integral  over  wavenumber  by 

S4  =  4 /S  (Kj^)  sin2  (K2z*/2k)  dK^^  (3-71) 

••00  ^ 

(elongated  irregularities) 

^4  =  [(K^^2)2*/2k]  dK^dKy  (3-72) 

(isotropic  irregularities) 


no 


■fK'*\-^%’^ji:^jt-^-^  '-'  "<  ~- 


where 


L-L^L 


(3-73) 


The  geometry  is  shown  In  Figure  3-1  where  the  layer  thickness  L 
is  taken  as  zero  in  the  thin  phase-screen  limit.  The  transmitter  emits  a 
spherical  wave  which  propagates  a  distance  in  free-space  until  it 
penetrates  the  scattering  layer.  The  wave  then  propagates  a  distance  z^ 
to  the  receiver.  The  phase  power  spectra  are  given  as 


of  f((n-2)/2]L 

s  (K  )  = _ ♦_ _ _2 _ 

*  ^  h  r[(n-3)/2]  (1+k2lJ)("-2)/2 


(3-74) 


(elongated  irregularities) 


.J  r[(n-l)/2lL„ 

S  {K  ,K  )  - _ y _ __2. _ 

*  .  r[(n-3)/2)  [lt(K|«2)L2l("-»/2 


(3-75) 


(isotropic  irregularities) 


where  r  is  the  gamma  function.  It  is  seen  that  these  two  PSD  follow  the 
relationship  (Rufenach,  1975) 


W  -  f  *y 


(3-76) 


so  that  in  both  cases  the  one-dimensional  phase  PSD  has  the  desired  power 
law  form 


3.4.1  Rayleigh  Signal  Component  for  Elongated  Irregularities 


Now  in  the  case  of  elongated  irregularities,  the  substitution 


V  =  K2z*/2k 

A 


(3-77) 


may  be  applied,  and  the  integral  expression  Equation  3-71  may  be  written 
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2  _  8  ^2  o2  r[(n-2)/2]  7  Sin2(v2)dv 


S^  * _ 1. 

A  r[(n-3)/2]5i/2  J  (l+2V2/e)("-2)/2 


f  sit 
J  (1+5 


(3-78) 


where 


5  =  z*/kL2  =  XZ*/2itL2 
0  0 


(3-79) 


Now  define  the  value  V  such  that  the  scintillation  index  is  0.5 


a  /2  o2  r[(n-2)/2] 


(0.5)2  «  - - ^ 

✓it  r[(n-3)/2]  5^/2  J  (1+2VV5)(”’2)/2 


¥  ■ 

/ 


sin2(v2)dV 


(3-80) 


V  is  the  cutoff  value  chosen  so  that  wavenumbers  smaller  than  ✓!  kV'/z* 
contribute  a  value  of  0.5  to  the  scintillation  index  in  the  weak  scatter 
approximation.  Since  this  weak  scattering  expression  is  valid  for 
$4  £0.5  (Crane,  1977),  V‘  is  a  well  defined  quantity.  Now  define  the 
outer  scale  cutoff  value,  Lr,  such  that  scale  sizes  larger  than  or  equal 
to  Lr  contribute  a  value  of  0.5  to  the  S4  integration,  or 


U  =  (z*/2kV')^/2  =  (xz*/4itV')^/2 


(3-81) 


Equation  3-81  can  be  rewritten  as 

Lr/Lq  =  (xz*/4^L2V')^/2=  (5/2V')1/2 


(3-82) 


Hence  Lr  is  the  scale  size  which  separates  large  scale  sizes  from  the 
small  scale  sizes  that  cause  the  major  scintillation  effects.  Now  it  is 
apparent  that  the  value  of  Lr/Lq  is  dependent  upon  the  spectral  index 
of  the  one-dimensional  phase  fluctuations,  2-n,  the  outer  scale,  Lq,  the 
parameter  5  and  the  phase  standard  deviation,  In  the  case  of  inter¬ 
est  here  the  one-dimensional  phase  PSD  has  a  spectral  index  of  -3,  so  n  is 
taken  as  5. 


Figure  3-2  shows  Lr/Lq  as  a  function  of  with  5  treated 
parametrically.  Values  of  5  ranging  from  10~®  to  1  are  shown.  For  a  con¬ 
stant  value  of  ?,  there  is  a  minimum  value  of  the  phase  standard  deviation 
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below  which  no  solution  is  found  for  Values  of  a  below  the  minitnutn 

value  correspond  to  the  situation  where  all  irregularities  are  close  to 
the  phase-screen  and  cause  only  dispersive  effects.  Recall  that  the  focal 
length  for  an  irregularity  of  strength  (|)q  is 


F  =  TrL^/X(|> 

0  0 

(3-83) 

F/z*  =  itLVxz*,^„ 

0  0 

(3-84) 

So  as  0^  or  equivalently  (jiQ  increases  for  a  given  xz*/Lq^  or  5  value, 
the  focal  length  becomes  smaller  than  the  effective  propagation  distance. 
Thus  as  increases,  more  and  more  irregularities  contribute  to  angular 
scattering  and  thus  add  contributions  to  the  value  of  S^. 

As  0^  continues  to  increase  for  a  fixed  5  =  xz*/2itLQ^,  the 
focal  distance  for  more  and  more  of  the  ionization  irregularities  becomes 
smaller  than  the  effective  propagation  distance  z*.  Eventually,  for  0^ 
large  enough,  all  the  ionization  scale  sizes  contribute  to  the  S4  calcula¬ 
tion  so  that  the  received  signal  is  said  to  have  no  dispersive  signal  com¬ 
ponent.  In  that  case  all  the  original  phase  irregularities  contribute  to 
the  Rayleigh  signal  component  and  Lr  =  Lq. 


Now  that  a  scale  size  Lr  has  been  chosen  to  characterize  those 
sizes  which  contribute  to  the  Rayleigh  part  of  the  scattered  signal,  let 
us  choose  an  entire  phase  PSD  corresponding  to  the  Rayleigh  signal  compon¬ 
ent  with  outer  scale  size  Lr.  For  this  Rayleigh  PSD  to  match  the  origi¬ 
nal  PSD  at  small  scale  sizes,  as  it  should,  the  two  PSD's  must  have 
identical  limits  for  large.  In  the  limit  for  large 


(K,L„) 


0 

n^ 


2-n 

X 


(3-85) 


which  must  be  identical  to  the  Rayleigh  signal  component  PSD.  Now  for  the 
case  of  a  K"^  one-dimensional  phase  PSD,  n  =  5  in  Equation  3-74  and 
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■>5  <  \:<f- - ' 


-21 -2  -  -2  1-2 


Thus  the  choice 


%R  " 


(3-86) 


(3-87) 


yields  a  new  Rayleigh  signal  PSD  with  scale  size  dependence  identical  to 
the  original  phase  PSD  scale  size  dependence  at  small  scale  sizes  (large  K 
values). 

With  the  aid  of  Equation  3-87  it  is  possible  to  place  the  value 
of  the  phase  standard  deviation  of  the  Rayleigh  signal  component,  o^r,  on 
the  curves  shown  in  Figure  3-2.  Figure  3-2  can  also  be  cast  into  the  form 
of  Figure  3-3  which  shows  the  resulting  Rayleigh  phase  standard  deviation, 
O.R,  as  a  function  of  the  original  phase  standard  deviation,  o^.  For  a 
fixed  value  of  the  parameter  5  there  is  a  minimum  value  of  a.  below  which 
angular  scattering  is  negligible  and  does  not  reach  the  value  of  0.5. 
For  this  case  no  Rayleigh  signal  component  is  observed.  As  increases 
with  constant  5  the  value  of 

F/z*  =  iTLj/xz*a^ 

decreases  so  that  all  scale  sizes  eventually  contribute  to  the  Rayleigh 
signal  component  and  there  is  no  dispersive  component.  For  that  case  the 
Rayleigh  spectrum  is  identical  to  the  original  spectrim  with 

%R  ^  ^R  '  ^0’ 

This  case  corresponds  to  the  upper  portion  of  the  plot  where  would 
be  larger  than  if  any  curves  appeared  there. 

Returning  to  Figure  3-2  it  is  seen  that  for  a  constant  value  of 
5,  there  is  a  minimum  value  of  0^  below  which  does  not  reach  0.5. 

These  minimum  values  appear  as  vertical  lines  on  the  figure.  In  this  case 
of  weak  scattering  there  is  no  Rayleigh  signal  component.  For  intermediate 
values  of  0^  larger  than  the  minimum,  there  is  both  a  dispersive  signal 
component  and  a  Rayleigh  signal  component;  here  Lr  and  may  be  read 


114 


PHASE  STANDARD  DEVIATION  a,  (RADIANS) 


Figure  3-2.  Rayleigh  signal  component  outer  scale  versus  o  and 
0^1^,  elongated  Irregularities.  ^ 


Figure  3-3.  Rayleigh  signal  component  n  pversus  o,,  elongated 
irregularities. 


from  the  figure.  For  still  larger  values  on  any  constant  5  curve, 
angular  scattering  effects  dominate  and  the  Rayleigh  signal  may  be  gener¬ 
ated  directly  from  the  original  values  of  and  with  no  dispersive 
signal  component. 

As  discussed  in  Section  4.3,  an  important  signal  parameter  which 
characterizes  whether  dispersive  or  diffractive  effects  dominate  is 
defined  as 


(3-88) 


For  a  greater  than  unity  angular  scattering  or  diffraction  effects  domi¬ 
nate  the  signal  structure.  For  o  less  than  unity  dispersive  effects  pre¬ 
dominate.  Therefore,  it  is  necessary  that  the  value  of  a  for  the  Rayleigh 
signal  component  be  greater  than  unity  or  else  the  heuristic  basis  for  the 
choice  of  the  Rayleigh  PSD  loses  merit.  For  a  K"^  one-dimensional  phase 
PSD  the  value  of  a  is  given  in  Equation  4-106  as 


a 


^n(Lo/^i)X2*a^ 


(3-89) 


Now  for  the  Rayleigh  signal  component 


(3-90) 


Now  Equation  3-87  may  be  used  in  Equation  3-90  to  obtain 

“R  ’  {3-91) 

Since  results  already  exist  for  as  a  function  of  with  5  as  a 
parameter  it  is  easy  to  compute  or  as  a  function  of  0^  with  5  taken 
parametrically. 
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Figure  3-4  shows  the  value  of  oj^,  for  the  Rayleigh  signal  compo¬ 
nent,  as  a  function  of  for  the  composite  signal  for  values  of  the 
parameter  c  ranging  from  10“®  to  100.  The  ratio  of  outer  scale  to  inner 
scale  Lr/ar  is  taken  as  100  in  the  figure;  this  value  is  a  reasonable 
minimum.  Values  of  the  phase  standard  deviation  for  the  Rayleigh  signal 
component,  oxn,  are  also  shown  parametrically  ranging  from  1  to  1000  radi¬ 
ans.  As  can  be  seen  in  the  figure,  the  minimum  value  of  is  2.5  for 
0^1^  greater  than  1.0  radian.  Since  values  of  less  than  1.0  radian  are 
prohibited  (see  the  discussion  of  Figure  3-3),  the  value  of  2.5  is  a  true 
minimum  for 

The  fact  that  or  for  the  Rayleigh  signal  component  is  always 
greater  than  unity  is  important  for  two  reasons.  The  first  reason,  already 
discussed,  is  that  the  Rayleigh  signal  component  is  thus  characterized  by 
angular  scattering  or  diffractive  effects.  The  second  reason  is  that  the 
time  delay  jitter  for  the  Rayleigh  signal  component  is  then  independent  of 
o  and  is  simply  related  to  the  coherence  bandwidth,  o^qi,,  defined  in 
Section  4. I. 5. 2. 


From  Equation  4-135  the  time  delay  jitter  for  the  composite 
signal  is  given  by 


M2  oV 


(3-92) 


Now  it  will  be  shown  that  the  time  delay  jitter  is  always  inversely  pro¬ 
portional  to  the  coherence  bandwidth  if  the  signal  has  a  Rayleigh  com¬ 
ponent.  Furthermore,  if  the  signal  has  a  Rayleigh  component,  the  value  of 
the  coherence  bandwidth  is  the  same  as  it  is  for  the  composite  signal. 

Now  if  a  is  large,  diffraction  is  the  dominant  effect  and  the 
time  delay  jitter  becomes 


•C> 

r 


If  a  is  small,  dispersive  effects  dominate  arid  the  Rayleigh 
signal  component  must  be  generated  to  isolate  the  diffractive  processes. 
If  the  Rayleigh  signal  component  exists,  then  or  is  greater  than  unity 
and 


1 


'tR  -  — 


(3-94) 


/2  (I), 


cohR 


where  the  additional  R  subscripts  are  used  to  denote  the  Rayleigh  signal 
component. 


Now  for  a  K"^  one-dimensional  phase  PSD,  the  value  of  o^oh 
given  by  Equation  4-74 


"coh 


2iru 
0  0 


iin(L„/t,)Xz*o2 


(3-95) 


where  Ag  and  A2  are  taken  from  Equations  4-44  and  4-45.  Equation  3-95 
is  valid  for  the  composite  signal.  Now  since,  from  Equation  3-87 


'■r/'+r  “  '•o''”? 


(3-96) 


the  value  of  the  coherence  bandwidth  for  the  Rayleigh  signal  component  is 
approximately  identical  to  the  value  of  the  coherence  bandwidth  for  the 
composite  signal. 


‘’cohR  “  “coh 


since  the  only  difference  is  the  logarithmic  term  which  will  have  little 
effect. 


At  this  point  it  is  useful  to  consider  the  important  signal 
parameter  Xg  defined  in  Equation  4-73  as  the  signal  decorrelation  dis¬ 
tance.  For  a  K"^  one-dimensional  phase  PSD  the  values  of  Aq  and  A2  are 
given  by  Equations  4-44  and  4-45  and  the  value  of  the  signal  decorrelation 
distance  is  given  by 
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The  value  of  the  signal  decorrelation  distance  for  the  Rayleigh  signal 
component  may  be  written  as 


"oR 


‘"(''■p/»i>“jR"t 


(3-98) 


Equation  3-96  may  be  substituted  into  Equation  3-98  with  the  result  that 


A  _  a  t 

oR  0 


(3-99) 


Exact  equality  is  not  obtained  because  of  the  differences  in  the  logarithm 
term.  This  important  result  states  that,  if  a  Rayleigh  signal  component 
exists,  its  signal  decorrelation  distance  is  approximately  equal  to  the 
signal  decorrelation  distance  of  the  original  composite  signal. 


3.4.2  Rayleigh  Signal  Component  for  Isotropic  Irregularities 

In  the  case  of  isotropic  irregularities  the  transformations 

(3-100) 


and 


V  =  Kj^Z*/2k 


(3-101) 


may  be  introduced  into  Equation  3-72  and  the  identical  definition  of  V 
established  as  was  used  in  the  case  of  elongated  irregularities  so  that 
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defines  V.  The  calculation  of  the  Rayleigh  signal  component  parameters 
0(|)R’  ^  “R  accomplished  with  the  same  procedures  used  in 

the  elongated  irregularity  case. 

Figures  3-5,  3-6,  and  3-7  show  the  results  for  Lp/Lg,  and  op 
as  functions  of  the  value  of  the  phase  standard  deviation  of  the  com¬ 
posite  signal  and  various  values  of  the  parameter  xZ*/2TrLo^.  As  can  be 
seen  the  general  behavior  of  these  results  is  similar  to  the  results  for 
the  case  of  elongated  irregularities.’ 

3.5  EXAHPLES  OF  STATISTICALLY  GENERATED  SIGNALS 

3.5.1  Example  1:  Diffraction  Important 

In  this  section  examples  of  two  statistically  generated  signals 
are  presented  and  the  results  compared  to  the  MRS  results  presented  in 
Section  2.  These  two  examples  present  interesting  tests  of  the  signal 
generation  algorithms  since  one  case  is  dominated  by  diffractive  effects 
and  the  other  is  dominated  by  dispersive  effects. 

First  consider  the  MRS  calculation  described  in  detail  in 
Section  2.6.  In  this  example  the  carrier  frequency  is  2.25  GHz,  the  phase 
standard  deviation  is  350  radians,  the  outer  scale  is  1  km  and  the  inner 
scale  is  10  m.  The  incident  wave  is  plane  (infinite  z^)  and  the  scaUer-  - 
ing  layer  is  1000  km  thick  with  an  additional  free-space  propagation 
distance  of  1000  km  from  the  exit  of  the  layer  to  the  receiver  plane. 

This  MRS  calculation  is  applicable  only  to  a  two-dimensional  geometry  and 
thus  corresponds  to  the  theoretical  development  for  the  infinitely  elon¬ 
gated  irregularity  case. 

Now  since  the  computed  values  of  the  generalized  power  spectrum 
S(K,t)  have  been  specialized  to  the  thin  phase-screen  case  for  conveni¬ 
ence,  it  is  necessary  to  find  the  thin  phase-screen  geometry  equivalent  to 
the  MRS  geometry  just  described.  For  the  case  at  hand  the  single  thin 
phase-screen  is  centered  in  the  scattering  layer  and  is  located  a  distance 
of  1500  km  from  the  receiver  plane.  The  value  of  xzj,/2tiLo^  is  then 
0.031831.  Now  using  the  results  shown  in  Figures  3r2  and  3-3  it  is  seen 
that  for  the  Rayleigh  signal  component 


K-^  PHASE  PSO 
ISOTROPIC  IRREGULARITIES 
£  >  XZ*/2iL? 


o.o  =  194  rad 


Lp  =  564  m 

Now  in  order  to  obtain  results  for  comparison  to  the  MPS  results 
presented  in  Section  2.6,  it  is  convenient  to  numerically  generate  the 
solution  to  the  parabolic  wave  equation  U(x,Zy.,v4-a^)  by  taking  the 
Fourier  transform  of  the  impulse  response  function 


U(x,z^,v+o3q)  =  /  h(x,z^,T)e"‘‘'’''dT  (3-103) 

Equivalently  U  may  be  directly  obtained  from  Equation  3-51.  Equation  3-103 
is  identical  to  Equation  3-47  where  the  z^  and  Uq  dependence  are  included 
for  completeness.  Note  that  the  impulse  response  function  h  is  also 
dependent  upon  the  carrier  frequency.  Realizations  of  the  impulse  response 
are  obtained  in  discrete  form  using  Equation  3-70.  The  generalized  power 
spectrum  S(K,t)  is  given  for  the  case  of  elongated  irregularities  by 


S(K.,)  =  exp  j -  1  -  1  an  i  (3-m) 


where  the  parameters  cxr,  and  are  given  by  Equations  4-106, 

4-73,  and  4-74,  respectively  and  of  course  refer  to  the  Rayleigh  component 
of  the  received  signal  which  is  the  only  component  one  can  generate  with 
this  technique.  Ag  and  A2  which  are  required  in  Equations  4-73  and  4-74 
may  be  obtained  from  Equations  4-44  and  4-45.  The  carrier  frequency 
dependence  of  the  impulse  response  function  is  evident  in  the  expressions 
for  the  parameters  o,  (*toh»  *o* 


The  Fourier  transform  given  by  Equation  3-103  may  be  implemented 
discretely  by  using  Equations  3-51  and  3-52  from  the  earlier  development. 
In  order  to  directly  compare  to  the  MPS  case,  32  discrete  frequencies  are 
used  over  a  bandwidth  of  4.096  MHz.  The  same  number  (32 j  of  delay  cells 
are  used  with  delay  spacing 
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Figure  3-8.  Nonnalized  power  impulse  response  function,  a  =  100. 

At  =  244.1  X  10"^  sec 
and  thus  a  total  time-delay  interval 
T  =  32  At  =  7.812  x  10'®  sec 

Values  of  the  important  Rayleigh  signal  parameters  wtohR  and  or  are  given 
in  the  plane  wave,  thin  phase-screen  case  as 

“cohR  "’'“oV^^22r®5R  rad/sec 

“R  "  V®^R“cohR  =  91.7 

To  perform  the  integration  over  t  which  is  required  in  Equation  3-51  or 
Equation  3-103,  the  variable  t  ranged  discretely  from  -2.44x10"^  to 
7.32x10'®  sec  so  that  values  of  a)^,Qj,pT  ranged  from  -0.19  to  5.82.  As  can 
be  seen  from  Figure  3-8  which  shows  the  normalized  power  impulse  response 
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function  defined  in  Section  4.3,  the  range  of  non-zero  values  of  the  power 
impulse  response  is  about 

0  1  “coh  ^  i  3 

for  a  =  100.  In  the  case  at  hand  a  is  close  to  100  and  so  the  range  of  _ 
values  taken  in  the  numerical  calculation  includes  all  important  t  values; 

Another  way  to  check  the  accuracy  of  the  numerical  calculation 
is  to  form  the  integral 

N-1  M-1 

5]  S(mAK,sAT)  aKAt 

m=0  s=0 

This  integral  is  the  discrete  equivalent  of 

eo  00 

J  JS(K,T)dKdT 

•  00  mqo 

and  should  be  equal  to  unity  since  the  continuous  integral  is  equal  to 
r(Ax=0,Af=0)  according  to  Equation  4-84  or  4-87.  This  important  quantity 
is  always  formed  and  checked  in  every  statistical  signal  generation. 

For  this  example  there  are  1024  cells  in  the  x-direction,  with  a 
total  length  of  208.33  m  and  a  spacing  of  0.20345  m  which  is  precisely 
three  times  the  grid  spacing  used  in  the  MRS  calculation. 

In  order  to  obtain  results  for  statistical  averaging  purposes, 
the  statistical  signal  generation  algorithm  is  performed  five  times  with 
five  different  random  nunber  sequences.  As  described  above,  for  each 
individual  signal  generation  calculation,  32  frequencies  are  generated 
with  1024  x-points  for  each  frequency.  Ihis  signal  is  referred  to  as  a 
single  statistical  realization  and  is  similar  to  a  MRS  frequency  selective 
signal  realization. 

Figures  3-9  to  3-14  show  a  comparison  of  the  mutual  coherence 
function  r(Ax,Af)  for  a  single  statistical  realization  versus  the  strong 
scattering  theory  described  in  Section  4.  Figure  3-9  shows  the  amplitude 
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of  r  corresponding  to  values  of  ax  of  0.0,  0.61  m,  1.22  m,  1.83  m  and  2.44 
m  versus  theory.  These  values  of  ax  are  chosen  to  correspond  to  the  MPS 
values  presented  in  Figures  2-26  to  2-31.  The  theoretical  results  pre¬ 
sented  here  are  similar  to  those  presented  in  Figures  2-26  to  2-31  except 
that  here  the  thin  phase-screen  theoretical  limit  is  taken,  it  is  seen 
that  the  difference  in  the  phase  of  r  between  the  thick  layer  theory  and 
the  thin  phase-screen  theory  is  simply  a  term  that  is  linearly  proportion¬ 
al  to  the  frequency  deviation  from  the  carrier  frequency.  In  Section  4 
the  received  time-domain  signal  is  shown  to  be  related  to  the  Fourier 
transform  of  r  so  that  this  phase  term  that  is  proportional  to  frequency 
yields  only  an  additional  time  delay.  As  can  be  seen  in  Figure  3-9  the 
agreement  between  statistical  results  and  theory  is  good  but  is  somewhat 
worse  than  the  agreement  between  the  MPS  results  and  theory. 

Figures  3-10  to  3-14  compare  the  phase  of  r  computed  from  the 
single  statistical  realization  versus  the  thin  phase-screen  theory. 
Agreement  is  very  good  in  the  most  important  region  close  to  zero  Af  where 
the  power  in  the  spectrum  of  the  modulation  waveform  is  concentrated,  but 
elsewhere  the  phase  of  r  diverges  somewhat  from  the  analytic  results. 

Figures  3-15  to  3-20  compare  the  value  of  r  obtained  by  averag¬ 
ing  the  results  of  five  independent  statistical  realizations  to  the  ana¬ 
lytic,  thin  phase-screen  results.  It  is  apparent  in  Figure  3-15  that  the 
average  amplitude  of  r  compares  quite  well  to  the  theory  and,  in  fact, 
does  better  than  the  MPS  calculation  at  matching  the  theoretical  thin 
phase-screen  results.  The  agreement  with  theory  is  quite  good  for  the 
entire  range  of  Af  for  all  five  of  the  ax  values  shown. 

Figures  3-16  to  3-20  show  a  comparison  of  the  phase  of  r  for  the 
average  of  five  statistical  realizations  versus  the  analytic  thin  phase- 
screen  value.  As  expected  from  the  averaging  process,  the  variations  in 
the  phase  are  reduced  in  comparison  to  the  variations  observed  for  the 
case  of  a  single  realization.  In  all  five  figures  the  average  statistical 
phase  and  the  analytic  phase  agree  quite  well  near  zero  Af  but  disaqree 
somewhat  for  |Afl  greater  than  about  1  MHz.  This  minor  disagreement  .i  the 
phase  of  r  does  not  appear  to  be  important  to  the  resulting  time-domain 
signals  to  be  presented  next.  The  most  important  aspect  of  these  compari¬ 
sons  is  the  good  agreement  between  numerical  results  and  theory  near  zero 
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Figure  3-15.  Comparison  of  amplitude  of  r,  average  of  five  statisti 
cal  realizations  versus  theory. 


Figure  3-16.  Comparison  of  phase  of  r,  average  of  five  statistical 
realizations  versus  theory,  ax  =  0. 
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Figure  3-20.  Comparison  of  phase  of  r»  average  of  five  statistical 
realizations  versus  theory,  ax  -  2.44  m. 


Af.  This  region  is  the  most  critical  since  zero  Af  corresponds  to  the 
carrier  frequency  where  most  of  the  energy  in  the  modulation  waveform  is 
concentrated. 

Figure  3-21  shows  an  example  of  the  envelope  of  the  received 
triangular  waveform  as  determined  by  the  statistical  signal  generation 
technique.  The  received  complex  envelope  is  given  by 

e{x,z^,T)=L  /M(v)U(x.z^.vfa>o)e'''''dv  (3-105) 

fcIT  mCO 


where  v=2Trf  and  T^.  is  488.3  nsec,  twice  the  duration  of  the  transmitted 
triangular  pulse.  U(x,z^,v+(Uo)  is  determined  from  the  discrete  form 
of  Equation  3-103  where  the  impulse  response  function  h(x,z^,T)  is  the 
statistically  generated  impulse  response  function  obtained  using 
Equation  3-70.  In  Figure  3-21  samples  of  the  amplitude  |e(x,Zr,T)l  are 
shown  for  101  values  of  x  separated  in  distance  by  0.61  m,  the  MRS  grid 
interval.  Thus  Figure  3-21  is  directly  comparable  to  Figure  2-32  in 
Section  2.  The  time  domain  results  shown  in  these  two  figures  are  both 
independent  realizations  of  a  stochastic  process  and  hence  cannot  be 
identical.  However,  the  similarity  in  overall  structure  indicates  that 
the  s'^atistically  generated  realization  is  an  adequate  approximation  to 
the  exact  MRS  calculation. 


Figure  3-22  shows  the  average  pulse  shape  of  the  received  time 
domain  signal  which  is  obtained  by  averaging  all  the  5120  realizations  of 
|e(x,Zr,T)|  generated  in  five  runs  with  1024  x-samples  per  run.  This 
average  pulse  shape  is  found  to  be  quite  similar  to  the  result  for  the  MRS 
calculation  shown  in  Figure  2-33.  The  overall  appearance  of  the  two  pulse 
shapes  is  very  similar  with  a  noticeable  concentration  of  power  at  small 


133 


TIME  DEUY  (CHIPS) 


Figure  3-21.  Envelope  of  the  received  time-domain  waveform- 
statistical  signal  realization. 
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Figure  3-22.  Average  received  time-domain  envelope— statistical 
signal  realization. 


delays  for  the  MPS  case  as  compared  to  a  slightly  more  spread  out  situa¬ 
tion  for  the  statistical  signal  generator  case.  For  the  MPS  case  shown  in 
Figure  2-33,  the  mean  time  delay,  <t>,  is  1.1  chips  and  the  time  delay 
jitter,  0^,  is  1.4  chips.  For  the  statistical  signal  shown  in  Figure 
3-22  the  mean  time  delay  is  1.2  chips  and  the  time  delay  jitter  is  1.7 
chips.  In  both  cases  a  chip  is  the  unit  of  time  here  equal  to  488.3 
nsec.  The  theoretical  result  in  the  thin  phase-screen  case  is  1.3  chips 
for  the  mean  time  delay  and  1.8  chips  for  the  time  delay  jitter. 

Thus  rather  good  agreement  is  obtained  for  the  MPS  results  and 
the  statistical  signal  results  when  compared  to  theory.  One  should  keep 
in  mind  that  the  MPS  calculation  is  a  direct  solution  to  the  parabolic 
wave  equation  and  involves  no  approximations.  The  theoretical  calculation 
requires  a  number  of  approx imations^to  obtain  solutions  for  average  quan¬ 
tities  and  hence  suffers  in  accuracy  relative  to  the  very  general  MPS 
calculation.  The  theoretical  calculation  includes  both  dispersion  and 
diffraction  effects  but  the  statistical  signal  generation  technique  is 
only  applicable  to  the  diffraction  effects.  The  rather  good  agreement 
shown  in  this  example  indicates  that  the  heuristic  method  devised  to 
calculate  the  Rayleigh  or  diffraction  part  of  the  composite  signal  is 
adequate  to  the  task. 

3.5.2  Example  2:  Focusing  Important 

As  another  example  of  statistical  signal  generation  consider  a 
case  where  large  scale  structure  leads  to  focusing  and  to  other  effects  of 
large  phase  irregularities.  As  in  the  previous  example  the  MPS  results 
are  compared  to  the  statistically  generated  results. 

The  scattering  geometry  is  identical  to  that  of  the  previous 
example.  A  plane  wave  propagates  through  a  1000  km  thick  scattering  layer 
and  then  propagates  in  free-space  an  additional  1000  km  to  the  receiver 
plane.  The  carrier  frequency  here  is  7.5  GHz,  the  phase  standard  deviation 
is  300  radians,  the  outer  scale  is  10  km  and  the  inner  scale  is  10  m.  In 
the  MPS  calculation,  ten  phase-screens  are  uniformly  spaced  over  1000  km 
to  represent  the  thick  scattering  medium.  A  grid  of  length  100  km  is 
divided  into  16384  points  to  represent  the  signal  amplitude  and  phase. 
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For  this  example,  32  discrete  frequencies  are  used  over  a  two-sided  band¬ 
width  of  80  MHz.  This  bandwidth  is  chosen  as  the  null-to-null  bandwidth 
of  a  PN  spread  spectrum  code  with  chip  rate  of  40x10®  chips/sec  or  to 
represent  a  single,  band-limited,  triangular,  modulated  waveform  with  T(. 
equal  to  25x10"®  sec. 

Figures  3-23  to  3-30  show  the  resulting  time-domain  waveform  and 
some  related  measurements  from  the  MPS  realization.  Figure  3-23  shows  101 
samples  of  the  magnitude  of  the  envelope  of  the  received  time-domain  wave¬ 
form.  Each  of  these  101  samples  is  taken  from  equally  spaced  points  along 
the  entire  MPS  calculation  grid.  Thus  each  curve  is  separated  in  distance 
by  990  m  and  very  little  resolution  along  the  x-axis  is  shown.  More  reso¬ 
lution  is  shown  in  Figure  3-24  which  shows  the  received  time-domain  en¬ 
velope  for  the  first  101  MPS  grid  points.  The  spacing  between  each  grid 
point  is  6.1035  m  corresponding  to  16384  cells  equally  spaced  over  a  100 
km  distance  so  that  the  101  curves  shown  in  the  figure  correspond  to  a 
total  distance  of  616  m  along  the  MPS  grid. 

The  effect  of  the  large  scale  phase  irregularities  is  apparent 
in  Figure  3-23  and  3-24  in  the  appearance  of  many  noticeable  instances  of 
signal  focusing  where  the  peak  amplitude  is  greater  than  unity.  This 
focusing  results  in  a  measured  value  of  the  S4  scintillation  index  of 
1.15.  Values  of  S4  greater  than  unity  are  indicative  of  focusing  ef¬ 
fects.  The  focal  length  for  the  largest  (10  km)  irregularities  is  given  by 

'^0 


which  is  26,180  km.  Thus  the  receiver  plane  is  well  inside  the  focal 
distance  of  the  largest  irregularities  and  this  example  is  dominated  by 
large  phase  effects.  These  large  phase  effects  are  apparent  in  Figure 
3-23  where  the  received  signal  is  seen  to  vary  in  mean  time  delay  with  no 
dispersive  signal  spreading.  Oispersive  effects  are  not  strong  enough 
here  to  cause  pulse  spreading. 

Figures  3-25  and  3-26  further  illustrate  the  point.  These 
figures  show  the  mean  time  delay  <t>  and  the  time  delay  jitter  a-  as  a 
function  of  distance  along  the  MPS  grid.  As  can  be  seen  the  mean  time 
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I  3-23.  Envelope  of  the  received  time-domain  waveform— MRS 
signal  realization. 
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Figure  3-24.  Envelope  of  the  received  time-domain  waveform— MPS 
signal  realization. 
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Figure  3-25.  Mean  tl 
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delay  varies  slowly  with  distance  and  corresponds  to  the  location  of  the 
signal  peaks  shown  in  Figure  3-23.  Notice  the  negative  mean  time  delay 
(early  arrival)  at  the  beginning  of  the  realization  and  the  increasingly 
late  time  of  arrival  (positive  delay)  around  40  km  in  the  realization. 
Both  of  these  features,  are  evident  in  Figure  3-23. 

Figure  3-26  shows  the  value  of  obtained  for  each  individual 
MRS  cell  in  the  MRS  grid  as  a  function  of  distance  along  the  MRS  grid. 
For  the  most  part  the  value  of  0,^  is  constant  and  equal  to  the  spread  of 
a  triangular  signal,  0.316  T^.  (see  Equation  2-73).  Thus  the  major  dis¬ 
persive  effect  in  this  MRS  calculation  is  increased  time  delay  that  is 
accompanied  by  very  little  pulse  broadening. 

However  in  the  range  from  40  to  50  km  along  the  MRS  grid  there 
is  some  very  interesting  phenomena  occurring  with  evidence  of  isolated, 
high  measurements  of  <t>  and  o^.  Figure  3-27  shows  the  peak  amplitude 
of  the  received  time-domain  envelope  as  a  function  of  distance  along  the 
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Figure  3-27.  Peak  signal  amplitude  as  a  function  of  distance  along 
the  MRS  grid. 
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MPS  calculation  grid.  Notice  that  there  is  an  extended  region  of  about  10 
dB  fading  at  around  40  km  and  a  region  of  very  deep  fading  or  defocusing 
near  45  km.  These  deep  fading  regions  are  responsible  for  the  large 
observed  values  of  <t>  and  o^. 

Figures  3-28  and  3-29  show  closeups  of  the  plots  of  <t>  and 
for  the  region  from  40  km  to  50  km  on  the  MPS  grid.  It  can  be  seen  that 
the  large  values  of  <t>  and  near  45  km  coincide  with  the  region  of  the 
deepest  signal  fades  shown  in  Figure  3-27.  It  is  also  apparent  that  the 
isolated  large  <t>  and  measurements  coincide  with  small  regions  where 
the  signal  level  is  very  low.  These  isolated,  large  <t>  and  admeasure¬ 
ments  are  not  artifacts  or  errors  in  the  simulation  but  rather  they  are 
caused  by  large  phase  decorrelation  across  the  signal  spectrum  which  has  a 
much  higher  probability  of  occurrence  during  deep  fading  conditions.  The 
reason  for  the  higher  than  normal  phase  decorrelation  during  deep  fades  is 
that,  in  general,  fades  are  caused  by  cancellation  of  signals  coming  from 
different  directions.  For  very  deep  fades  to  occur,  such  cancellation 
must  occur  over  much  of  the  signal  bandwidth.  Under  these  circumstances 
very  little  change  in  frequency  (or  direction  of  arrival)  is  required  to 
obtain  quite  a  large  change  in  signal  phase. 

Figure  3-30  shows  the  received  time-domain  signal  for  the  range 
from  44.5  km  to  46.1  km  on  the  MPS  grid.  Notice  the  increased  mean  time 
delay  of  about  0.4  chips  at  the  front  of  the  figure  caused  by  propagation 
through  a  large  striation.  At  the  back  of  the  figure,  corresponding  to 
the  very  deep  fade  location,  phase  decorrelation  across  the  signal  spec- 
trim  causes  a  very  extended  received  time-domain  waveform  so  that  some 
aliasing  is  apparent.  This  type  of  deep  fading  behavior  is  associated 
with  the  occurrence  of  focusing  and  defocusing  and  is  not  modeled  by  the 
statistical  signal  generation  of  Rayleigh  fading  signals.  However  this 
behavior  is  handled  correctly  in  the  MPS  propagation  simulation. 

Figure  2-37  in  Section  2  shows  a  comparison  of  the  magnitude  of 
the  two-frequency  mutual  coherence  function  calculated  from  the  MPS  reali¬ 
zation  versus  the  theoretical  strong  scatter  solution  given  in  Section  4. 
As  seen  in  the  figure,  there  is  some  disparity  between  the  two  calcula¬ 
tions.  The  source  of  the  differences  lies  in  the  fact  that  the  theoretical 
calculation  is  a  strong  scatter  approximation  not  exactly  valid  for  the 
focus ing/defocusing  conditions  characteristic  of  this  MPS  calculation. 
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Figure  3-30.  Received  envelope  for  MRS  calculation  from  44.4  to 
45.0  km  along  the  MPS  grid. 

In  order  to  generate  the  Rayleigh  portion  of  this  time-domain 
signal  it  is  first  necessary  to  determine  the  values  of  the  phase  standard 
deviation  <j^r  and  the  outer  scale  Lr  of  the  Rayleigh  signal  compon¬ 
ent.  These  values  may  be  obtained  from  Figures  3-2  and  3-3  for  the  value 
of  C  of  9.5x10"®  as 

®(|)R  ■ 

Lfj  =  168  in 

Values  of  the  important  signal  parameters  or  and  ucohR  found  as 
“cohR  =  =  5.97x10®  rad/sec 

“R  "  “o°<t)R''“cohR  " 
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where  the  Inner  scale  for  the  Rayleigh  signal  component  is  taken  as  1.68 
m.  Since  the  inner  scale  always  appears  in  the  logarithm  for  a  K”^  one¬ 
dimensional  phase  power  spectrum,  small  changes  in  its  value  have  little 
effect. 

As  is  the  case  in  the  previous  example,  it  is  desired  to  obtain 
a  statistical  realization  of  the  received  complex  envelope  for  canparison 
to  the  MRS  results.  Thus  realizations  of  the  function  U(x,Zp,v+ato) 
are  required. 

Again  U(x,Zj.,v+a)^j)  is  generated  at  32  discrete  frequencies 
over  the  null-to-null  80  MHz  bandwidth  which  corresponds  to  a  chip  dura¬ 
tion,  T(.,  of  2.5x10"®  sec  and  a  PN  code  rate  of  40  ftps.  Thus  32  delay 
cells  are  used  with 

At  =  1.25  X  10"®  sec 
and  a  total  delay  interval  of 

T  =  32  At  =  4.0  x  10"^  sec 

To  generate  the  impulse  response  function  h(x,z^,T)  the  dis¬ 
crete  values  of  the  delay  t  range  from 

-2.5  X  10"®  sec  £  t  £  3.625  x  10"^  sec 
which  corresponds  to 

-14,93  <  a)^^^T<  216.5 

Figure  4-3  for  ot  equal  to  10  shows  that  all  the  power  in  the  impulse 
response  function  is  contained  within  the  limits 

-0-5  <  <  4 

SO  that  the  range  cited  above  is  more  than  sufficient  to  obtain  all  the 
power.  In  fact  the  samples  of  S(K,t)  are  negligible  outside  the  limits 
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Thus,  in  this  example,  the  impulse  response  function  h(x,Zp,T)  has  pro¬ 
perties  similar  to  a  delta  function  since  it  is  non-zero  at  only  a  few  ’ 
values  of  delay.  The  transform  of  the  impulse  response  function 
U(x,Zj.,v+ajQ)  is  then  relatively  constant  with  frequency  v  and  thus 
exhibits  high  correlation  properties  with  respect  to  v. 

As  in  the  previous  example,  1024  cells  are  used  to  represent  the 
variation  in  x  with  a  value  of  ax  of  2.03  m  which  is  precisely  3  times  the 
value  used  in  the  MPS  calculation. 

Figure  3-31  shows  the  amplitude  of  the  received  complex  envelope 
|e{x,Zr,T)|  as  a  function  of  x  and  t.  101  curves  are  shown  each  repre¬ 
senting  the  received  signal  as  a  function  of  t  and  separated  in  distance  x 
by  6.10  m.  This  figure  should  be  compared  to  Figure  3-24  which  shows 
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Figure  3-31.  Received  envelope  for  statistical  signal  realization. 
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similar  results  for  the  MPS  calculation.  It  is  apparent  that  this  reali¬ 
zation  of  the  Rayleigh  component  of  the  received  signal  has  none  of  the 
focusing/defocusing  properties  or  the  large  slow  variations  in  mean  time 
delay  evident  in  the  MPS  calculation.  This  behavior  is,  of  course,  a 
consequence  of  the  definition  of  the  Rayleigh  signal  component  which  is 
that  portion  of  the  signal  with  Rayleigh  fading  statistics. 

Comparison  of  Figures  3-24  and  3-31  show  that  the  statistical 
calculation  for  the  individual  time-domain  signals  are  similar  to  the  MPS 
calculation.  However  the  large  scale  focusing  and  defocusing  effects  are 
not  represented.  In  addition,  the  large  slow  changes  in  the  mean  arrival 
time  caused  by  large  scale  phase  variations  are  also  not  included  in  the 
statistical  realization. 

Numerical  calculation  of  the  two-frequency  mutual  coherence 
function  for  this  statistical  signal  realization  shows  constant  correla¬ 
tion  with  frequency  for  both  amplitude  and.  phase  of  r  for  all  values  of 
AX.  This  correlation  property  has  already  been  discussed  and  is  another 
indication  of  the  absence  of  severe  frequency  selective  scattering  and 
diffractive  effects  in  this  example. 

It  should  be  noted  that  the  comparisons  shown  in  this  example 
emphasize  the  differences  in  the  signal  realizations  formed  by  the  MPS 
simulation  and  the  statistical  signal  generation  technique.  This  approach 
was  taken  because  it  yields  an  interesting  discussion  of  the  different 
aspects  of  electromagnetic  wave  propagation.  For  the  study  of  the  effects 
of  strong  scintillation  on  receiver  performance,  it  is  always  sufficient 
to  consider  only  Rayleigh  fading,  which  is  the  worst  case  for  a  one-way 
propagation  path  through  strongly  turbulent  ionized  media. 


SECTION  4 

TWO-FREQUENCY  MUTUAL  COHERENCE  FUNCTION 


In  this  section  an  analytic  solution  is  obtained  for  the  two- 
position,  two-frequency  mutual  coherence  function  for  spherical  wave  prop¬ 
agation.  It  is  assumed  that  strong  scattering  conditions  prevail  and  that 
the  quadratic  approximation  to  the  phase  structure-function  is  therefore 
valid. 

This  approximation  was  used  by  Sreenivasiah,  et  al.,  (1976)  for 
the  case  of  plane  wave  propagation.  Furthermore,  Sreenivasiah  and 
Ishimaru  (1979)  have  recently  generalized  their  previous  results  to  obtain 
the  two-frequency  mutual  coherence  function  for  beam  wave  propagation  in 
homogeneous  turbulence  using  the  quadratic  structure-function  approxima¬ 
tion.  More  recently  the  two-position,  two-frequency  mutual  coherence 
function  was  obtained  for  spherical  wave  propagation  using  the  extended 
Huygens-Fresnel  principle  (Fante,  1981).  Although  the  quadratic  struc¬ 
ture-function  approximation  can  sometimes  lead  to  difficulties  (Wandzura, 
1980)  it  is  appropriate  for  calculation  of  the  two-frequency  mutual  coher¬ 
ence  function  but  not  for  calculation  of  higher  moments  of  the  field 
(Fante,  1980).  Fante  discusses  the  accuracy  of  the  quadratic  structure- 
function  approximation  (Fante,  1981)  for  the  case  of  atmospheric  turbu¬ 
lence  with  a  Kolmogorov  power  spectriin  of  irregularities.  He  has  found 
that  the  accuracy  is  a  function  of  the  irregularity  power  spectrum  and  of 
the  strength  of  the  turbulence  (Private  Communication,  1982),  with 
accuracy  increasing  for  stronger  scattering. 

The  strong  scatter  solution  presented  in  this  section  forms  the 
basis  of  the  statistical  signal  generation  technique  described  in  Section 
3  of  this  report.  Here  the  turbulent  ionized  medium  occupies  a  thick 
finite  layer  with  transmitter  and  receiver  located  in  free-space  on  oppo¬ 
site  sides  of  the  layer. 
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General  analytic  solutions  are  derived  for  two  cases.  In  the 
first  case  the  random  slab  is  represented  by  a  one-dimensional  power  spec¬ 
trum  of  electron-density  fluctuations  which  corresponds  to  propagation 
through  elongated  irregularities  as  would  occur  for  an  equatorial  satel¬ 
lite  link  to  a  ground  station.  In  the  second  case  the  random  slab  is 
represented  by  isotropic  ionization  irregularities  which  would  occur  in 
the  ionosphere  for  propagation  roughly  along  the  direction  of  the  earth's 
magnetic  field.  Both  cases  taken  together  represent  the  extremes  of  the 
range  of  results  to  be  expected  for  propagation  through  ionospheric  fluc¬ 
tuations,  or  ionization  irregularities  caused  by  bariun  cloud  instabili¬ 
ties  or  by  nuclear  detonations. 

The  complex  general  analytic  results  obtained  for  the  two  cases 
of  interest  are  simplified  by  the  use  of  the  thin  phase-screen  approxima¬ 
tion  (zero  slab  thickness)  to  obtain  tractable  expressions  for  the  two- 
frequency  mutual  coherence  function  as  well  as  its  Fourier  transform, 
referred  to  as  the  generalized  power  spectrum. 

The  accuracy  of  the  thin  phase-screen  approximation  is  a  matter 
of  some  importance  to  this  wrk.  This  approximation  is  considered  in  some 
detail  in  this  section  where  comparisons  of  the  case  of  a  thick  layer  to 
the  thin  phase-screen  case  are  presented  for  the  S4  scintillation  index, 
the  mean  time  delay,  and  the  time  delay  standard  deviation. 

4.1  FORMULATION 

Consider  a  monochromatic  spherical  wave  E{'p,z,co,t)  which  origi¬ 
nates  from  a  transmitter  located  at  (0,0, -z^)  and  propagates  in  free- 
space  in  the  positive  z  direction  where  it  is  incident  on  an  ionization 
ir'^egularity  layer  which  extends  from  0  £  z  £  L.  After  emerging  from  the 
layer  at  z  =  L,  the  wave  then  propagates  in  free-space  to  a  receiver 
located  at  (0,0,z^).  This  geometry  is  shown  in  Figure  3-1  in  the 
preceding  section.  As  the  wave  propagates,  its  phase  substantially 
behaves  as  (-i<k>z+ia)t)  so  write 

E(p,z,(D,t)  =  U('p,z,u))  exp{i(u)t-J<k(z' )>dz' ) }  (4-1) 
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(4-2) 


where  <k(z)>  is  the  iiean  wave  number  given  by 

<k{z)>  =  (1  -  =  k  (1  -  kj/k^)^/^ 

where  c  is  the  speed  of  light  in  a  vacuum,  Ne  is  the  mean  ionization 
density,  nc  is  the  critical  electron  density  and  is  related  to  the  clas¬ 
sical  electron  radius  re  by  nc  =  Tt/(x^re),  (re  =  2.82x10"^^  m). 


It  has  been  shown  that  U{p,z,a»)  satisfies  the  parabolic  wave 
equation  under  the  Markov  approximation,  (Tatarskii,  1971;  Yeh  and  Liu, 
1977) 


»ll  o  <Np>  0)2 

V2u  -  2i<k(z)>  ^  -  <k2(z)>  U  =  0  (4-3) 

3Z  (l-0)2/o)2) 

where  ANe  is  a  small  variation  in  the  ionization  level  and  o)p  *  kpC 
is  the  circular  plasma  frequency  of  the  background  ionization.  The 
exponential  time  dependence  has  been  suppressed  and  it  has  been  assumed 
that  o)p  «  0),  otherwise  signal  attenuation  would  be  the  dominant  effect. 

4.1.1  Power  Impulse  Response  Function 

Now  consider  the  case  where  the  transmitted  waveform  is  no 
longer  a  monochromatic  wave,  but  can  be  expressed  as  a  waveform  modulated 
on  a  carrier 


m(t)e 


icoot 


where  m(t)  is  the  transmitted  modulation  waveform  or  transmitted  complex 
envelope  and  uq  is  the  carrier  angular  frequency.  Now  if  M(a))  is  the 
spectrum  of  the  transmitted  complex  envelope. 


M(u))  =7  m(t)e"''“^dt 

»CP 


(4-4) 
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then  the  received  complex  envelope  can  be  expressed  as 


CO 

e(p,z^,t)  =L  jM{w)UCp,z^,a5<-a)Q)e^“^di 


(4-5) 


where  U(p,z^,aj)  is  the  solution  to  the  parabolic  wave  equation  at  the 
receiver  plane  z,,.  This  result  is  proven  in  Section  3  of  this  report. 

The  average  received  power  is  given  by 
I(z^,t)  =  <ee*>  =  f  M(u))M*(aj')  <U('p,Z^,a5<-ii»Q)U*('p,Z^,a)'  +  (DQ)> 

4u2  _oo 


X  e^^“"“'^^daxia)'  (4-6) 

The  two-frequency  mutual  coherence  function,  r,  is  the  quantity  in  the 
angle  brackets.  Under  the  assumption  (to  be  justified  later)  that  r  is  a 
function  only  of  the  transverse  spatial  and  frequency  differences,  Equa¬ 
tion  4-6  may  be  transformed  using  the  sun  and  difference  transformation 

0).  =  U  -  0)' 

d 

^  (i^w' )  (4-7) 

=  A//"<  a>s+o)d/2)M*(  V  Ap=0,Z  ,0)^) 

4tT^  -00 

iu)  jt 

X  e  (4-8) 

Since  the  input  signal  is  m(t),  the  input  power  is  given  by 

I(t)  =  m(t)m*(t)  =  JL  f  rM(u))M*(a)'  )e^(“'‘^'  ^^dajdai'  (4-9) 

4^2  ly 
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Using  the  same  sum  and  difference  transformation,  the  input  power  may  be 
written 

00  ^ 

I(t)  =  J_//M(aj5+ajj/2)M*(«5-toy2)  e  dw^du^  (4-10) 

Now  if  the  input  power  I(t)  is  a  delta  function 


1  i  Jiyt . 

I(t)  =  6(t)  =  ^  /e  dw^ 


(4-11) 


it  is  apparent  that 


00 

1—  M((Og+(o^/2)M*(wg-a>^/2)da^  =  1 

2v  .o» 


(4-12) 


So  for  the  special  case  of  an  input  power  delta  function,  the 
output  power  at  the  receiver  is 


00  Q 

G(Zr’^^  “  T-  /r(AP=0,z^,u)j)e  da,^ 

cTT  «*oo 


(4-13) 


where  the  new  symbol  G  has  been  introduced  for  this  important  function. 

Now  to  find  the  output  power  response  to  a  general  input  I(t), 
first  multiply  Equation  4-10  by  exp(-iu^t)  and  integrate  with  respect  to  t 

ri(t)e  dt  =  J_jf^M(a,5+a,j/2)M*(a)5-u)^/2)e  du^dw^dt 

47[2 


—  /  M(a)5+u)y2)M*((05-(oy2)da)5 

2tr  _oo 


(4-14) 
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The  above  result  conies  about  because  the  t  Integration  yields  a  delta 
function,  2Tr6(a)c|-W(i') .  Now  use  Equation  4-14,  just  derived,  in 
Equation  4-8  to  obtain 


I(z^.t)  =-//l(f)e  r(Ap=0.z^,u^)e  doi^df 

2ir 


= /Kf)  —  7 r(Ap=0,z^,a,^j)e^“‘‘^^  ^  L^df 


=  /l(t' )G(z  ,t-t')dt'  (4-15) 

The  second  expression  above  is  recognized  as  a  convolution.  Thus  the 
output  power  due  to  a  general  input  power  is  given  by  the  convolution  of 
the  input  power  with  the  power  impulse  response  function  G(zr,t). 

The  derivation  above  follows  the  basic  outline  given  in  Ishimaru 
(1978)  and  is  included  in  this  report  to  serve  as  an  aid  to  the  notation 
and  to  emphasize  the  importance  of  the  power  impulse  response  function 
G(zr,t).  It  will  be  seen  that  the  moments  of  G(Zf.,t)  serve  as  useful 
definitions  of  the  mean  time  delay  and  time  delay  jitter  which  may  be  used 
to  characterize  the  severity  of  the  propagation  medium, 

4.1.2  Two-Frequency  Mutual  Coherence  Function 

Thus  the  two-frequency  mutual  coherence  function  r  is  important 
for  the  calculation  of  pulse  propagation  in  a  random  medium  and  it  serves 
as  a  basis  from  which  to  calculate  the  important  power  impulse  response 
function  and  its  moments.  Under  the  Markov  approximation,  r  satisfies  the 
following  equation  (Tatarskii,  1971;  Yeh  and  Liu,  1977). 


IE 

3z 


+ 


i 

2kik2 


(k2v|i-kiv|2)r 


■  [2k i0ik282A(pi-P2)  ”  (k i^8i^+k2^82^)A(0) ]r  =  0  (4-16) 
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where 


r  ~  »yj  ,z,(i)j  )U*(x2  »y2  »z,(i)2)^ 

and  ki  and  k2  are  the  wavenumbers  at  two  frequencies  ui  and  U2> 
respectively,  and  vlj  is  the  two-dimensional  Laplacian 

axi  ayi 

with  a  similar  definition  for  ix2» 

Following  Yeh  and  Liu  (1977)  let 


2/..  2 


tOp  /t^l 


2/..  2 


1  -  OJjj  /(Ol 


(4-17) 


with  a  similar  definition  for  e2» 

The  function  A(p)  is  the  integral  of  the  autocorrelation  func¬ 
tion  of  electron  density  fluctuations,  in  the  direction  of  propa¬ 
gation 


A(pi“P2)  =  y  B5(pi-P2»z' )dz' 


(4-18) 


where 


so  that 


P1-P2  =  (xi-X2,yi-y2) 


5  = 

<Ne> 


A(prP2)  =  2iT//e  ♦.(Kj^.K  =0)d\ 


(4-19) 
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where 


I 


! 


u 

S' 


iL. 

32 

3n2 

32 

3X2 

3Y2 

3C3X 

3n3Y 

I 

I 


1 


1 

I 

I 


r  is  now  written  in  terms  of  the  above  sin  and  difference  arguments  as 
r(?,n.z.ay)  where  ay  =  ckj. 


The  unknown  two-frequency  mutual  coherence  function  may  be 
written  as  r  =  FiTo  where  Fq  is  the  exact  free-space  solution  in  the 
parabolic  approximation 


-iki(xi+yi)  +  ik2{x2+y|) 
2(z+z^) 


(4-22) 


which  under  the  sum  and  difference  transformation  used  previously  becomes 

-ik  (?X+nY)  -  ik  [(XM2)/2  +  {c2+n2)/8]) 

_J - ^  -  (4-23) 

(z+z^)  1 

Substituting  F  =  F^Fg  into  Equation  4-21  and  neglecting  near  zone  terms  of 
the  order  of  ks?^/(2+zt)  and  smaller,  one  obtains 


m 

I 


1. 


i. 
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Mil*- i - - - ... 


+  terms  with  =  0  (4-24) 

3X  3Y 

Equation  4-24  is  valid  in  the  region  0  £  z  £  L  with  boundary 

condition 

ri(5»n»z=0,a)jj}  =  1  (4-25) 


Since  the  boundary  condition  is  independent  of  X  and  Y,  and  the  equation 
has  no  terms  other  than  the  derivatives  with  respect  to  X  and  Y,  it  is 
apparent  that 


Hi  =  ill  =  0 

9X  3Y 


Now  the  substitutions 


z'  =  z  +  z^ 

0  =  ?/z'  =  ?/{z+Z^) 
(^  =  n/z'  =  n/(z+z^) 
yield  the  following 


(4-26) 

(4-27) 

(4-28) 
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(4-29) 


The  additional  substitution  Pi  =  r2r3  where 


r3  =  expj.iA(0)(z.-z,).*(L-yV 


yields  the  simplified  equation  for  T2 


(4-30) 


3r2 

32' 


2  k2  z'2  [ae2  3(|.2j 


-  k** 

[A(z'(e+,j>))  -  A{o)]r2  =  0 

4  k2 


(4-31) 


where  kj  has  been  neglected  with  respect  to  kg.  The  effect  of  this 
assumption  is  to  restrict  the  validity  of  the  solution  to  a  small  range  of 
wavelengths  centered  about  kg. 

Again  following  Sreenivasiah  (1976)  and  Yeh  and  Liu  (1977), 
expand  the  function  A(z'(9+(|»))  in  a  Taylor  series  and  neglect  all  but  the 
first  two  terms  so  that 


A(2'(0+(|)))  =  Aq  +  2' 2(0+,j,)2A2 


(4-32) 


Equation  4-32  is  the  analytic  expression  of  the  quadratic  phase  structure- 
function  approximation. 
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4.1.3 

Aq  and  A2  CALCULATION 

4.1. 3.1 

Isotropic  Irregularities 

The  Ao  and  A2  coefficients  may  be  calculated  for  isotropic 
irregularities  by  performing  the  angular  integration  in  Equation  4-19  to 
obtain 

00 

A(prP2)  -  4Tr2/jo(Kj^p)  =0)  K^^dKj^ 

0  4  ^ 

(4-33) 

where 

=  k/  +  k/ 

P  =  |P1-P2| 

The  Jq  Bessel  function  may  be  expanded  in  a  power  series 
of  equal  powers  of  p  equated  to  obtain 

and  coefficients 

00 

Ao  =  4n2/MK_^,K  =0) 

0  ^ 

(4-34) 

CO 

A2- 

(4-35) 

a)  Gaussian  PSD 

For  a  Gaussian  power  spectrum  of  the  form 

♦  (K)  =  — L  Jiilexpl-  k2Lo^/4} 

?  <Ne>^  8iiV2 

(4-36) 
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with  corresponding  autocorrelation  function 


B(r)  =  exp  {-rVlo^} 
e 


(4-37) 


it  can  be  shown  that 


Aq  =  LoTr^/2  — e^ 

<N  >2 


A,  =  -  I 


1/2  % 


Lo  <V' 


b)  Power  law  PSD 

Consider  a  power  law  spectrim  of  the  form 

^  _  '’Ne(2TT)“3/2)i.(£./L^)(P-i)/3  +  l/L^) 


(4-38) 


(4-39) 


(4-40) 


which  has  a  K"P  behavior  for  1/K  ranging  from  Lq,  the  outer  scale  to 
1^,  the  inner  scale.  The  corresponding  autocorrelation  function  is 


\(P-3)/2 


°Ng  *  *<)  *  *1  ) 

<Np>^(»,/g(p-3)/^K,p.3,/,(«,/g 


(4-41) 


Here  the  exponential  may  be  expanded  as  the  Taylor  cosine  series  since  iji 
is  real.  The  results  for  the  coefficients  Aq  and  A2  are  identical  to 
Equations  4-42  and  4-43  with  the  exception  that  the  three-dimensional  PSD 
for  electron  density  fluctuations  is  replaced  by  the  two-dimensional  PSD. 
The  two-dimensional  PSD  is 
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s„  (K^,K^)  = 

e 


,2  (2.)-ljf(t,/^)(p-3)/2  K(p.,,/  t  K|  t  l/L^) 

— 2 - i _ /M- 


<Ne>2K 


(P-3)/2' 


(*i 


iP-i;/2 


(4-47) 


^2  +  k2  +  1/L2 

X  Z  '  0 


The  tvro-dimensional  autocorrelation  function  corresponding  to  this 
two-dimensional  PSD  is  given  by  Equation  4-41.  It  is  now  easy  to  show 
that  the  values  of  Ag  and  ^2  for  infinitely  elongated  irregularities  are 
identical  to  the  values  for  isotropic  irregularities. 

In  the  work  that  follows,  the  ratio  of  the  constants  Aq  and  A2 
appears  often.  At  this  point  it  is  noted  that  for  a  three-dimensional 
K"P  power  law  electron  density  spectrum 


^  =  -  Lo2(2p-8) 
A2 


p  >  4 


(4-48) 


Remember  that  the  corresponding  one-dimensional  phase  PSD  which  is  often 
used  in  this  report  has  a  power  law  form. 

4.1.4  oj  CALCULATION 


It  is  useful  to  establish  the  relationship  between  the  coeffi¬ 
cient  Aq  and  the  variance  of  phase  fluctuations,  For  a  layer  of 
ionization  az  thick  the  incremental  phase-shift  is  given  by 


k  an  (x,y,z)dz 
A,(x,y)  =  _  J  - - 

--d/2 


(4-49) 


where  we  have  assumed  that  the  critical  electron  density,  n^,  is  much 
greater  than  the  mean  electron  density  and  have  expanded  the  formula  for 
the  index-of-refraction  (neglecting  collisions)  in  a  power  series  and 
retained  only  the  first  two  terms.  Now  from  the  correlation  function  of 
the  phase  fluctuations 
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<A(|)(Xi,yi)A(j)(X2,y2)  = 


k2<N^>2 
0  e 


<aN  (xpyi,z)AN  (X2,y2.z')>  dzdz' 


<N  >2 
e 


(4-50) 


The  double  integration  can  be  reduced  to  a  single  integration  by  a  change 
of  variables  (Papoulis,  1965,  p.  325),  so  that 

Azk2<N  >2  AZ  _ 

<A({i(xi,yi)A<|)(x2»y2)>  = - - —  f  B  (^j-P2,n)  (1  -  in|/Az)  dn 

^  "^2  (4_51) 


If  AZ  is  greater  than  the  correlation  length  of  electron  density  fluctua¬ 
tions,  then  the  contribution  of  the  second  term  is  negligible  and  the 
integration  limits  may  be  changed  to  plus  and  minus  infinity.  Evaluation 
of  the  resulting  expression  at  ■pi-'p2  *  0  yields 


k:<N„>2 

al  =  _L-L_  LA 
4n2  0 

c 


4  p  0  0 


(4-52) 


where  Equation  4-18  is  used  for  Aq  and  L  is  the  total  thickness  of  the 
layer  of  ionization  irregularities.  For  a  K"**  three-dimensional  ioni¬ 
zation  power  spectrum,  Ag  is  given  by  Equation  4-44  and 


ct^2  =  2(xrg)2LLQa2 


(4-53) 


where  is  the  variance  cf  the  electron  density  fluctuations. 


4.1.5  r  SOLUTION 


Now  that  general  expressions  have  been  derived  for  the  coeffi¬ 
cients  Aq  and  A2,  the  solution  for  the  two-frequency  mutual  coherence 
function  may  be  continued.  Introduce  the  additional  substitutions 
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(4-54) 


b=i  kVk/  (4-55) 

4  P  ^ 

V  =  (abA2)^/^  z'  =  a^z'  (4-56) 

„  =  {^]>/*  e  =  a,e  (4-57) 

e  =  (^y/«  4  =  a,*  (4-63) 

With  these  substitutions  Equations  4-31  and  4-32  may  be  combined  to  yield 
II2  -  i  L  [ii.  +  iLl  r2  -  (w2+e2)v2r2  =  0  (4-59) 

3v  v2  L3p^  3e^J 
The  boundary  condition  becomes 

r2(u»e,v=(abA2)  ^/^z^.Ujj)  '•*  1 


An  analytic  solution  of  the  form 
T2  =  f{v)exp(-g(v)(y^e2) } 

may  be  substituted  into  Equation  4-59  to  obtain  the  following  two 
equations,  the  first  consisting  of  terms  independent  of  (p^e^) 

£  +  =  0 

3v 

and  the-  second  consisting  of  terms  with  the  factor  (u^+e^) 

ii  <•  ii?!  t  =  0 

3v 


(4-60) 


(4-61) 


(4  .62) 
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This  pair  of  equations  may  be  solved  exactly  with  the  result 


f(v)  = 


av^cosho>(v-v^)  +  sinha(v-v^) 


(4-63) 


g(v)  =  _  iv  1  -  ovi 
4 


(av^sinha(  v-v^)+cosha(v-v^)\ 
aVj.  cosh  «( V- Vj.  )+s  i  nh  a(  V- Vj. )  / 


(4-64) 


where 


a  =  2e 


i3ir/4 


K-  ' 


=  (abA2) z^=  a^z^ 


Equations  4-30  and  4-60  may  be  multiplied  to  obtain  the  desired 
solution  for  Ti  in  the  ionized  layer  where 

0  <  z  <  L 


Zt  1  1  L  +  \ 

In  transcribing  Equation  4-30,  k^j  has  again  been  neglected  with  respect 
to  kg  to  obtain 

ri(0,(l),z',o)^)  =  exp  [-Ao(z'-z^)kJk2/8kJ] 

X  f(aiZ')  exp  [-a2^(e^(|)2)g(a^z' )  ]  (4-65) 

In  accordance  with  the  definitions  given  by  Equations  4-33  and  4-34,  Aq  is 
used  to  replace  A(0)  in  the  expression  for  Fj. 


To  complete  the  solution  it  is  necessary  to  solve  for  in  the 

region 

L  <  2  < 

or 

L  +  2^  <  2'  <  2^  +  Zt 

Equation  4-65  serves  as  the  boundary  condition  at  z'  =  L+z^.  Since  the 
region  z'  L+z^  corresponds  to  free-spece  with  no  ionization,  Equation 
4-29  is  appropriate  where  the  last  term,  which  is  dependent  on  the  func¬ 
tion  A,  may  be  neglected  since  A  is  zero  in  the  absence  of  ionization 
irregularities.  Also,  in  keeping  with  previous  assumptions,  the 
term  is  ignored  with  respect  to  the  kg^  term  in  Equation  4-29. 

The  Fourier  transform  pair 

"  i(K  6+K  (k) 

ri(0.4.,z‘,o)j)  =//e  ®  ri(KQ.K^,z',a>^)dKgdK^  (4-66) 

,  -i(K  e+K  6) 

fjvKg.K^.z' ,uy)  =  (l/2Tr)^J Je  e  r  ri(e,({>,z',ay)ded(|)  (4-67) 

may  be  substituted  into  the  suitably  modified  Equation  4-29  to  obtain  the 
algebraic  equation 

A 

-  +  --^-1-  (K|+Kpfi  -  0  (4-68) 

3z'  2  k2  z'2  ®  ^ 

Equation  4-68  may  be  solved  and  the  boundary  condition  4-65 
applied  at  z'  =  L+z^  with  the  result 

ri(K9.K^,Zt+z^,u)^)  =  ri(Kg,K^,L+Z^,a)jj) 

X  exp  [-iY(K2+Kp]  (4-69) 
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where 


1  Kd 

2  (L«j)(Zj*z^) 


(4-70) 


The  final  result  may  then  be  obtained  by  taking  the  Fourier  transform  of 

A 

Equation  4-65  to  obtain  ri(Kg,K^,z'+zt ,uy )  and  then  taking  the  inverse 
Fourier  transform  of  Equation  4-69.  The  appropriate  Fourier  transform, 
pairs  are  given  by  Equations  4-66  and  4-67.  The  required  integrals  are 
easily  found  in  Gradshteyn  and  Ryzhik  (1965)  and  the  result  may  be  written 
as 


rj( ~  ri( n»z' =z^+z^, (»j^) 

=  exp  {-o2aj|/2u^}  f[ai(L+z^)] 

( ?^n^)a2^g  [a  i(  L+z^)  ]/  { z^+z^ }  2 
1  +  i4Ya2^g[ai(L+z^)] 

X  {l+i4Ya2^g[ai{L+z^) ]}’^  (4-71) 

where  the  e  and  (j)  coordinates  have  been  transformed  into  the  c  and  n  co¬ 
ordinates  according  to  the  transformation  of  Equations  4-27  and  4-28.  The 
subscript  I  on  r  has  been  introduced  to  denote  that  the  solution  is  appli¬ 
cable  to  the  case  of  isotropic  ionization  irregularities.  The  case  of 
elongated  irregularities  follows.  The  first  exponential  term  in  Equation 
4-71  comes  from  the  first  exponential  term  in  Equation  4-65  where  the 
relationship  between  and  Aq  given  by  Equation  4-52  has  been  utilized 
to  simplify  the  resulting  exponential. 

4. 1.5.1  Elongated  Irregularities 

For  transionospheric  propagation  the  case  of  elongated  irregu¬ 
larities  may  be  more  appropriate  than  that  of  isotropic  ionization  irregu¬ 
larities.  To  model  elongated  irregularities,  assume  that  the  striations 
are  long  tubes,  infinite  in  the  y-direction  shown  in  Figure  3-1. 


I 


\ 

\ 


I 
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For  this  case  the  splution  proceeds  in  the  sane  manner  as  described  for 
the  case  of  isotropic  irregularities. 

The  development  differs  since  there  is  no  y  dependence  in  the 
case  of  elongated  irregularities.  Thus  the  factor  of  4  in  Equation  4-61 
becomes  a  factor  of  2  for  elongated  irregularities.  The  resulting  solu¬ 
tion  for  the  function  f  is  then  the  square  root  of  the  result  obtained  for 
the  isotropic  irregularity  case.  Equation  4-62  for  g  remains  unchanged. 

In  the  succeeding  development,  the  only  difference  occurs  because  one¬ 
dimensional  Fourier  transforms  are  involved  rather  than  two-dimensional 
transforms  as  are  used  for  the  isotropic  irregularity  case.  The  final 
result  for  for  the  case  of  elongated  irregularities  is 

=  exp  |-oJa)|/2a)2}{f[ai(L+2^)]}l/2 

;^a2^g[ai(L+z^.)]/(z^+2^)2 

X  exp  { - 

1  +  i4Ya2^g[ai{L+2^)] 

X  (l+i4Ya22g[ai(L+2^)]}“i/2  (4-72) 

where  f  and  g  are  defined  by  Equations  4-63  and  4-64.  The  subscript  E  on 
r  refers  to  the  elongated  irregularities. 

4. 1.5. 2  Thin  Phase-Screen  Approximation 

Further  simplification  is  obtained  if  the  thick  scattering  layer 
is  replaced  by  an  equivalent  thin  phase-screen  with  infinitesimal  thick¬ 
ness  and  the  same  overall  phase  variance,  As  will  be  seen  the 
accuracy  of  the  thin  phase-screen  approximation  is  a  function  of  the 
propagation  geometry,  the  wavelength  and  the  electron  density  PSD. 

The  thin  phase-screen  approximation  has  been  used  by  many 
authors  and  is  treated  in  detail  by  Mercier  (1962)  and  Salpeter  (1967). 
Wernik,  et  al.,  (1973)  and  Bramley  (1977)  have  examined  the  accuracy  of 
the  thin  phase-screen  approximation  for  the  calculation  of  the 
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scintillation  index.  Rino  (1980)  and  Rino  et  al.  (1981)  have  utilized 
this  approximation  to  obtain  numerical  propagation  results  and  to 
successfully  model  observed  data  from  the  DNA  Wideband  satellite 
experiment. 

Here  the  accuracy  of  the  thin  phase-screen  approximation  is 
examined  in  detail  as  it  affects  the  calculation  of  the  scintillation 
index  and  the  mean  time  delay  and  time  delay  jitter.  It  will  be  seen  that 
the  maximum  error  in  the  use  of  the  thin  phase-screen  approximation  is 
small;  therefore  this  useful  approximation  is  invoked  here  to  simplify  the 
mathematics  and  to  aid  in  understanding  the  results. 

Now  taking  the  limit  as  the  layer  thickness  approaches  zero  and 
utilizing  the  substitutions, 


-A  (z.+z  )2 
5  0  t  r 


(4-73) 


(4-74) 


simplified  expressions  are  found  for  the  two-position,  two-frequency 
mutual  coherence  function  in  the  thin  phase-screen  approximation 


rj(i;,n,z^,wj)  =  exp  {-aJw|/2(o2} 


X  exp 


-{?2+r,2)/il2 

^•^’V“coh 


(4-75) 


i 


(4-76) 


Thus  the  thin  phase-screen  approximation  allows  the  complicated 
expressions  4-71  and  4-72  to  be  written  in  terms  of  the  simplified  param¬ 
eters  Aq  “coh*  following  it  will  be  seen  that  Hq*  signal 

decorrelation  distance,  and  ucoh*  coherence  bandwidth,  are  two  very 
important  parameters  which  describe  the  characteristics  of  the  signal 
fluctuations. 


Figure  4-1  shows  the  thin  phase-screen  model  of  the  general 
propagation  geometry  shown  in  Figure  3-1.  The  distance  from  the  trans¬ 
mitter  to  the  phase-screen  is  and  the  additional  distance  to  the 
receiver  is  z,..  The  total  free-space  propagation  distance  is 
Z|.+Zf..  It  is  apparent  from  Equation  4-74  that  the  coherence  bandwidth 
is  a  reciprocal  quantity.  That  is,  an  interchange  of  transmitter  and 
receiver  does  not  affect  its  value.  However  the  decorrelation  distance  jIq 
is  not  reciprocal.  Zq  is  a  measure  of  the  average  distance  between  fades 
at  the  receiver  and  depends  on  the  path  geometry. 

A  more  useful  measure  is  the  signal  decorrelation  time,  tq* 
which  is  a  measure  of  the  inverse  fading  rate  or  inverse  fading  band¬ 
width.  tq  is  a  time-domain  measure  of  fading  whereas  to  ®  spatial 
measurement.  The  signal  decorrelation  time  is  given  by  the  equation 


(4-77) 


where  Vgff  is  the  relative  velocity  of  the  1 ine-of-sight  at  the  receiver 
location.  Although  this  velocity  is  a  function  of  transmitter  and  re¬ 
ceiver  motion  as  well  as  irregularity  motion,  it  is  sufficient  to  consider 
the  case  where  the  transmitter  and  receiver  are  stationary  and  the  phase- 
screen  is  in  motion.  Assume  that  the  scattering  layer  consists  of 
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Figure  4-1.  Thin  phase-screen  propagation  geometry. 


striations  moving  in  unison  at  a  velocity  V  perpendicular  to  the  line-of- 
sight.  Then  the  projection  of  these  striations  at  the  receiver  location 
due  to  a  signal  originating  at  has  velocity 


V(z  +zj 
'  r  t' 


'eff 


(4-78) 


Combining  Equations  4-73,  4-77  and  4-78,  the  received  signal  decorrelation 
time  is  obtained  as 


’j' 

% 

'I 


i 


I 


(4-79) 


It  is  apparent  that  tq  is  only  a  function  of  the  layer  velocity  and  is 
reciprocal  as  it  should  be. 


For  the  often  used  in  situ  power  spectrum  of  one-dimensional 
electron  density  fluctuations,  Aq/A2  may  be  obtained  from  Equations  4-44 
and  4-45  and  the  signal  decorrelation  time  is  given  by 


4.2  GENERALIZED  POWER  SPECTRUM 


(4-80) 


At  this  point  it  is  convenient  to  introduce  the  generalized 
power  spectrim  which  is  defined  as  the  Fourier  transform  of  the  two- 
frequency  mutual  coherence  function.  The  generalized  power  spectrum 
serves  two  useful  purposes.  First,  it  is  the  basic  function  which  is  used 
in  statistical  signal  generation  as  discussed  in  Section  3.  Second,  it 
acts  as  an  intermediate  function  from  which  to  easily  determine  the  power 
impulse  response  function.  This  latter  aspect  is  considered  in  detail  in 
the  following  subsection. 

For  the  cases  of  isotropic  and  elongated  irregularities,  the 
generalized  power  spectra  are  defined  as 


Si(K  ,K  ,Zr,T)  = 

X  exp  [-i(  ^+rK  ^)+i(»)jT]d5dnduj  (4-81) 


Se(Kj.,z^,t)  =  (2Tr)-2jr /rE(c»V“d^  (4-82) 


The  inverse  transforms  are  given  by 


rj(c,n,2^.(0jj)  =  jr//ST(K  ,K  ,z^,t)  exp  [+i(i;K  +nK  )-iu)^T]dK  dK  dx 

Si  si  si 


(4-83) 
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(4-84) 


r£( t)  exp  [+i?K  -tUjjt]  dK  dx 

oOBsaA  *•  *>  S 


In  the  case  of  elongated  irregularities,  the  one-dimensional 
generalized  power  spectrun  is  given  by  Equation  4-82.  For  the  case  of 
isotropic  irregularities,  the  one-dimensional  generalized  power  spectrum 
is  given  as  the  integral  over  the  two-dimensional  spectrum  (Rufenach, 
1975). 

Si{K  ,z^,t)  =  JSj(K  ,K  ,z^,T)dK  (4-85) 

Equation  4-81  is  easily  integrated  with  respect  to  to  obtain 

Si(K  ,2r,T)  =  rT(4»’i=0,z  ,t)  exp  [-i^K  +iu)jT]dcdio.  (4-86) 

where  rj  is  evaluated  at  zero  n  because  the  integration  yields  a  delta 
function.  The  inverse  transform  of  Equation  4-86  is 

00  00 

rj(c,n=0,z^,T)  =/J'Sj(K  ,z^,t)  exp[icK  -iu.xjdK  dx  (4-87) 

The  integrals  of  Equations  4-81  and  4-82  can  be  computed  by 
recognizing  that  the  interior  integrals  in  these  two  equations  with 
respect  to  ?  and  n  are  related  to  quantities  computed  above.  That  is, 
from  Equation  4-81 

{2ir)"^jfjr  rj(i:,n,z^,u)^)  exp[-i(!;K  +tiK  )]dcdn 

•?oe«?oo  **  ' 

=  ri(Kj,.K^,z^,u)j) 

=  (z^+zp2pj(K^^K^^z  =  (4.88) 
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where  Tj  is  given  by  Equation  4-69.  The  (zt+z^)^  term  results  from  the 
transformation 

e  =  c/(VV 

♦  =  'l/(Zj+2^) 

relating  angles  and  distance  perpendicular  to  the  direction  of 
propagation. 

In  the  case  of  infinitely  elongated  irregularities  it  can  simi¬ 
larly  be  shown  that  the  interior  integral  of  Equation  4-82  may  be  written 
as 

00 

(l/2Tt)  /r£(?,z^,uj)  exp[-i;K  Id? 

-00  ^ 

A 

where  Pi  in  Equation  4-91  is  the  one-dimensional  equivalent  of  the  two- 
dimensional  result  given  by  Equation  4-69. 

Equations  4-88  and  4-91  may  be  substituted  into  Equations  4-81 
and  4-82  to  obtain  simplified  expressions  for  the  generalized  power 
spectra 

00 

Si(K^,K^,z^,t)  =  (1/2ti)  /Ti(Kj..K^,z^,u)j)  exp( iw^jTldw^  (4-92) 

00 

SglK  ,Zr,T)  =  (l/2it)  f  r^CK  exp(ia)^T)du^  (4-93) 


(4-89) 

(4-90) 
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The  Fourier  transforms  of  Equations  4-92  and  4-93  cannot  be  per 
formed  analytically  except  for  the  case  of  a  thin  phase-screen.  In  this 
special  case  it  is  possible  to  show  directly  from  Equation  4-69  and  its 
equivalent  for  elongated  irregularities  that 


p  exp|-o^»|/2»|)  exp  i  )*o*l"“d/"coh>  I 


(thin  phase-screen) 


(4-94) 


*o<‘"V“coh>l 

(thin  phase-screen)  (4-95) 

In  the  thin  phase-screen  approximation  the  preceding  four  equations  can  be 
combined  to  yield  for  the  generalized  power  spectra 


Si(K  ,K  ,z  ,t)  exp{-(K2+K2)^2/4 

^  5  ^  25/2^3/2  (  ^  ° 


(thin  phase-screen) 


(4-96) 


S  (K  .z  ,t)  =  exp  j  -  K2£2/4  .  1 

‘x  2^/2^  (  2 


"coh"^ 


K2t2/4]- 


(thin  phase-screen) 


(4-97) 


where  Hq  and  cocoh  ^re  given  by  Equations  4-73  and  4-74.  The  parameter  a 
is  given  by  the  equation 


a 


“coh 


(4-98) 


and  will  be  shown  to  be  an  important  measure  of  the  ratio  of  the  propaga¬ 
tion  distance  to  the  focal  distance  of  the  largest  irregularities. 
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The  one-dimensional  generalized  power  spectrum  for  isotropic 
irregularities  is  given  by  Equation  4-85;  for  the  thin  phase-screen 
approximation  of  Equation  4-96  the  integral  is  obtained  from  the  tables  as 

.  -Pi-  -  I  -  “coh^l^l 


X  exp(s2/4a2)Kj^^(g2/4jj2) 


(4-99) 


where 


6  =  a2[K2£2/4  .  +  1 


(4-100) 


Equations  4-97  and  4-99  are  the  one-dimensional  forms  of  the 
generalized  power  spectra  which  are  used  in  the  statistical  generation  of 
realizations  of  the  signal  received  after  propagation  through  a  strongly 
turbulent  medium.  Appendix  A  shows  plots  of  these  functions  for  a  range 
of  values  of  the  parameter  o.  It  is  important  to  note  that  the  one-dimen¬ 
sional  generalized  power  spectrum  is  required  for  the  case  of  isotropic 
irregularities  since  statistical  realizations  of  the  received  signal  are 
desired  that  are  a  function  of  one  spatial  dimension  only.  Generation  of 
these  realizations  is  discussed  in  detail  in  Section  3  of  this  report. 

4.3  POWER  IMPULSE  RESPONSE  FUNCTION  SOLUTION 

In  this  subsection  the  impulse  response  function  due  to  a  trans¬ 
mitted  power  delta  function  is  determined  in  the  thin  phase-screen  approx¬ 
imation.  Equation  4-13  gives  the  power  impulse  response  function  as 


G(z  >t)  =  (1/2it)  r(4p=0,z  ,u).  )e  ^  dt 


(4-13) 


Now  it  is  easy  to  prove  from  Equation  4-86  that  for  zero  4p 


Gi(Zr,T)  =jrSj(K  ,z^,T)dK^ 


(4-101) 
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Equation  4-101  is  true  because  the  integration  on  the  right-hand  side 
of  Equation  4-86  yields  the  delta  function  2tt6(?)  which  is  then  integrated 
with  respect  to  ?  to  yield  the  right  hand  side  of  Equation  4-13.  An  ex¬ 
pression  identical  to  Equation  4-101  is  found  for  the  case  of  elongated 
irregularities.  Thus  Equation  4-101  may  be  applied  to  the  generalized 
power  spectra  given  by  Equations  4-97  and  4-99  to  analytically  determine 
the  power  impulse  response  functions: 


Gl(Zr»'f)  = exp  {-(a(0(.Qj,T)  V2 }  exp  {(1-o2u)|.q^t) 


(4-102) 


®E^^r’^^  =  ^  exp  (-(aa>^Qj^T)2/2}(l-a2a)^Q^x)l/2 
ztt 

X  exp  ((l-a^w^^^x) V4a2}  <1/4  ((l-a^u^^^T) V4o2}  (4-103) 

where  4i(.(x)  is  the  complementary  error  function 

CO 

4c(x)  =  4:/exp  (-y2)dy 

/tt  X 

and  Ki/4(x)  is  the  modified  Bessel  function  (Abramowitz  and  Stegun,  1965). 
4.3.1  Limiting  Values 

It  is  apparent  that  the  power  impulse  response  function  G 
depends  upon  the  parameter  o  defined  by  Equation  4-98.  Using  Equation 
4-74  for  it  can  readily  be  shown  that 


-A2Xz^z^. 


(4-104) 


Thus  a  is  a  reciprocal  quantity  and  does  not  depend  on  the  relative  place¬ 
ment  of  the  transmitter  and  receiver. 


For  a  Gaussian  irregularity  spectrun  the  values  of  Aq  and  A2  are 
obtained  fron>  Equations  4-38  and  4-39  and 


(Gaussian  PSD) 


(4-105) 


For  a  K"**  in  situ  three-dimensional  electron  density  spectrum  or  a 
one-dimensional  phase  PSD,  Equations  4-44  and  4-45  give  Ag  and  A2  with  the 
result 


(K-**  Ng  PSD) 


(4-106) 


It  is  apparent  that  a  is  a  measure  of  the  severity  of  the  scattering  since 
it  is  directly  proportional  to  the  phase  standard  deviation.  Recalling 
from  Section  2.3  that  the  focal  distance  of  a  Gaussian  lens  is  given  as 


where  Lq  is  the  scale  size  and  (jjg  is  the  peak  phase,  it  is  seen  that  o  is 
proportional  to  Z/F  where  Z  is  the  effective  propagation  distance.  Here 
the  peak  phase  for  a  single  lens,  (j>o,  and  the  phase  standard  deviation, 
0^,  are  taken  as  roughly  equivalent.  Thus  large  values  of  o  correspond 
to  situations  where  Z  »  F  and  phase  fluctuations  have  been  fully  con¬ 
verted  to  amplitude  fluctuatic":,  through  the  process  of  diffraction. 

Small  values  of  a  correspond  to  situations  where  phase  fluctuations  have 
not  yet  been  converted  to  amplitude  fluctuations,  as  would  occur  close  to 
a  random  phase-screen.  Thus  the  value  of  a  is  a  measure  of  the  relative 
importance  of  scatter. ng  and  dispersion  with  large  values  indicating 
strong  scattering  effects  and  small  values  indicating  strong  dispersive 
effects. 


176 


In  order  to  calculate  the  power  impulse  response  function  in  the 
large  a  limit  it  is  convenient  to  write  Equation  4-101  as 


«bOO«*OQ  ^ 


(4-107) 


where  Sj  is  given  in  the  thin  phase-screen  approximation  by  Equation 
4-96.  For  large  values  of  a  it  can  be  shown  (Stakgold,  1967,  p.21) 

lim^exp  {-a^x^l  =  6(x) 

a-^oo 

Thus  for  large  a  the  generalized  power  spectrum  of  Equation  4-96  is  given 
by 

^,2 

lim  S,(K  ,K  ,z  .t)  = 

1  Vi  n  • 


a-xo 


I- 


In  this  case  the  integration  specified  by  Equation  4-107  can  easily  be 
performed  with  the  result 


lim  G,(z  ,t)  = 


a-xio 


“coh  “P  <-“coh’>  •  ’  >  8 


T  <  0 


(4-109) 


For  the  case  of  infinitely  elongated  irregularities  it  is  recognized  that 
the  generalized  power  spectrum  can  be  written  as 


lim  Sp(K  ,z  ,t)  =  — ■ 


exp  {-  k2ji2/4} 


^  (toh^^  “  K5«o^/4) 

(  /2  ) 


(4-110) 


With  this  simplification  in  mind  the  power  impulse  response 
function  is  easily  found  from  Equation  4-101 

(4-111) 

0  ,  T  <  0  ■ 

Equations  4-109  and  4-111  give  the  power  impulse  response  func¬ 
tion  for  the  case  of  large  a  where  scattering  and  diffraction  effects  dom¬ 
inate  dispersive  or  phase  effects. 

In  the  limit  of  small  a,  corresponding  to  the  case  that  disper¬ 
sion  dominates  over  diffraction,  the  generalized  power  spectra  of 
Equations  4-97  and  4-99  take  on  the  identical  form 

lim  S  (K  ,z  ,t)  =  lim  S  (K  ,z  ,t) 

=  exp  1-  K|«|/4  -  <.2»c|hTV2)  (4-112) 

Equation  4-112  follows  directly  from  Equation  4-97  for  S^.  For  the  case 
of  isotropic  irregularities.  Equation  4-112  follows  from  Equation  4-99 
where  use  is  made  of  the  fact  that 


lim  (exp  (B2/4a2)  Kj/JgV4a^)}  =  /2n 

Equations  4-101  and  4-112  may  then  be  combined  to  obtain  the  power  impulse 
response  function  in  the  limit  that  dispersion  is  dominant.  That  is. 


Gt(2  s 

r  r 


=  lim  Gp(z  ,t) 
a-K) 


(4-113) 
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It  is  noted  that  Equation  4-13  gives  a  Fourier  transformation 
relationship  between  r  and  6.  The  inverse  transform  may  be  evaluated  at 
zero  wd  with  the  result 


r{Ap=0,z^,a)^=0)  =/G(z^.T)dT  (4-114) 

•CO 

Since  r(0,z^,0)  is  unity,  the  above  equation  provides  a  check  on  the 
derivation.  It  is  easily  shown  that  all  the  expressions  for  the  impulse 
response  given  here  obey  the  normalization  condition  of  Equation  4-114. 

Figures  4-2  and  4-3  show  the  power  impulse  response  function  for 
values  of  a  of  1,  3,  10,  and  100  for  the  case  of  isotropic  irregularities 
and  for  o  values  of  1,  3,  and  10  for  the  case  of  elongated  irregularities. 
These  curves  are  plotted  directly  from  Equations  4-102  and  4-103;  conse¬ 
quently  the  figures  strictly  apply  to  cases  where  the  thin  phase-screen 
approximation  is  valid.  However,  as  will  be  seen  in  the  following,  this 
restriction  is  a  minor  one. 


is: 
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Figure  4-3.  Power  impulse  response  function  for  elongated  irregu¬ 
larities  in  the  thin  phase-screen  approximation. 


-3  -2 


Figure  4-4.  Power  impulse  response  function  for  isotropic  and  elongated 
irregularities  in  the  thin  phase-screen  approximation. 


In  order  to  simplify  the  presentation  of  these  results,  the 
power  impulse  response  function  Is  normalized  by  dividing  by  the  coherence 
bandwidth  and  choosing  ucoh^  the  abscissa.  In  both  figures  it  is 
apparent  that  the  power  impulse  response  function  changes  from  a  roughly 
Gaussian  shape  at  unity  a  to  an  exponential  shape  for  larger  a.  In  the 
case  of  isotropic  irregularities,  the  power  impulse  response  function  is 
very  close  to  its  asynptotic  exponential  form  for  a  value  of  a  of  100. 

Figure  4-4  shows  the  power  impulse  response  function  for  the 
case  of  ct  £  0.01.  Note  that  a  different  convention  is  used  for  plotting 
purposes  than  is  used  in  the  previous  two  figures.  For  this  range  of  a 
dispersive  effects  are  dominant  and  the  isotropic  and  elongated  irregular¬ 
ity  cases  are  identical.  For  a<  0.01  it  is  evident  that  the  small  a, 
asymptotic  form  of  the  power  impulse  response  function,  given  by  Equation 
4-113,  is  reached  for  a  =  0.01.  Equation  4-98,  which  defines  a,  may  be 
rewritten  as 

=  laJ 

coh  0  <|) 

For  most  cases  of  interest  involving  lossless  propagation  through  ioniza¬ 
tion,  the  carrier  frequency  is  large  with  respect  to  the  phase  standard 
deviation  and  the  right-hand  side  of  the  above  expression  is  large.  In 
that  case  the  limit  of  small  a  corresponds  to  the  limit  of  large  coherence 
bandwidth,  wcoh* 


1  im  G,(z  ,t)  lim  G,(z  ,t) 
r  r’  '  r  r 

“coh"” 


(4-115) 


The  expression  on  the  "ight-hand  side  of  Equation  4-115  is  recognized  from 
Equation  4-113  as  a  Dirac  delta  function  (Stakgold,  1967,  p.  21),  which 
may  be  written  as 


lim  Gj(z^,t)  =  6(t) 
“coh"“ 


(4-116) 


Thus  in  the  small  a  limit,  dispersion  or  phase  effects  dominate,  and  the 
power  impulse  response  function  is  a  delta  function. 


4.4  MEAN  TIME  DELAY  AND  TIME  DELAY  JITTER 

Two  important  measures  of  the  effect  of  scintillation  on  wide 
bandwidth  signals  are  the  mean  time  delay  and  the  time  delay  jitter. 

The  mean  time  delay  and  time  delay  jitter  of  the  received  signal  may  be 
obtained  (Yeh  and  Liu,  1977)  or  these  measurements  may  be  applied  directly 
to  the  power  impulse  response  function  to  obtain  results  which  are  indeed 
measurements  of  the  propagation  channel.  With  the  latter  approach  in  mind 
define 

<T>  =  /TG(2^,T)dT  (4-117) 


0^  =  y'  T^G(z^,T)dT  -  <T>^  (4-118) 


The  mean  time  delay  and  the  time  delay  jitter  are  functions  of  the  re¬ 
ceiver  location  but  this  functional  dependence  is  omitted  here  for  economy 
of  notation.  The  mean  time  delay  and  the  time  delay  jitter  are  simply 
related  to  the  two-frequency  mutual  coherence  function  since  from  Equation 
4-13 

G(z^,t)  =  L  jf  r(c=0,ti=0,Z  ,a)jj)e  ^  du),  (4-13) 

2t(  .00 


Both  integrals  4-117  and  4-118  may  be  performed  directly  since  the  t  inte¬ 
grations  yield  Dirac  delta  functions  or  a  derivative  of  a  Dirac  delta 
function.  The  results  are 


<T> 


ar 

3o)  j  0)  gO 


0 


2 


T 


au)^2  ^^=0  aa)jj  u)jj=0 


(4-119) 

(4-120) 
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Equations  4-119  and  4-120  are  general  results  and  apply  to  both  elongated 
and  isotropic  irregularities  and  are  not  limited  to  the  case  of  a  thin 
phase-screen. 

4.4.1  Mean  Time  Delay  <t> 


Equation  4-119  may  be  easily  applied  to  the  two-frequency  mutual 
coherence  functions  given  by  Equations  4-71  and  4-72  and  results  obtained 
that  are  valid  for  the  thick  layer  geometry  shown  in  Figure  3-1.  After 
much  algebra  the  results  are  found  to  be 


<T,>  = 


3I<«cA„ 
0  0 


{Zt+z^)(z^+L) 


(4-121) 


<Tp>  =  i  <Tt> 

E  2  ^ 


(4-122) 


An  intuitive  explanation  for  the  factor  of  one-half  in  the  relationship 
between  the  mean  time  delay  for  isotropic  and  elongated  irregularities  is 
that  the  elongated  irregularities  cause  scattering  in  only  one  plane 
whereas  isotropic  irregularities  cause  scattering  in  both  planes 
perpendicular  to  the  propagation  path. 

It  is  noted  that  the  expression  given  in  Equation  4-121  is 
reciprocal.  Interchange  of  transmitter  and  receiver  is  accomplished  by 
the  substitutions 


z'  =  z  -L 
r  r 


(4-123) 


=  Z..+L 


(4-124) 


where  the  prime  quantities  are  the  values  after  the  interchange.  These 
equations  apply  to  the  thick  scattering  layer  geometry  shown  in  Figure 
3-1. 


For  a  transmitted  plane  wave  z^  increases  to  infinity  and 
the  mean  time  delay  simplifies  to 


(4-125) 


lim  <Tj>  = 

Zt^- 


k^cA 
0  0 


{2z^-L) 


This  result  can  be  shown  to  be  identical  to  Equation  54  of  Yeh  and  Liu, 
1977  neglecting  contributions  from  the  transmitted  signal  in  the  referenc¬ 
ed  work.  Now  for  a  transmitted  plane  wave  and  the  thin  phase-screen 
approximation,  L  in  Equation  4-125  is  neglected  and 


lim  <T,>  (4-126) 

z.>«  ^  k2cA^ 
t  0  0 


where  the  receiver  location  in  the  thin  phase-screen  approximation  is 
written  as  z^^^^  to  indicate  that  the  receiver  location  may  depend  on 
the  thickness  of  the  scattering  layer  to  be  represented  by  an  equivalent 
thin  screen.  One  logical  location  for  a  thin  phase-screen  chosen  to 
represent  a  thick  layer  is  the  center  of  the  original  layer.  In  this  case 


^reff  ^r  ’ 


(4-127) 


For  this  choice,  and  the  case  of  an  incident  plane  wave,  the  phase-screen 
approximation.  Equation  4-126,  gives  a  result  identical  to  the  exact 
calculation.  Equation  4-125. 


In  the  more  general  spherical  wave  case  the  thin  phase-screen 
limit  is  obtained  by  setting  the  layer  thickness  L  equal  to  zero  in 
Equation  4-121  and  taking  z^  and  z^  as  and  Zj.g^^  to  indicate  that 

these  values  depend  on  the  original  layer  thickness.  The  result  may  be 
written  as 


<Ti> 


heff  heff 


^heff  ^reff^ 


(4-128) 
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which  is  reciprocal  as  it  should  be.  Now  to  locate  a  thin  phase-screen  at 
the  center  of  a  layer  of  thickness  L  let 


teff  =  ^ 

(4-129) 

reff  °  \ 

(4-130) 

It  is  quickly  shown  that  the  resulting  expression  for  the  mean  time  delay 

<,  >  =  IIM 

2l<|cAo  \ 

is  not  identical  to  the  original  thick  layer  result,  Equation  4-121,  in 
the  general  spherical  wave  case. 

In  summary,  the  thin  phase-screen  result  for  the  mean  time-delay 
is  exact  for  plane  wave  propagation,  but  not  for  the  more  general  case  of 
spherical  wave  propagation.  The  full  expression  given  by  Equation  4-121 
is  valid  for  a  thick  layer  and  spherical  wave  propagation. 

Another  choice  of  a  thin  phase-screen  location  is  discussed  by 
Fante  (1981)  where  it  is  shown  that  the  location  may  be  chosen  as  a 
function  of  the  PSD  to  give  the  exact  result  for  However  this 
choice  of  location  does  not  give  the  correct  result  for  the  two-position 
mutual  coherence  function. 


(4-131) 


4.4.2  <T>  and  -  Thin  Phase-Screen  Approximation 

In  the  thin  phase-screen  approximation  the  simplified 
expressions  for  Pj  and  given  by  Equations  4-75  and  4-76  may  be 
used  in  Equations  4-119  and  4-120  to  compute  the  mean  time  delay  and 
time  delay  jitter.  The  results  are  the  following  expressions 


(4-132) 
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1 


(4-134) 


where  the  parameters  and  o  are  defined  in  Equations  4-74  and  4-98 
respectively.  In  the  following  these  simple  expressions  for  the  thin 
phase-screen  approximation  are  compared  to  exact  results  for  a  thick 
scattering  layer.  It  will  be  shown  that  the  thin  phase-screen  calcula¬ 
tions  of  the  mean  time  delay  and  the  time  delay  jitter  are  valid  under  a 
wide  range  of  conditions  and  are  therefore  very  useful. 

4.5  VALIDITY  OF  THE  THIN  PHASE-SCREEN  APPROXIMATION, 

In  the  following  sections  of  this  report  the  validity  of  the 
thin  phase-screen  approximation  is  investigated  for  calculations  of  the 
scintillation  index  in  the  weak-scatter  approximation  (Salpeter»  1967)  and 
for  the  mean  time  delay  and  time  delay  jitter  in  the  strong  scatter 
limit.  The  scintillation  index  calculation  is  included  here  because  this 
result  is  important  for  the  statistical  signal  generation  technique 
described  in  Section  3  of  this  report.  The  results  indue  !  here  are 
somewhat  similar  to  those  developed  by  Wernick,  et  al.  (1973)  with  the 
difference  that  the  results  shown  here  model  a  thick  layer  with  a  cen¬ 
trally  located  phase-screen  as  opposed  to  a  phase-screen  located  at  the 
layer  bottom.  It  can  easily  be  demonstrated  that  the  centrally  located 
phase-screen  gives  better  agreement  with  exact  weak-scattering  calcula¬ 
tions. 


The  geometry  under  consideration  is  shown  in  Figure  3-1  with  the 
single  exception  that  the  transmitted  wave  is  assumed  plane  so  that  the 
transmitter  is  located  at  infinity.  This  transmitter  placement  is  solely 
for  analytical  convenience  in  this  subsection. 
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4.5.1 


S^-Rytov  Approximation 

From  Tatarskii  (1967,  Equation  7.21)  the  perturbation  of  the 
logarithm  of  the  electric  field  in  the  Rytov  approximation  is  given  as 


i|;(r)  = 


W 


-ik|r-r' I 

An(r’)  dV 

|r-r' I 


(4-136) 


The  unperturbed  incident  field  u^(r)  is  here  taken  as  a  plane  wave  travel¬ 
ing  in  the  positive  z  direction 

UQ(r)  =  exp(-ikz}  (4-137) 


An(r)  is  the  variation  in  the  index-of-refraction  and  the  integration  is 
over  the  entire  extent  of  the  fluctuations. 


As  usual  the  factor  jr-r|  in  the  exponential  is  expanded  in  the 
Fresnel  approximation  as 


|r-r' I 


(4-138) 


and  as  z-z'  in  the  denominator.  Expression  4-138  is  valid  as  long  as  the 
wavelength,  a,  is  less  than  the  inner  scale  size,  z,  so  that  the 
scattering  of  even  the  smallest  eddies  is  in  the  forward  direction  (Fante, 
1975). 


Substituting  the  incident  field  into  the  integral  expression 
given  by  Equation  4-136,  one  obtains 

♦(0,0,z)  =  fmiv)  dv  (4-139) 

271  J  (z-z  ) 


The  mean-square  log-amplitude  fluctuation  is  given  as 
=  <(Reii))^> 


(4-140) 


and  depends  upon  the  correlation  function  of  the  index-of-refraction 
fluctuations  which  can  be  written  as 


<An(r 


■')An(r'')>  =  j 


S  (K)e 
n 


iK*(r'-r" 


(4-141) 


Now  assume  that  the  irregularities  are  uncorrelated  in  the  direction  of 
propagation  (this  simplifying  assumption  is  equivalent  to  the  Markov 
approximation;  Fante,  1975)  so  that 

<An(r')An(r")>  =  6(2'-z")  A(7'-p".z‘) 


*  _  iK_L*(p'-p")  _ 

=  2Tr6(z'-z“)/s^(Kj_.K^=0;z')e  dKj_  (4-142) 

•»eo 


Here  p  and  Kj.  are  the  vector  and  associated  wavenumber  in  the  plane  niormal 
to  the  direction  of  propagation.  It  is  evident  that  the  spectra  Sj^(K)  and 
S  (i<,,K  =0;z')  must  be  both  real  and  even  for  Equations  4-141  and  4-142  to 
bS  valii 


Thus  from  Equation  4-139,  the  log-amplitude  fluctuation  x  or  the 
real  part  of  ij)  may  be  written  as 


k2  AM(r')  cos[k(x'^+y'^)/2(z  -z')] 
x(0.0.z,)  =  I-  j  - dV  (4.143) 


Utilizing  Equation  4-142  for  the  correlation  function  of  the  index-of- 
refraction  fluctuations,  one  obtains  the  mean-square  log-amplitude  fluctu¬ 
ation 


2tt 


I  I  dK^S„(K„K^=0;z-)e 


iK,-(p'-?') 


X  cos[k(x'^-y'^)/2(z^-z')]  cos[k(x"^+y"^)/2(z^-z") ]  (4-144) 


188 


The  exponential  factor  may  be  expanded  as  cos[Kj^(x'-x")+K  (y'-y") ]  and  all 
integrations  may  be  easily  done  with  the  aid  of  a  table  bf  integrals  to 
obtain 


(p-  =  2itk^  J* dz'  jf 

X  Q  ^00 


dKj^  S^(K 


0;z’)  sin2[(z  -z')(K2+K2)/2k]  (4-145) 


X  y' 


Now  a  general  mathematical  description  for  three-dimensional 
index-of-refraction  fluctuations  subject  to  a  Gaussian  power  spectrum  is 
given  by 


Sn(K;z')  =  e  exp{-  K2+e2K2+K2)L2/4}  (4-146) 

where  e  represents  an  elongation  factor.  The  limit  as  c  approaches 
infinity  corresponds  to  the  case  of  elongated  irregularities.  A  value  of 
unity  for  e  corresponds  to  isotropic  fluctuations  perpendicular  to  the 
direction  of  propagation.  Note  that  the  normalization  here  is  such  that 

7  S  (Kiz')dK  »  (4-147) 

•  00 

For  the  case  of  isotropic  irregularities  the  integral  can  be 
performed  by  transforming  to  polar  coordinates.  For  elongated  irregulari¬ 
ties  the  integration  is  more  complex  and  the  limit  as  e  tends  to  infinity 
must  be  taken  after  performing  the  K  integration.  Then  perform  the 
z'  integration  and  finally  the  integration  to  obtain 

S4E  =  2o2  |l  +  (1+16(C2-?l)^)^''‘*  sin  (i  tan-^  4(?2-^l)) 

-  (l+iec^)^/**  sin(i  tan-^  452)]!  (4-149) 
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where  the  weak  scattering  approximation 

Si  =  4o2 
4  X 


(4-150) 


has  been  utilized  to  obtain  the  scintillation  index  and  the  substitution 


a2  =  ✓?  k2LL  o2 

(|>  on 


(4-151) 


which  is  valid  for  a  Gaussian  PSD  has  been  utilized.  The  subscripts  I  and 
E  refer  to  the  case  of  isotropic  and  elongated  irregularities, 
respectively.  For  ease  of  notation  the  additional  substitutions 


=  ^L/2,tL2 


(4-152) 

(4-153) 


have  been  incorporated  into  Equations  4-148  and  4-149.  Equations  4-148 
and  4-149  are  the  expressions  for  the  scintillation  index  for  a  thick 
layer  in  the  Rytov  approximation. 

4.5.2  S*-Thin  Phase-Screen  Approximation 

The  thin  phase-screen  results  for  the  scintillation  index  may  be 
obtained  directly  from  Salpeter  (1967)  as 


X  y  e 


(4-154) 


5|e  =  "  fW  dK 


(4-155) 


where  Zg  is  the  distance  from  the  phase-screen  to  the  observer  and  the 
phase  power  spectral  density  is  given  by 


(4-156) 
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(4-157) 


Equation  4-157  Is  the  one-d1fnensiona1  equivalent  to  the  poMer  spectral 
density  given  by  Equation  4-156.  The  final  results  for  the  scintillation 
Index  In  the  weak-scatter,  thin  phase-screen  approximation  are  then 

Sjj  .  2o2  [1  -  (l+16c|)-^]  (4-158) 

SJe  -  2o2  [1  -  (1+165|)-i/-  cos(itan-i  45g)]  (4-159) 

Mhere 

-  nJZtll  (4-160) 


In  the  following  subsection,  the  results  for  the  scintillation 
Index  In  the  thick  layer  Rytov  approximation  are  compared  to  the  above 
results  In  the  thin  phase-screen  approximation.  For  the  purposes  of  the 
comparison  It  Is  assumed  that  the  phase-screen  Is  located  at  the  center  of 
the  thick  layer,  that  Is 

l,  ■  If  -  I  L  (4-161) 


4.5.3  Comparison  of  Results 

Figures  4-5  and  4-6  show  jthe  thick  layer  results  for  the 
normalized  scintillation  Index  S^//2a  as  a  functio"}  of  the  normalized 
free-space  propagation  distance  x(z^-f)/2nL^  for  values  of  the  normalized 
layer  thickness  ranging  from  0.0001  to  10.  Figure  4-5  Is  a  direct  plot 
of  the  results  of  Equation  4-148  and  Is  valid  for  Isotropic  Irregulari¬ 
ties.  Figure  4-6  Is  a  plot  of  the  results  of  Equation  4-149  and  Is  valid 
for  Infinitely  elongated  Irregularities. 

As  shown  In  these  two  figures  the  value  of  S.//2  o  peaks  at 
unity  for  x(z^-L)/2hL2  large. For  values  of  normalized  layer  thickness, 
smaller  than  the  normalized  free-space  propagation  distance  by  a 
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factor  of  ten  or  more,  the  actual  layer  thickness  is  unimportant  and  the 
result  is  independent  of  layer  thickness  for  very  small  layer  thickness. 
One  must  be  careful  not  to  use  the  results  shown  in  these  two  figures  as 
evidence  that  the  thin  phase-screen  approximation  is  useful  only  for 
situations  where  the  normalized  free-space  propagation  distance  is  ten  or 
more  times  the  normalized  layer  thickness. 

Figures  4-7  and  4-8  show  the  differences  between  the  scintilla¬ 
tion  index  calculation  in  the  Rytov  approximation  and  in  the  thin  phase- 
screen  approximation  where  the  phase-screen  is  located  at  the  center  of 
the  original  thick  layer.  T^^e  abscissa  is  the  normalized  free-space 
propagation  distance,  ordinate  is  the  error  quantity 

S^(Rytov)  -  S|( phase-screen) 

Sj( Rytov) 


plotted  in  percent.  Values  of  the  normalized  layer  thickness,  shown 
here  range  from  0.0001  to  10.  In  all  cases  the  error  is  small  even  when 
the  normalized  frae-space  propagation  distance  and  the  normalized  layer 
thickness  are  identical.  For  a  fixed  normalized  layer  thickness  the  dif¬ 
ference  between  the  two  approximations  increases  with  decreasing  normal¬ 
ized  free-space  propagation  distance  until  a  fixed  maximum  difference  is 
reached.  This  fixed  maximtm  difference  corresponos  to  the  thick  layer 
geometry  where  the  receiver  is  located  at  the  edge  of  the  original  thick 
layer,  that  is  at  z^  =  L.  Note  that  the  value  of  the  maximum  difference 
depends  on  the  normalized  thickness  as  is  shown  further  in  the  next  two 
figures. 

Figures  4-9  and  4-10  show  the  values  of  the  percent  difference 
between  the  Rytov  thick  layer  calculation  of  and  the  thin  phase- 
screen  calculation  of  S4  for  a  thick  layer  geometry  in  which  the 
receiver  is  located  at  the  edge  of  the  scattering  layer,  that  is 
5z  =  Equations  4-148  and  4-149.  In  the  thin  phase-screen 
approximation  the  phase-screen  is  located  at  the  center  of  the  layer  so 
that  5e  "  Equations  4-158  and  4-159.  In  both  the  case  of 

isotropic  irregularities  (Figure  4-9),  and  the  case  of  infinitely 
elongated  irregularities  (Figure  4-10),  the  maximum  difference  is  fifty 
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\(  Zr-L)/2nLl 

Figure  4-7.  Percent  difference  In  between  the  Rytov  and  the  thin 
phase-screen  approxlnatlon.  Isotropic  Irregularities. 
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Figure  4-8.  Percent  difference  In  Su  between  the  Rytov  and  the  thin 
phase-screen  approxlwatlon,  elongated  Irregularities. 
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percent.  This  difference  reflects  the  worst  case  accuracy  of  the  thin 
phase-screen  approximation  in  the  weak  scattering  regime.  Note  that  the 
difference  has  a  minimum  in  both  figures  corresponding  to  points  where  the 
thin  phase-screen  approximation  gives  a  result  identical  to  the  Rytov 
approximation. 

In  Section  3  of  this  report  the  thin  phase-screen  approximation 
is  used  to  determine  the  onset  of  severe  scintillation  defined  as  the 
point  at  which  the  scintillation  index  has  a  value  of  0.5.  For  this 
our pose,  a  maximin  difference  of  fifty  percent,  when  compared  to  the  more 
complicated  Rytov  approximation,  is  considered  acceptable. 

These  brief  results  illustrate  the  following  points  regarding 
the  use  of  the  thin  phase-screen  approximation  for  weak  scattering. 

First,  the  accuracy  of  the  thin  phase-screen  approximation  for  a  thick 
layer  is  dependent  on  the  relative  lengths  of  the  free-space  propagation 
path  and  of  the  propagation  path  through  the  turbulent  ionized  layer. 
Second,  for  a  free-space  propagation  distance  ten  or  more  times  larger 
than  the  layer  thickness,  the  thin  phase-screen  approximation  is  an 
accurate  model  of  a  thick  layer  with  five  percent  or  less  error  for  S^. 
Third,  for  equal  propagation  distances  in  free-space  and  through  the  thick 
layer,  the  error  in  the  thin  phase-screen  approximation  is  small.  Fourth, 
for  a  thin  phase-screen  centrally  located  to  replace  a  thick  layer,  the 
accuracy  of  the  thin  phase-screen  approximation  is  a  function  of  the  nor¬ 
malized  layer  thickness  The  maximun  error  in  is  fifty  percent. 

Therefore  the  thin  phase-screen  approximation  is  a  useful, 
accurate  tool  to  use  to  determine  the  onset  of  severe  Rayleigh  scintilla¬ 
tion  for  the  purposes  of  statistical  signal  generation. 

4.5.4  <T>  and  -  Thick  Layer  Versus  Ihin  Screen 

In  this  section  the  simplified  results  for  <t>,  the  mean  time 
delay,  and  a^,  the  time  delay  jitter,  both  in  the  thin  phase-screen 
approximation,  are  compared  to  exact  results  in  the  strong  scatter 
approximation  obtained  from  Equations  4-117  and  4-118  by  integrating  the 
exact  power  impulse  response  function.  For  the  exact  calculation  the 
power  impulse  response  function  is  obtained  by  numerically  taking  the 
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Fourier  transform  of  the  two-frequency  mutual  coherence  function  as 
specified  by  Equation  4-13.  The  required  two-frequency  mutual  coherence 
function  is  evaluated  numerically  in  the  general  case  of  a  thick 
scattering  layer  from  Equations  4-71  and  4-72  for  isotropic  and  elongated 
irregularities,  respectively.  Then  <t>  and  are  obtained  by 
numerically  integrating  according  to  Equations  4-117  and  4-118. 

In  all  cases  in  this  subsection  the  geometry  under  consideration 
is  shown  in  Figure  3-1.  For  convenience  the  transmitter  distance  is 
taken  as  infinity  so  that  the  incident  wave  is  plane.  The  layer  thickness 
is  L  and  the  free-space  propagation  distance  outside  the  layer  is  z^-L. 

A  K”**  three-dimensional  irregularity  spectrum  is  used  with  a  ratio  of 
outer  scale  to  inner  scale  of  100. 


In  the  thin  phase-screen  approximation  the  important  parameters 
(D^oh  and  a  have  simple  expressions  for  the  geometry  of  interest, 

Vn^  <*-“> 


"coh 

a  «  ■ 


(4-162) 


(4-163) 


Equations  4-162  and  4-163  are  obtained  directly  from  Equations  4-74  and 
4-106.  For  steeper  spectra  the  constants  in  these  expressions  take  on 
different  values  but  the  basic  algebraic  dependence  on  the  important  quan¬ 
tities  xz^/Lq^  and  remains  unchanged.  The  expression  for 
may  be  substituted  into  Equations  4-132  and  4-134  to  obtain  for  the  thin 
phase-screen  case  and  isotropic  irregularities 


\  w 


(4-164) 

(4-165) 


For  elongated  irregularities.  Equations  4-133,  4-134,  and  4-162  yield 

(4-166) 


o2  \2  a2)  2«L2 


(4-167) 


To  model  a  layer  of  thickness  L  with  transmitter  at  infinity  and  receiver 
located  a  distance  from  the  entrance  to  the  layer,  a  thin  phase- 
screen  located  at  the  layer  cente**  has 


Z  rf  “  z  -  L/2 
reff  r 


in  Equations  4-164  to  4-167. 


(4-168) 


Therefore,  based  on  the  thin  phase-screen  results  of  Equations 
4-164  to  4-167,  it  is  convenient  to  plot  the  expressions  on  the  left-hand 
side  of  the  four  equations. 

Figure  4-11  shows  the  normalized  mean  time  delay  for  isotropic 
irregularities 


<Tt>(»)  /o2 

I  0  i 


as  a  function  of  the  normalized  free-space  propagation  distance 
for  values  of  the  normalized  layer  thickness  Cl  >^anging  from  10“®  to 
10.  As  defined  previously 


198 


ISOIROPIC  IRRERULARITIES 
K-"  30  IRRERULARITY  SPECIRIIM 


30  IRRERIILARItY  SPECIRUM 


•W.  £.  ft* 


A»y^  z^.’;.-,-  i'-*'»v-:.V'i<.y;^^^^v--^*»y^'v 


«2  =  «r/2,L2 

5|_  =  XL/21.L2 

Results  applicable  to  the  case  of  elongated  irregularities  may  be  obtained 
by  dividing  the  ordinate  by  two.  It  is  apparent  from  the  figures  that  for 
any  given  free-space  propagation  distance,  the  value  of  the  normalized 
mean  time  delay  converges  to  the  thin  phase-screen  result  as  Cl  de¬ 
creases. 


For  a  constant  value  of  the  normalized  layer  thickness,  5^^,  the 
normalized  mean  time  delay  becomes  a  constant  function  of 
values  of  the  normalized  free-space  propagation  distance  less  than  about 
one-tenth  of  the  normalized  la)«r  thickness.  This  situation  corresponds 
to  a  thick  scattering  layer  with  very  small  free-space  propagation  dis¬ 
tance  so  that  changes  in  the  free-space  propagation  distance  are  unimpor¬ 
tant  and  do  not  affect  the  geometry. 

Figure  4-12  shows  the  normalized  time  delay  jitter 


as  a  function  of  the  normalized  free-space  propagation  distance, 

Values  of  the  normalized  layer  thickness  ranging  from  10"®  to  10  are 
shown.  The  general  appearance  of  this  figure  is  similar  to  that  of  the 
previous  figure.  The  limiting  behavior  described  above  for  decreasing 
layer  thickness  at  a  constant  free-space  propagation  distance  and  for  de¬ 
creasing  free-space  propagation  distance  at  a  constant  layer  thickness 
also  occurs  in  Figure  4-12.  That  is,  for  a  given  normalized  free-space 
propagation  distance,  the  normalized  time  delay  jitter  converges  to  the 
thin  phase-screen  result  as  the  normalized  layer  thickness  decreases.  For 
a  constant  value  of  normalized  layer  thickness,  the  normalized  time  delay 
jitter  reaches  a  constant  value  for  normalized  free-space  propagation  dis¬ 
tance  values  less  than  about  one-tenth  of  the  normalized  layer  thickness. 
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Figure  4-12  is  also  valid  for  the  case  of  infinitely  elongated 
irregularities  where  the  ordinate  in  this  case  is 


At  this  point  consider  the  calculation  of  the  difference  between 
the  exact  thick  layer  results  and  the  thin  phase-screen  results  for  values 
of  the  mean  time  delay  and  time  delay  jitter.  It  was  shown  previously  in 
Section  4.4.1  that  both  calculations  yield  the  s^e  value  of  the  mean  time 
delay.  It  is  also  apparent  that  this  comparison  is  also  available  for  the 
time  delay  jitter  for  a  few  values  of  normalized  layer  thickness  from  the 
curves  of  Figure  4-12.  To  make  such  a  comparison  it  is  important  to  modi¬ 
fy  the  effective  free-space  propagation  distance  shown  on  the  abscissa  so 
that  the  thin  phase-screen  is  located  at  the  center  of  the  original  thick 
scattering  layer. 


Figure  4-13  shows  values  of  the  difference  between  the  exact 
strong  scatter  calculation  and  the  thin  phase-screen  approximation 


201 


a^j(exact)  -  a^j(phase-screen) 
a^j(axact) 

plotted  in  percent  as  a  function  of  normalized  layer  thickness  and  nor¬ 
malized  free-space  propagation  distance.  Figure  4-13  applies  to  the  case 
of  isotropic  irregularities  and  is  restricted  to  a  value  of  a  of  100. 
Values  of  the  normalized  layer  thickness  5|_  range  from  10"**  to  10.  It 
can  be  seen  that  for  values  of  the  normalized  layer  thickness  less  than  or 
equal  to  the  normalized  free-space  propagation  distance,  the  difference  is 
at  most  about  10  percent.  The  maximum  difference  is  independent  of  layer 
thickness  and  is  22.5  percent.  This  maximun  difference  corresponds  to  the 
case  of  a  receiver  located  at  the  exit  of  the  thick  layer  with  zero  free- 
space  propagation  distance.  As  in  the  preceding  calculation  of  the  thin 
phase-screen  approximation  to  the  scintillation  index,  this  geometry  is 
the  most  difficult  to  model  with  a  centrally  located  phase-screen. 

Figure  4-14  shows  another  comparison  of  the  accuracy  of  the  time 
delay  jitter  calculation  in  the  thin  phase-screen  approximation.  The  cal¬ 
culations  presented  in  this  figure  differ  from  the  previous  calculations 
only  in  that  the  value  of  the  parameter  a  is  unity  for  Figure  4-14.  It  is 
apparent  that  the  thin  phase-screen  approximation  is  more  accurate  for 
this  decreased  value  of  the  parameter  a  with  a  maximum  difference  of  about 
9.5  percent  for  unity  a. 

Figure  4-15  shows  the  difference  between  the  exact  strong  scat¬ 
ter  results  and  the  thin  phase-screen  results  for  the  time  delay  jitter 
for  elongated  irregularities  and  unity  a.  The  general  behavior  noted  in 
the  two  previous  figures  is  again  observed  here.  The  maximum  time  delay 
jitter  difference  is  now  decreased  to  about  6  percent  for  this  elongated 
irregularity  case. 

The  previous  three  figures  show  that  the  maximum  difference  in 
the  time  delay  jitter  calculation,  caused  by  modeling  a  thick  scattering 
layer  as  a  thin  phase-screen  located  at  the  center  of  the  layer,  is  a 
function  of  the  parameter  a. 
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Figure  4-14.  Relative  time  delay  jitter  accuracy  of  the  thin  phase-screen 
approximation,  isotropic  irregularities,  o  *  1. 


ELON^TED  IRREGULARITIES 


r^25 

HO 

Q) 

O 

b  20 

CD 

A 


15 

I 

I 

q  5 


..\f*  -  lA 

p- V2  10 


»V  IW  IV  lU  lU  IQ 

Mzr-L)/2nLl 

Figure  4-15.  Relative  time  delay  jitter  accuracy  of  the  thin  phase-screen 
approximation,  elongated  irregularities,  a  =  1. 


Figures  4-16  and  4-17  show  the  maximum  error 


a^(thick  layer)  -  o^(phase-screen) 

(j^(thick  layer) 

in  percent  for  the  isotropic  and  elongated  irregularity  cases,  respective¬ 
ly.  These  curves  were  obtained  by  taking  a  value  of  the  normalized  layer 
thickness,  ^l,  of  10  and  a  very  small  value  of  the  free-space  propaga¬ 
tion  distance,  52-5l,  of  10“®.  This  case  yields  the  desired  maximum 
difference.  It  can  be  seen  in  both  figures  that  the  error  decreases  to 
zero  for  values  of  a  less  than  unity  and  reaches  the  maximun  value  of  22.5 
percent  for  values  of  a  of  about  10  and  greater.  It  was  noted  in  Section 
4.3.1  that  a  is  proportional  to  the  quantity  z^/F  where  z^  is  the 
propagation  distance  and  F  is  the  focal  distance.  Hence  for  a  small,  the 
receiver  is  "close  to"  the  phase-screen  and  dispersive  effects  dominate 
diffractive  effects.  In  this  case  the  time  delay  fluctuations  are  caused 
by  large-scale  phase  fluctuations  and  are  well  modeled  by  a  thin  phase- 
screen.  For  a  large,  diffraction  is  more  important  than  dispersion  and  a 
single  phase-screen  model  of  a  thick  layer  is  relatively  less  accurate. 

It  should  be  emphasized  that  the  maximum  error  in  caused  by 
modeling  a  thick  layer  by  a  thin  phase-screen  is  22.5  percent.  This  error 
appears  to  be  modest  relative  to  uncertainties  in  the  knowledge  of  the 
propagation  environment  in  the  case  of  strong  scattering. 

As  a  final  comparison  of  the  accuracy  of  the  thin  phase-screen 
approximation  to  calculate  the  time  domain  properties  of  a  wide  bandwidth 
waveform  subject  to  strong  scattering,  consider  the  effects  of  the  approx¬ 
imation  on  the  power  impulse'  response  function  Gi(zr,T).  Figures  4-18 
to  4-20  show  the  power  impulse  response  function  for  isotropic  irregular¬ 
ities  for  a  thick  layer  and  its  thin  phase-screen  equivalent.  The 
abscissa  and  ordinate  are  normalized  by  multiplying  and  dividing  by 
ojcoh*  three  figures  the  results  correspond  to  a  value  of  a  of 

3  where  scattering  effects  slightly  dominate  dispersive  effects.  In  all 
three  figures  the  normalized  free-space  propagation  distance,  5z“Cl» 
is  unity;  the  normalized  layer  thickness,  ?[_,  takes  on  values  of  0.1, 

1.0,  and  10  in  the  three  figures.  The  thin  phase-screen  equivalents  thus 


?.n&- 


have  zero  layer  thickness  and  values  of  the  effective  normalized  propaga¬ 
tion  distance  between  phase-screen  and  receiver  of  1.005,  1.5,  and  6 
respectively.  These  values  are  computed  directly  from  Equation  4-168. 

In  Figure  4-18  the  layer  thickness  is  small  relative  to  the 
free-space  propagation  distance  and  the  impulse  response  function  is  exact 
for  the  thin  phase-screen  approximation  to  this  layer. 

In  Figure  4-19  the  layer  thickness  and  the  free-space  propaga¬ 
tion  distance  are  identical.  The  power  impulse  response  function  for  the 
thick  layer  is  quite  similar  to  that  of  the  thin  phase-screen  approxima¬ 
tion  in  this  case  with  no  difference  in  the  mean  time  delay  and  a  9.6 
percent  difference  in  the  time  delay  jitter. 

In  Figure  4-20  the  layer  thickness  is  ten  times  the  free-space 
propagation  distance  with  a  resulting  19  percent  error  in  the  result  for 
0^  obtained  in  the  thin  phase-screen  approximation. 


Figure  4-18,  Comparison  of  power  impulse  response  function,  thin 
phase-screen  versus  thick  layer. 


Results  of  this  subsection  show  that  the  thin  phase-screen 
approximation  gives  an  adequate  approximation  to  the  exact,  strong  scatter 
result  for  time  delay  jitter.  It  has  also  been  shown  that  the  thin-phase 
screen  approximation  provides  an  accurate  approximation  to  the  power 
impulse  response  function  and  thus  may  be  used  as  the  basis  for  the  sta¬ 
tistical  generation  of  signals  as  discussed  in  Section  3  of  this  report. 
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APPENDIX  A 


A.l  EVALUATION  OF  THE  GENERALIZED  POWER  SPECTRUM  S(K,t) 
A. 1.1  Numerical  Method 


In  this  appendix  the  generalized  power  spectrum,  as  calculated 
in  Section  4  in  the  thin  phase-screen  approximation,  is  evaluated  for  a 
number  of  interesting  conditions.  The  evaluation  of  S(K,t)  is  somewhat 
tedious  as  is  seen  here  but  is,  of  course,  quite  necessary  for  the  calcu¬ 
lation  of  statistical  realizations  of  the  impulse  response  function.  In 
the  thin  phase-screen  approximation  the  generalized  power  spectrum  is 
given  by  Equations  4-97  and  4-99  repeated  here  for  convenience 


SJK.t)  ■  exp  •  1 

^  23/2,  4  “ 


t 

i  «3 


(A-1) 


3  ^  “coh^o 


3/2 

^  4nV2 


exp  I  -  —  -  i.  Fi.  K^4^-u)  .  t1 

I  4  0  2  0  coh  J 


X  exp  (3V4o2)Ki/i^(02/4o2) 


where 


B  =  a2  j"!  + 


(A-2) 


(A-3) 


Equation  A-1  gives  the  generalized  power  spectrum  for  elongated 
irregularities  and  is  appropriate  for  the  two-dimensional  geometry  solved 
by  the  MPS  propagation  simulation.  This  equation  is  easy  to  evaluate 
numerically. 
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Equation  A-2  gives  the  generalized  power  spectrum  for  the  case 
of  isotropic  irregularities  and  is  somewhat  more  difficult  to  evaluate. 

It  is  convenient  to  consider  four  cases  and  use  different,  efficient  eval¬ 
uation  algorithms  for  each  case. 

Case  1:  0  <  e2/4a2  <  3.5 

For  32/40,2  between  zero  and  3.5,  Sj(K,t)  is  easily  evaluated 
since  in  this  range  the  Kj/j,  Bessel  function  may  be  directly  cal¬ 
culated  using  the  expressions 


K  (z)  = 
v'  ' 


/I  (z)-I  (z) 

1 1  V 

2  \  sinviT 


(A-4) 


I  (z)  =  (z/2)^  I 

^  k=0 


(2^4)^ 
k!r(  vfk+1) 


(A-5) 


which  are  taken  directly  from  Abramowitz  and  Stegun,  1964.  For  z  £  3.5 
sufficient  accuracy  is  obtained  using  at  most  twelve  terms  in  the 
convergent  series  given  by  A-4.  The  r  function  required  in  Equation  A-5 
is  available  on  most  computer  systems  and  is  also  given  in  Abramowitz  and 
Stegun. 

Case  2:  32/4a2  >  3.5 

For  b2/4o2  >  3.5  it  is  necessary  to  use  the  asynptotic 
expansion  for  the  factor 

exp  (bVAo^)  Ki/4  (bVAo^) 
which  has  the  form 
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where  p  =  4v^.  Equation  A-6  is  taken  directly  from  9.7.2  in  Abramowitz 
and  Stegun. 

Cases  3  and  4:  $  less  than  zero 

For  3  negative  it  is  necessary  to  carefully  evaluate  the  term 
S*/"  K,/,(62/4„2) 

This  is  easily  done  with  the  aid  of  Equation  9.6.31  from  Abramowitz  and 
Stegun.  The  result  is  found  as 

3^/^  Ki/4  = 

l3l^/^{Kl/4(eV4a2)  +  n/2  Ii/4(3V4a2)}  .  3<0  (A-7) 


Also  one  needs  the  identity 


(A-8) 


Substituting  Equation  A-7  into  the  expression  for  Sj(K,t)  and  replacing 
the  appropriate  factors  in  the  exponent  by  use  of  A-8,  one  obtains  for 
3  <  0 


SjCK.t) 


“coh'olsl 

4.V4 


exp  -  i  k2£2 


X  jexp  |- 1  ^  I  exp  (  3V4a^)Ki/4(  3^/4o^ 

+  2^2^  exp  (^  '  2 )  (-B^/4a^)Il/4(3^/4a2)|  (A-9) 
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Case  3:  e  <  0,  0  <  32/45,2  ^  3.5 

For  3  negative  and  32/4o2£  3.5  accurate  results  can  be  obtained 
using  the  convergent  series  for  Kj/^  and  Ij/^  given  by  Equations  A-4  and 
A-5  to  evaluate  the  expression  for  the  generalized  power  spectrum  given  by 

A-9. 

Case  4:  3  <  0,  3^40^  >  3.5 

For  3  negative  and  3^/4 greater  than  3.5  use  of  the  convergent 
series  expression  becomes  inefficient  and  causes  inaccurate  numerical 
results.  Instead  use  the  asymptotic  expansion  given  by  Equation  A-6  for 
the  term 


exp  (3^/40^)  K1/4  (32/4a2) 


to  evaluate  the  term  with  the  K  Bessel  function.  For  the  remaining  term, 
exp  (-  3^40^)  11/4(3^40^) 

use  the  asymptotic  expression  9.7.1  in  Abramowitz  and  Stegun. 


e"^I  (z)  =»4=!i  -  id.  -t-  ( 

''  /2Ttz  (  8z  2!(8z)2 

-  (u-l)(v-9)(y-25)  ^  ) 

3!(8z)3  ”  ) 


(A-10) 


where  y=4v2. 


The  numerical  techniques  described  in  this  section  enable 
accurate  and  efficient  evali.c'tion  of  the  generalized  power  spectra.  In 
the  following  subsection  results  for  these  functions  are  shown  for  a 
number  of  different  values  of  the  parameter  a. 
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A.1.2  Examples  of  S(K,t} 


In  this  section  plots  of  the  generalized  power  spectra  S£(K,t) 
and  Sj(K,t)  are  shown  as  a  function  of  the  parameter  a.  This  parameter 
is  described  in  Section  4  where  it  is  seen  that  a  much  less  than  unity 
corresponds  to  the  case  v4iere  large  scale  phase  effects  and  dispersion  are 
more  important  than  angular  scattering  or  diffraction,  a  greater  than 
unity  corresponds  to  the  case  where  diffraction  provides  the  dominate 
scattering  mechanism. 

Figures  A-1  to  A-20  show  plots  of  the  function  S£(K,t)  and 
Sj(K,t)  for  values  of  the  parameter  o  of  10.0,  5.0,  1.0,  0.1,  0.01.  In 
all  cases  the  generalized  power  spectra  are  plotted  as  a  function  of  the 
normalized  variables  AqK  and  WQoh'^*  ^  convenient  normalization 

as  can  be  readily  seen  from  Equations  A-1  and  A-2.  Two  types  of  figures 
are  shown  here.  In  one  case  contours  of  the  generalized  power  spectra  are 
shown  as  functions  of  the  two  normalized  variables  and  w^obT*  A 
second  plot  of  S(K,t)  is  provided  which  gives  a  three-dimensional  view  of 
this  function  to  provide  better  understanding  of  the  functional  depend¬ 
ences. 

Figures  A-1  to  A-4  for  a  value  of  a  of  10  shows  the  behavior  of 
the  generalized  power  spectra  for  large  o.  The  delta  function  behavior 
for  the  case  of  elongated  irregularities  as  given  by  Equation  4-110  is 
shown  in  Figures  A-1  and  A-3.  Recalling  from  Fourier  analysis  that  the 
variable  K  is  related  to  the  sine  of  the  scattering  angle,  e  in  the  para¬ 
bolic  approximation,  for  small  e 

K .  M  (A-ii) 

X  X 

Since  the  generalized  power  spectrin  takes  on  the  form 
6  (I  K2to^-«coh^| 
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when  diffractive  or  angular  scattering  effects  are  dominant,  there  is  a 
one-to-one  correspondence  between  scattering  angle  and  time  delay.  That 
is. 


K  =  +  —  “  —  ■—  ♦  a»l  (A-12) 

£o  X 

Of  course,  the  resulting  time  domain  power  impulse  response 
function  is  the  integral  of  the  generalized  power  spectrun  over  all  K  or 
equivalently,  over  all  scattering  angles  (e.g..  Equation  4-107). 

The  effect  of  the  difference  between  isotropic  irregularities 
and  infinitely  elongated  irregularities  is  also  seen  from  a  comparison  of 
Figures  A-1  and  A-3  versus  A-2  and  A-4.  Notice  that  for  the  case  of  iso¬ 
tropic  irregularities,  the  characteristic  horseshoe  shape  given  by  the 
equation 


“coh^4 

has  a  tendency  to  be  filled  in.  For  elongated  irregularities,  less  energy 
is  present  inside  the  horseshoe. 

The  reason  for  this  difference  is  that  angular  scattering  can 
occur  in  only  one  plane  for  the  case  of  infinitely  elongated  irregular¬ 
ities.  For  large  o  and  elongated  irregularities  a  single  time  delay  can 
result  from  scattering  at  only  two  angles  as  specified  by  Equation  A-12. 
For  isotropic  irregularities,  angular  scattering  occurs  in  two  planes  and 
a  large  number  of  different  scattering  angles  can  yield  the  same  time 
delay  even  for  large  o.  It  is  the  angular  scattering  from  the  additional 
plane  which  causes  the  horseshoe  shaped  figure  to  be  filled  in  for  the 
case  of  isotropic  irregularities. 

Figures  A-5  through  A-16  show  the  development  of  the  generalized 
power  spectra  as  the  parameter  a  decreases  from  5.0  to  0.1. 
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Figures  A-17  to  A-20  for  the  case  of  a  of  0.01  show  the  gener¬ 
alized  power  spectra  for  the  case  where  angular  scattering  is  unimportant 
and  large  phase  effects  and  dispersion  are  the  dominant  factor;  As  dis¬ 
cussed  earlier,  the  case  of  small  a  corresponds  to  the  basic  geometry 
where  the  focal  length  of  the  irregularities  is  much  larger  than  the  prop¬ 
agation  distance  to  the  receiver.  Equivalently  the  receiver  is  located 
very  close  to  the  scattering  medium.  In  this  case,  the  phase  irregulari¬ 
ties  comprising  the  scattering  layer  appear  as  very  large,  frequency 
selective  lenses  which  cause  little  angular  scattering  but  which  may  cause 
different  frequency  components  of  a  wide  bandwidth  waveform  to  propagate 
at  different  velocities.  These  large  phase  excursions  in  the  medium  may 
also  cause  large  scale  changes  in  the  mean  time  delay.  In  either  case  the 
generalized  power  spectrum  for  isotropic  or  elongated  irregularities  has 
the  Gaussian  form  given  by  Equation  4-112  and  shown  in  Figures  A-17  to 
A-20. 
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Generalized  power  spectrum  for  isotropic  irregularities 


coh 


h 


3 


6 


4 


2 


0 


CONTOURS 
A:  1.0E+00 
B:  5.0E-01 
C:  3.0E-01 
0:  1.0E-01 

E:  1.0E~02 
F:  1.0E-03 
G:  1.0E-04 


0.0 


J _ I 

3.5  7.0 


Figure  A-6.  Generalized  power  spectrum  for  isotropic  irregularities, 
a  =  5. 
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Generalized  power  spectrum  for  elongated  irregularities 
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